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INTERIOR  BALLISTICS  OF  FORCEMENT  PERIOD 


EVGENY  V.  CHURBANOV 
Doctor  of  Technical  Sciencies, 
Professor, 

State  Baltic  Technical  University, 
St.  Petersburg 


The  onward  movement  of  projectile  inside  the  barrel  takes  place  due  to 
pressure  force  coursed  by  powder  gas  towards  the  bottom  of  the  projectile,  and 
the  drag  forces  caused  by  interconnection  between  the  projectile  and  the  barrel, 
i.e.  drag  force  caused  by  gripping  of  projectile  by  rifling,  by  driving  band 
sliding  along  the  rifles,  by  resistance  that  appears  in  course  of  interaction 
between  positioning  band  and  the  barrel, 

A  special  ballistic  parameter  -  engraved  pressure  -  is  used  in  interior 
ballistics  for  indirect  consideration  of  the  process  of  gripping  of  projectile  by 
rifling.  In  this  case  the  gripping  drag  force  is  not  considered  in  die  projectile 
onward  movement  equation.  At  the  same  time  at  the  beginning  period  of  the 
projectile  movement,  which  is  called  the  period  of  forcement,  when  driving 
band  grips  into  rifles,  the  drag  force  is  comparable  to  powder  gas  pressure  force. 
The  period  of  forcement  is  the  most  complicated  stage  of  the  shot  due  to  the  fact 
that  besides  the  main  processes,  i.e.  powder  burning,  powder  gas  formation, 
projectile  onward  movement,  gripping  of  projectiles  by  rifles,  additional 
processes  are  to  be  also  considered,  i.e.  ignition  of  explosive  propellant,  powder 
gas  heat  emission,  outflowing  of  powder  gas. 

The  process  of  gripping  of  projectile  by  rifles  to  a  large  extent  depends 
on  the  gun  loading  method.  In  case  of  separate  loading,  the  gripping  process 
starts  at  the  beginning  of  the  onward  movement  of  projectile,  and  the  process 
itself  is  relatively  slow.  In  case  of  unitary  loading  the  projectile  gets  out  of  the 
projectile  case  before  the  gripping  process  begins,  so  the  gripping  process  starts 
under  dynamic  conditions,  and  it  is  fast. 
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In  the  process  of  gripping  there  appears  contact  surface  between  driving 
band  and  inner  wall  of  the  barrel.  This  contact  surface  is  characterized  by 
development  of  normal  tension  and  tangent  tension.  As  the  result  the  gripping 
resistant  force  becomes  a  factor.  This  force  can  be  determined  by  experiment 
and  by  calculation. 

Without  simplification  the  main  task  of  forcement  period  interior 
ballistics  can  only  be  solved  numerically.  One  will  need  to  integrate  the  system 
of  equations  that  includes  9  simple  differential  equations  and  18  algebraic 
equations  containing  28  variables  and  more  than  100  ballistic  parameters.  For 
engineering  purpose  the  system  of  equations  can  be  simplified,  and  analytic 
dependencies  of  different  phases  of  projectile  onward  movement  on  time  v/ithin 
the  period  of  forcement  can  be  determined. 

The  period  of  forcement  can  be  considered  indirectly  with  accuracy 
needed  with  the  use  of  three  forcement  parameters,  not  one  parameter  as  is 
usually  done.  Investigation  of  those  parameters  leads  to  understanding  the 
physical  essence  of  engraved  pressure,  and  helps  in  finding  methods  of  its 
theoretical  and  experimental  determination. 

In  the  author’s  book  “Interior  Ballistics  of  Forcement  Period”  issued  in 
1997  in  State  Baltic  Technical  University  one  can  find  detailed  description  of 
the  theory,  experiments  and  methods  of  consideration  of  period  of  forcement  in 
solving  ballistic  problems. 
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ESTIMATION  OF  INFLUENCE  INPUT  PARAMETERS  ON 
SERVICEABILITY  OF  BALLISTIC  SYSTEMS 


ZAKHARENKOV  V.F.,  MOROZ  S.  A.,  Baltic  State  Technical  University, 
198005,  Saint-Petersburg,  Russia 

The  estimation  of  serviceability  of  ballistics  launchers  demands  to  take 
into  account  the  influence  of  input  random  information  at  output  characteristics 
such  as  maximal  or  muzzle  gas  pressure,  the  muzzle  projectile  velocity,  position 
of  the  powder  burning  end  and  so  on.  All  that  demand  the  modernization  of 
existing  methods  of  ballistic  calculations  and  design  of  sets  in  the  direction  of 
researching  the  random  relationships,  escorting  the  shot  phenomena.  Random 
spread  of  output  parameters  (the  response  functions)  is  wholly  and  fully 
determines  by  relationships  of  input  parameters  spreads. 

The  response  functions  can  be  install  as  theoretically  so  experimentally. 
The  using  of  the  models  with  that  or  other  accuracy  degree  describing  real 
ballistic  processes  are  demand  at  theoretical  solution.  The  connection  between 
probability  values  of  input  and  output  data  may  be  ascertain  or  on  the  basis  of 
extreming  deviation  method,  or  by  the  Monte-Carlo  statistical  tests'  method. 
The  first  method  is  a  more  simple,  but  he  does  not  reflect  actual  relationships  of 
random  deviations  of  output  parameters,  and  mostly  base  on  their  apriori 
assumption  on  the  normal  (Gauss)  law  of  distribution.  The  more  reliability 
results  give  the  Monte-Carlo  statistical  test  method  that  has  not  limitations  on 

complexity  systems  but 
demand's  availability  of  great 
statistical  sampling. 

The  module  kla.exe  of 
the  block  ballistic  calculations 
KOBRA  of  Baltic  State 
Technical  University  (BSTU) 
was  chose  as  a  working 
program  for  ballistic 
calculations  in  the  frames  of 
present  paper.  These  program 
modules  were  modify  in 
direction  of  inclusion  in  it  the 
procedures  for  forming 
random  pseudonumbers 

generator  RAV  and  NORM  for  modeling  random  input  parameters  distributed 
on  uniform  Gauss  and  laws  (Table  1).  The  choice  of  distribution  laws  was  base 
on  analysis  the  statistical  information  of  a  number  ballistic  systems.  The 
extreme  deviations  of  parameters  (constructive,  physical  and  of  load 
conditions)  chosen  out  from  condition  ^=±5%  from  nominal  values  and 
corresponded  to  3  a -rule  for  the  level  of  confidence  probability  as  0,  997.  The 


Table  1 


Characteristics 

Symbol 

Distribution  law 

Burning  rate 

u, 

Gauss 

Covolume 

a 

Uniform 

Charge  density 

A 

Gauss 

Energy  of  powder 

f 

Gauss 

Initial  pressure 

Po 

Uniform 

Ignition  pressure 

PB 

Uniform 

Powder  density 

6 

Uniform 

Projectile  mass 

q 

Gauss 

Charge  mass 

a> 

Gauss 

Initial  chamber  volume 

Wo 

Gauss 

Bore  cross-section  area 

S 

Gauss 

Flame  temperature 

Ti 

Uniform 

Caliber 

d 

Gauss 
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integral  function  of  distribution  F(x)  calculated  as:  P-ck— (x-a)  .  for 

Jb-a  b-a 

a 

X 

uniform  distribution  law;  F(x)=  J(p(z)dz=0f+0icD(x).  for  nonnal  (Gauss) 


distribution  law,  where  ‘l(x)= 


The  calculation  procedure  for  random-numbers  generation  for  making  input 
random  information  given  on  Fig.l.  The  system  of  serviceability  conditions 

may  be  write  as:  N>ri;  Ps<[Ps]; 

MP3];  d^<[dj. 

The  restrictions  on  values  of  parameters 
must  be  give.  The  probability  estimation 
demands  the  knowledge  distribution  of 
reliability  density  f(n),  which  for  scheme 
with  accumulation  of  faults  (for  example, 
decreasing  ballistics  due  to  erosion  wear) 
have  view  [5]_^f(n)=ir(^^r'^^n^‘e>g<-An),  if 
n>  1,  where  r(N)  -  gamma-function.  If  the 
level  of  the  aimfunction  is  knov/  that 
probability  P(N>n)  may  be  calculates  as 
P(N>n)=l-F(n).  The  calculation  of  the 
total  probability  with  allowance  of  other 
serviceability  conditions  may  be  carry  out 
in  compliance  with  relation,  which  takes 
into  account  interaction  effects  and  three-parametrical  conditions  of 


Input  data 


Procedure  for  Gauss 

Procedure  for 

raiidorn-numbers 

Unifonu  random- 

generation 

numbers  generation 

^  . -  -  ▼ 

Calculation  procedure  for  forming 
_ random  input  data 


Procedure  for  calculation  of  random 
values  of  ballistic  parameters 

Calculation  procedure  for  estimation  of  randoiW 
_ output  values  of  ballistic  parameters  ] 


Fig.l  Block  -scheme  for  forming  of  random 
values  of  ballistic  parameters 


serviceability  [5,6]:  ^n)  =  n(P,„.  )j;  (P„„,)j = 

where: 


TO  m  n 


(Py.yry,)]-  probability  of  serviceability  conditions  for  j-loading;  m-  number 
of  intervals;  i-  elementary  interval  of  serviceability;  y,;-  probability  Yi  (V- 
velocity)  realizing  inside  i-intep/al;  -  conditional  probability  Y2  (Pm),  realizing 

inside  i-interval;  Py3i/Y,y,j-  condition  of  probability  Ys  (P^  ),  realizing  inside  i 
and  j -intervals. 

The  calculation  of  distribution  character  parameters  characterizing 
serviceability  and  stability  of  ballistic  set  ( gas  pressure  and  projectile  velocities), 
in  characteristic  time  moments  corresponding  to  achievement  of  initial  pressure 
(to),  maximal  pressure  (tm),  finish  time  of  powder  burning  (Ik)  and  muzzle  time 
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(til)  at  the  example  of  repeated  many  times  random  shot-loading  of  122  mm 
caliber  gun  system  reflected  on  Fig.  2. 


Fig. 2  Histograms  of  pressure  and  projectile  velocity's  distributions  in 
characteristic  moments  of  time 

Given  gi'aphs  allow  to  speak  about  the  character  of  variation  not  only  of 
output  values  of  ballistic  parameters  v  and  p,  but  about  their  dynamics  of 
changing  in  different  moments  of  the  shot  period.  The  treatment  of  got  results 
have  let  to  build  the  variation  field  of  velocity  and  pressure.  The  results  of 
calculation  represented  on  Fig,  3. 


Fig.  3  Variation  fields  of  gas  pressure  and  projectile  velocities 


The  values  of  variation  fields  of  parameters,  calculated  on  extreme  deviations  of 
input  data,  given  on  the  same  figure.  One  can  see  that  for  100  and  more  trials 
the  deviation  fields  on  both  methods  are  slightly  differ.  So  approximate 
estimation  of  variation  parameters,  according  to  the  extreme  deviation  method, 
may  be  use  for  preliminary  analysis  of  the  opportunities  of  ballistic  systems  with 
sufficient  accuracy  for  practice. 

For  replying  to  question  about  statistical  significance  of  got  results  and 
the  choice  function,  adequately  describing  Monte-Carlo  method  results,  it  is 
necessary  to  make  estimation  of  coefficients  in  the  law  describing  probability 
distribution  of  random  values  output  parameters  on  limited  sample  from 

x<|3 

general  population  P(a<X<p)  =  ^P, .  For  distribution,  for  level  of 

X>a 

confidence  probability  a=QQ5,  the  average  significance  ^  =  56.3  MPa  with 
value  of  confidence  dispersion  Sx  =0.09746  MPa  (  for  n=300  )  was  gotten.  That 


for  the  execution  of  condition  p|t„„  <^<t,_„„|  =  l-a  gj^es  for  probability 

P-0,  975  crictical  significance  973  =2,  26.  Then  confidence  interval  for  muzzle 
pressure  will  compose  56.05  56.50.  Analogically  for  muzzle  projectile 

velocity  it  will  have  gotten:  680.3 <  <682.3  m/s. 

Practical  application  of  got  distributions  may  be  use  for  solution  of 

whole  row  tasks:  the  calculation  of 
xn.s  dispersion  or  hit  probability  in  a  target, 

'  serviceability  -estimation  of  recoil 
mechanisms  and  other  aggregates  of 
guns,  bore  stress  and  the  dynamics  of  it 
vibrations  and  so  on.  As  example  of 
using  got  results  on  Fig.4  are  given 
epures  of  tangential  stresses  in  4  cross- 
calculated  according  to  formula 

P,  =0,  =  0  and  the  random 


Fig.  4.  The  distribution  of  tangential 
stresses  in  different  bore  cross-sections 


sections 
2 


of 

„2 


bore  wall, 
^  rf2r?+r^ 


E  2pb'2r'+r'  1 

3  r'  r'-rf  3  r'  r'-r'  3  ^ 


at 


assigm-nent  of  gas  pressure  on  walls  inside  a  bore. 
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BURNING  OF  7  AND  19-HALLS  POWDER  GRAINS  WITH 
TAKING  INTO  ACCOUNT  REAL  GEOMETRY  OF  SMALL 
PRIZMS  AFTER  COLLAPSE 

V.  F.ZAKHARENKOV,"^2i\X\c  State  Technical  University,  198005, 
Russia,  Saint-Petersburg 


The  calculation  of  burning  parameters  and  the  geometry  of  powder  grain 
with  n  halls  until  a  moment  of  collapse  does  not  present  difficulties  and  may  be 
carried  out  according  to  M.E.  Serebriakov  or  J.  Corner  methods  [  I  ],  [  2  ]. 
Geometrical  characteristics  of  powder  grain  before  the  start  of  burning  and  at 
the  moment  of  collapse  (when  halls  meet  each  other)  are  shown  on  Fig.l. 


Fig.l  Geometry  of  7>and  19-halls  powder  grain 

The  calculation  of  post-collapse  parameters  requires  the  accounting  of 
real  geometrical  shape  of  collapse  elements  -  curvilinear  small  prizms.  The 
complex  form  of  such  prizms  demands,  for  calculation  of  their  burning  surface, 
application  of  mathematical  apparatus  of  analytical  geometry  and  computer 
technique  for  providing  of  necessary  calculations.  The  existing  methods  [1],  [2] 
were  created  during  a  time  when  there  were  no  computer  technique  and 
demanded  the  elaborating  of  simple  calculation  methods.  As  a  result,  two 
schemes  of  such  calculations  have  been  developed.  In  accordance  with  the  first 
scheme  the  relative  thickness  of  powder  grain  at  the  moment  of  collapse  was 
accepted  equal  to  unit  as  to  the  second  one  -  less  then  a  unit.  In  the  practice  of 
powder  grain  parameters  calculation  in  Russia  the  first  scheme  (method  prof. 
M.E.Serebriakov  [1])  found  it's  application.  The  second  scheme  is  widly  used 
in  the  countries  of  West  Europe,  USA,  China  and  some  other  countries.  In 
Russia  the  second  scheme  was  developed  by  professor  prof.V.E.Slukhotsky  [  1  ]. 

In  accordance  with  the  first  scheme  under  z,=  1  the  bumning  outside 
and  inside  surfaces  are  connected  and  collapse  of  grains  on  several  groups  of 
curvelinear  prizms  takes  place:  for  7-hall  grains  -  2  groups;  for  19-hall  grain  -3 
groups  (Fig.2).  After  collapse  of  grain  a  local  groups  of  prizms  burn 
indepedently. 


4^ 

0 

ripHSMa  N  I 

npH3Ma  N  2 

npH3MaN  1 

ripMSMaN  2 

npH3Ma  N  3 

Fig. 2  Shape  of  prizms  of  7  and  19-halls  powder  grains  at  the  collapse  moment 


The  calculation  of  surface  burning  grain  during  this  period  may  be  given 
on  relative  burning  thickness  for  whole  groups  of  prizms  or  as  to  real  geometry. 
Classical  ballistics  is  based  on  the  first  approach,  theoretical  bases  of  which  in 
Russia  have  been  developed  by  professors  G.  V.  Oppokov,  M.E.  Serebriakov, 
K,  K.  Greten  and  I.P,  Grave  [1],  The  essence  of  the  given  approach  consists  in 
the  calculation  of  relative  burning  thickness  Zj  which  change  in  a  range  from  1 


to  z^  at  alteration  \|;  from  to  1.  By  this  the  value  z^ 


,  where  p,  -  the 


radius  of  circle  inscribing  into  outside  prizra  of  7-hall  powder  grain  and 
calculated  as  to  formula  I.P.  Grave  p,  =0,l772(d+2e|).  On  the  basis  of  geometrical 
approach  (Fig. 3)  the  following  formulas  for  calculation  of  radiuses  of  prizms 

inscribing  circles  were  got  by  us 


Pi  = 


P2  = 


f  3-S 


U-V3. 

(4"' 


(d  +  2ei); 


(d  +  2e,)  .  In  19-hall  grain 


the  third  group  of  prizms  appears  at  the 
moment  of  collapse.  The  radius  of 
circle  pj ,  inscribing  into  the  section  of 
these  prizms  is  equal  to 

25-40cosl5“ +  16cos']5Vj  ^  , 

p,  = - (d  +  2e,)  . 

'  ]2-8cosl5‘’ 

The  calculation  results  of 
geometrical  characteristics  of  standard 
(d  =  e|),  (2.^  =  2,5D)  7  and  19-halls 
powder  grain  and  burning's  parameters  before  and  after  collapse  for  both 
approaches  are  given  in  Table  1 . 


Fig.3  The  calculation  scheme  of 
the  grain  geometry 
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TaSjiHEta  1 


napaMCTp 

Moaejib  1 

MoflCJib  2 

4/7 

11/19 

4/7 

11/19 

ll.MM 

0,4 

1.16 

0,4 

1.16 

d3,  MM 

0,2 

5,8 

0,2 

5,8 

D,  MM 

2,2 

9,86 

2,2 

9,86 

1 

5,5 

24.65 

.5,5 

24.65 

HCl 

^ssam 

0,5804 

0,7343 

0,5804 

0,4162 

0,2337 

0,4162 

-0,02174 

-0,02162 

-0,02174 

-0,02162 

Vs 

0,8536 

0,8094 

0,8536 

0,8094 

Zk 

1,5316 

0,6529 

0,6529 

-0,8436 

-1,0982 

12 

-0,9406 

1,4406 

1,4406 

CTs 

1,4022 

1,7676 

1,7676 

For  calculation  of  current  burning  surface  of  prizms  under  the 
assumption  of  layers  parallel  burning,  i.e.  as  to  geometrical  law,  it  is  necessary 
to  know  the  thickness  of  burned  layer  in  every  moment  of  time.  This 
thickness  may  be  computated  on  the  basis  of  power  law  of  burning  u  =  UlP^  As 

u,  .  then  ei,  =  fuip''dt.  From  this  follows  that,  by  grains  burning  at 

at  dt  '  J 

! 

the  same  pressure  and  physical-chemical  and  kinetical  characteristics  of  powder 
material,  the  value  ei  for  all  prizms  will  be  alike  and  it  may  be  put  =e,  . 
Under  such  conditions  prizm  with  radius  p,  will  over  burn  the  first,  then  -  with 
radius  p-  and  the  last  -  with  radius  Pi.  Thus,  during  the  calculation  of  burning 
parameters  of  all  prizms  the  identical  thickness  step  value  dp  =  p, /n  may  be  set. 
The  last  one  is  determined  as  to  smallest  radius  p, ,  where  n  -  the  number  of 
radius  divisions  on  elementary  thicknesses.  Then  current  thickness  is  e  =  idp , 
where  i  =  i,2,...,N,  andN|=p, /dp.  At  such  conditions  full  burning  of  prizms  will 
be  determined  by  following  constraints:  for  prizms  with  radius  p.  :  e  >  p2 ;  for 
prizms  with  radius  P3 :  e  >  pj ;  for  prizms  with  radius  p,  :  e  >  pi . 

The  calculation  of  current  prizm's  surfaces  is  held  as  to  following  gene¬ 
ralized  formulas  ,  describing  coordinates  of  point  C  (y^,  =2RoCos30°;x,  :=2R„sin30°) 
and  radius  (x-x,)= +{y-Y,)- =  rU  and  areas  of  prizm's  faces  also 

S  =  Jy,dx-  Jy2dx ,  the  lengths  of  arcs  prizm's  surface  £,  =  }Vl  +  (y  )'dx,  lengths 

0  0  "‘I 

of  stright-lihe  segments  =  Y-yi  angle  of  crossing  between  them 

X, -X,  -y, 

tga  =  ,  where  a,  and  a,-  angular  coefficients  of  every  straight.  The 

l  +  aja, 

calculation  of  point  G  coordinates  (Fig. 3)  shows  that  circles  with  centers  in 
points  K  and  M  are  not  touch  in  this  point  as  it  takes  place  in  7-hall  grain. 


Subcontractory  points  have  coordinates  G(1,0782R,;3,85196R„),  G'(0,99216R,;3,875R„), 
that  speaks  about  presence  of  narrow  bridge  between  outside  and  middle  prizms. 
Difference  in  coordinates  is 


AX,,  =  0,08604R, ;  Ay  c  -  0,02304R,  .On  the  base 
of  expounded  above  it  may  be  drawen 
the  conclusion:  19-halls  powder  grain 
at  the  collapse  moment  of  inner  prizms 

does  not  provide  simultaneously  collapse  ^  ^ 

of  12  outside  prizms.  The  presence  of 

web  with  thickness  A=./fflJ+^=0,0891Ro  »  element  of  19-halls  powder 

that  constitute  about  9%  from  initial  die  moment  of  premier  collapse 

:  •  I  •  *  ■  r  -7  1  11  prizms  (a)  and  at  the  moment  of  secondary 
radius  up  to  join's  point  of  7-halI  collapse  prizms  (b) 
powder  gram  leads  to  form  of  prizms 

19-halIs  grain  shown  on  Fig.4-a.  The  view  of  the  same  prizm  at  the  moment  of 
it  secondary  collapse  (after  that  web  was  burned)  is  shown  on  Fig.4-b.  Full 
collapse  of  complex  prizm  occurs  at  the  moment  when  thickness  equal  to  A/2, 
i.e.  0,04555R,,  will  bum  over.  At  R,>1,04455R„  the  powder  grain  breaks  up  and 
points  G  and  G’  merge  into  one  overall  point  G  with  coordinates  =  l,0356Ro , 
yc  =  3,8636R„ .  Mentioned  peculiarity  of  collapse  were  considerated  at  the  deve¬ 
lopment  of  the  calculation  program  for  ECM.  Below  are  presented  the  results 
of  computer  calculations  of  current  burning  surface  and  relative  mass  of  gases 
for  7  and  19-hall  grains  versus  parameter  z. 

The  7-hall  grain.  The  calculations  of  perforated  tube  carried  out  as  M.E. 
Serebriakov  and  suggested  methods  for  powder  grain  with  burning  web 
thickness  2ej  =1,2  mm  and  type  sizes  of  hall  diameters  d=e,,  outside  diameter 
of  rosette  D=4ei-^3d=^6  mm  and  length  of  tube  2c  =  2,5D  =  16,5  mm.  The 
calculation  results  of  grain  geometric  form  function  coefficients  before  and  after 
the  halls  will  meet  each  other  and  parameters  of  progressivity  and  end  of 
burning,  computated  as  to  method  prof.  M.E.  Serebriakov  are  given  in  tabl.l. 
The  character  y  and  a  changing  versus  z  during  the  post-collapse  period 
collapse,  i.e  from  value  z=l  to  ,  is  shown  on  Fig.5. 


.9  1.1  1.2  1.4 


Fig,5  Dependence  of  a  and  v|;  for  7-halls  Fig.6  Dependence  of  o  and  v|/  of  19-halls 
powder  grain  versus  thickness  z  powder  grain  versus  thickness  z 

1-a  small  prizm;  2-abig  prizm;  3-sumraary  change  of  burning  surface;  4-cron 
M.E.Serebriakov's  method;  5-  vj/  on  offered  method;  6-  i|/  on  M.E.Serebriakov's  method;  B- 
the  burning  end  of  small  prizm;  A-  point  of  secondary  grain  collapse 
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19-halI  grain.  The  calculations  were  held  as  to  method  M.E.  Serebriakov  and 
the  supposed  method  for  powder  grain  with  thickness  of  burning  web  1 ,2  mm 
too.The  calculation  results  of  geometrical  form  function  coefficients  before  and 
after  collapse  and  also  parameters  of  progressivity  of  burning  computated  as  to 
M.E.Serebriakov  method  are  given  in  Table  1  and  on  Fig. 6. 

From  pictures  anyone  can  see  that  19-hall  grain  has  more  sugnificant 
burning  progressivity  in  comparison  with  7-hail  grain  and  hence  the  less  value  of 
maximal  pressure  in  comparison  with  the  last.  For  comparison  of  powders  of 
progressive  and  degressive  burning  ballistic  calculations  of  152  mm  caliber 
cannon  on  both  methods  were  carried  out  by  means  of  KOBRA  special 
program  package.  The  qualitative  results  of  influence  of  1,  7  and  19-halls 
powder  grain  on  pressure  and  projectile  velocity  are  shown  on  Fig.  7  and  8. 


A'  1 

/  Y 

/TV 

_ - 

>/ 

A/ 

t 

Fig.7  Influence  of  progressivity  of  burning  on  Fig.8  Influence  of  progressivity  of  burning  the 
pressure  projectile  velocity 

1- degressive  burning;  2- progressive  burning  of  7-hall  grain;  3- progressive  burning  of  19- 

halls  grain 

The  analysis  of  calculation  results  shows  that  more  degressive  powder 
(one-hall)  gives  more  gases  and  the  higher  level  of  maximal  pressure.The 
application  of  progressive  gun  powder  decreases  maximal  pressure  and  increase 
the  muzzle  pressure.  This  is  accordingly  leads  to  the  decreasing  of  muzzle 
velocity  (Fig.8).  Such  character  of  alteration  of  ballistic  parameters  permits  to 
provide  the  desired  muzzle  velocity  at  given  level  of  maximal  pressure  on 
account  of  selection  of  combined  (composed)  charge  characteristics.  For 
example,  the  division  of  charge  mass  between  degressive  and  progressive 
burning  components  in  ratio  0,4:0, 6  allows  to  decrease  maximal  pressure  on 
14%  while  the  muzzle  velocity  decrease  on  2,3%  only. 
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THE  TRANSVERSAL  FLUCTUATION  OF  PRESSURE  IN 
TRANS-PISTON  VOLUME  OF  BALLISTIC  LAUNCHER 


V.  A.  GOLUBEV,  V.  F.  ZAKHARENKOV,  Baltic  State  Technical  University, 
198005,  Russia,  Saint-Petersburg 


Knowledge  of  pressure  fluctuations  in  trans-piston  volume  and 
parameters,  defining  their  appearance  in  the  series  of  applied  tasks  is  more 
essential,  as  the  lasts  can  accord  noticeable  influence  on  the  character  of 
piston  motion  in  the  bore. 

In  check  cross-section  experimentally  recorded  (Fig.l)  the  presence  of 
sudden  change  local  transversal  pressure  gradients.  The  last  ones  as  to 
absolute  value  account  for  10-15%  from  maximal  value  of  average  ballistic 
pressure.  The  explanation  of  appearance  this  phenomena  may  be  given  on  the 
basis  of  fluctuations  modeling. 


0  2  4 


S  =  a  b  -  channel  cross-section 
EStEfflS  -  charge 
Mon  -  piston 


t ,  ms 


Fig.  1  Pressure  drop 


Fig.2  Structural  model 


Physical-mathematical  model  is  built  on  the  representation  of  trans¬ 
piston  volume  of  the  ballistic  set-up  as  a  set  of  longitudinal  layers  with  sub¬ 
volumes  Wi  (Fig.2).  Correspondingly  the  piston  is  represented  as  a  set  of  local 
sub-pistons  too.  Mass  of  each  sub-piston  is  qj  and  there  number  is  equal  to 
the  number  of  selected  layers.  Mass  and  energy  transfer  through  this  contact 
planes  takes  place  due  to  pressure,  temperature  or  density  gradients.  Gas 
parameters  inside  each  layer  is  taken  as  average  on  sub-volume.  Division  of 
powder  volume  on  a  number  of  rectangular  sub-volumes  demands  to  take 
account  of  disturbance  wave  propagation  rates  from  one  layer  to  another  [  1  ]. 
This  circumstance  leads  to  the  necessity  of  time  step  choice  At ,  for  integration 


of  ballistic  equations,  satisfying  to  condition  At  =!;.( 


Pi  NCi 


Velocities 


and  travels  of  sub-pistons  on  each  integration  time  step  is  calculated  from 
ballistic  parameters  in  every  i-layer.  However  stiff  connection  between  local 
pistons  make  to  tie  on  temporary  time-layer  values  of  velocity  and  paths  of 
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local  pistons  on  considered  transversal  distance  to  the  average  velocity  .and 
path  of  pistons  on  the  whole  by  means  of  equations; 


dVj  _  Sjpj 

dt  cpjq,. 


■  1  N  1  N 

^  =  ;  v:=-yv,  ;  ;  Vr‘=V; 

dt  ‘  N '  N ‘  ■ 


=  £■, 


P;=(f;O,V,.-0,5e.(p,q,V.^)/Wi;Ap.=P,-P^.,  ;  Ap,  =  -  p._,  ;  A®,  = 

Ci  p 

In  given  equations  are  admitted  following  signs:  -  the  relative  weight 
of  burned  part  of  powder  charge;  z  -  the  relative  thickness  of  burned  part  of 
powder  charge;  Jj^  -  the  burning  impulse  of  gunpowder  ;  p-  pressure  ;  A®  - 
mass  of  gas,  removing  from  layer  to  layer;  s  -  transversal  cross-section  area  of 
layer;  £  -  transient  way  of  piston(s);  V-  piston(s)  velocity;  9  -  coefficient  of  the 
accounting  of  secondary  works;  t  -  time;  At  -integration  time  step;  f-  powder 
energy  ("force"  of  gunpowder);  e-k-1-  adiabatic  coefficient  ;  q-  mass  of 
piston(s);  -  gas  density;  c  -  speed  of  sound  in  the  layer(s);  F-  contact  surface  of 
layers;  N-  total  number  of  layers;  i-number  of  layer. 

Worked  out  mathematical  model  was  realized  as  computer  program  for 
ECM  IBM  PC.  Below  are  presented  the  results  of  numerical  examinations  for 
the  case  of  two-layer  transversal  division  of  trans-piston  volume  for  set-up 
with  rectangle  channel  with  the  following  data: 

s,  =  s,  =61  sm- ,  ffl,  =  ®2  =  5  kg,  q,  =  Qj  =3  kg  and  time  step  less  then  lO"’  s. 

The  initial  stage  of  examinations  was  devoted  to  determination  of 
ballistic  parameters  in  trans-piston  volume  under  the  same  constructive  and 
technological  features  and  loading  conditions  in  sub-volumes.  As  a  result  the 
reference  variant  of  solution  without  the  accounting  of  transversal  pressure 
fluctuations  and  identical  burning  of  powder  charge  in  layers  was  found. 
Later  on  it  was  used  for  the  examinations  of  one-percent  deviations  mass  of 
charges,  burning  impulses,  force  and  density  of  powder  influence  on  values  of 
maximal  and  finite  pressure  drops. 

Fig.  3  shows  the  curves  p=f(t)  and  Ap=f(t)  for  the  case  of  one-percent 
deviation  of  "force"  of  powder  only.  Anyone  can  see  that  the  duration  of 
pressure  fluctuations  is  comparable  with  the  time  of  piston  moving  along  the 
bore. 


Fig.3  Gas  pressure  change  inside 
sub-volumes 


Fig. 4  Comparison  of  pressure  calculation 
with  experimental  data 
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The  particular  interest  presents  the  examination  of  transversal  pressure 
fluctuations  at  the  same  parameters  of  loading  in  sub-volumes  for  modeling  of 
charge's  destruction  in  layers.  Such  investigation  allows  to  wide  representation 
about  real  processes ,  taking  place  in  trans-piston  volume  of  ballistic  sets.  So, 
for  instance,  at  the  destruction  of  charge  in  layer  1  at  y ,  =0,45  and  in  layer  2  - 
at  \|/2=0,55  with  coefficients  of  crushing  charges  as  2  and  3  accordingly, 
calculation  curv^e  AP=f(t)  satisfactorily  consists  with  experimental  curve  on 
Fig.l  not  only  in  qualitative  but  and  quantitative  relation  (Fig.4  ).  At  t=l,8 
ms,  as  like  how  calculations,  experimentally  detected  the  beginning  of 
pressure  fluctuations. 

Thus  worked  out  model  allowed: 

•  to  bring  out  the  conditions  of  transversal  pressure  fluctuations  appearance 
within  the  trans-piston  volume  of  ballistic  launchers; 

•  to  give  quantitative  estimation  of  various  factors  influence  on  the  value  of 
transversal  pressure  gradients  in  the  trans-piston  volume,  which,  depending 
on  constructive-technological  parameters  and  of  charge  loading 
conditions,  can  achieve  lO-s-15  %  from  Pm  value  and  ought  to  be  taken  into 
account  during  the  designing  and  experimental  testing  of  ballistic  sets-up; 

•  to  bring  out  the  sensitivity  of  transversal  pressure  fluctuations  to  deviations 
of  charge  parameters  and  to  conditions  of  charge  loading,  nomenclature 
and  value  of  which  can  be  specified  during  further  researches  of  specific 
ballistic  sets. 
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BRAKE  DEVICES  OF  BALLISTIC  SET-UP  WITH 
A  RING  GAS  WITHDRAW 


V.A.  GOLUBEV,  V.M.  ALEKSEEV,  VA,  VASIN 
Baltic  State  Technical  University,  198005,  Russia,  Saint-Petersburg 

There  are  many  types  of  brake  devices  of  ballistic  sets-up  utilizing  the  residual 
energy  of  gases  with  the  purpose  of  the  recoil  pulse  reduction.  Such  brake  devices  begin 
to  work  after  the  piston  pass  through  the  muzzle.  The  variety  of  such  brake  devices  can 
be  explained  by  aspiration  to  find  the  optimal  design,  which  would  provide  the  least 
load  on  the  ballistic  set-up,  the  least  perturbation  on  the  piston  inside  the  brake  device 
chamber  and  the  least  mass  of  device. 

The  questions  connected  with  asymmetrical  action  of  gases  on  the  piston  inside 
the  device  chamber  become  more  urgent  as  power  of  both  set-up  and  brake  devices 
increasing.  Widely  used  in  practice  earlier  created  perfect  brake  devices,  have  the 
constructive  characteristic  a  >  -  0.26. 

In  the  present  paper  some  results  of  research  of  brake  devices  with  the  ring  gas 
supply  from  the  trans-piston  volume  are  submitted.  Using  such  devices  one  can  gain 
effective  reduction  of  the  recoil  pulse  and  minimize  the  piston  perturbation.  The 
fragment  of  set-up  with  the  brake  device  of  such  type  is  shown  at  the  below  drawing. 

Contouring  of  the  muzzle  brake  is  aimed  at 
obtaining  of  the  directed  change  of  the  gas  flow. 
In  the  brake  device  the  gas  flow  is  divided  on 
eentral  and  peripheral  flows  proportionally  to  the 
areas  of  central  window  and  ring  entrance  in  the 
peripheral  channel.  Thus,  up  to  the  entrance  in  the 
peripheral  channel,  the  gas  flow  is  not  exposed  to 
asymmetrical  transformations  and,  accordingly,  is  not  a  source  of  asymmetrical 
perturbation  on  the  piston  inside  such  chamber. 

The  distance  Lmin  (from  the  beginning  of  peripheral  expansion  of  gases  up  to 
forward  edge  of  central  hole)  is  rather  important  for  determination  of  the  rational  size 
and  mass  characteristics  of  the  brake  device. 

For  the  two-chamber  brake  device,  with  a  ratio  of  the  areas  of  central  and 
peripheral  channels  in  cross-sections  1-land  2-2  equal  to  1:4,  size  Lmin  should  be 
not  less  then  the  diameter  of  the  central  channel.  The  efficiency  of  such  scheme  of  brake 
devices  is  checked  up  experimentally  for  Ct  =  -  0.7  as  applied  to  two-chamber  design  of 
the  brake  device.  The  calculated  pressure  distribution  on  chambers  surface,  in 
accordance  with  channels  geometry,  is  proved  to  be  true  by  measurements  of  gas 
pressure  in  points  A  and  B,  which  make  10,5  MPa  and  2,5  MPa  accordingly,  i.e.  are  in 
the  ratio  1:4,2. 
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MATHEMATICAL  MODELING  OF  A  SHOT  WITH 
DISPERSIVE  MONOBLOCK  CHARGES 

A.N.  ISCHENKO,  J.P.  KHOMENKO 

Institute  of  Applied  Mathematics  and  Mechanics  at  Tomsk  State  University 

One  of  the  ways  to  increase  the  efficiency  of  an  artillery  shot  consists  in 
use  monoblock  powder  charges,  which  in  comparison  with  usual  granulated  and 
tubular  have  some  doubtless  advantages:  more  higher  density  and  additional 
concentration  of  energy,  opportunity  of  gradually  ignition  and  controlled 
progressive  combustion  during  a  shot,  redistribution  of  energy  in  burning  space 
and  other.  It  is  probably  to  use  monoblocks  in  various  circuits  of  a  shot.  The 
most  known  and  distributed  is  the  circuit  with  the  monoblock  in  the  chamber  of  a 
trunk  (MK).  Alongside  with  this  monoblock  can  be  used  as  travelling  charge  (PZ), 
and  also  in  the  various  combined  circuits  of  a  shot  together  with  granulated  or 
tubular  charges. 

It  is  necessary  to  mean,  that  if  for  monoblocks  manufacturing  to  use  usual 
ballistit  gunpowder,  with  burning  speed  about  centimeters  in  a  second,  than 
necessary  gasflow  is  possible  to  create  by  an  only  advanced  surface  of  burning.  It 
means,  that  the  charges  will  be  porous  and  the  porosity  can  be  as  gastransfer  and 
to  be  present  at  the  monoblock  in  an  initial  condition,  and  opened,  i.e.  appearing 
during  burning  fuel.  In  any  case  the  burning  of  such  charge  will  be  accompanied 
by  dispersion  of  a  porous  matrix  on  separate  particles.  In  a  common  case  the 
received  mix  will  include  a  few  factions,  each  of  which  is  moving  and  burning 
on  its  own  laws. 

The  similar  picture  is  observed  and  in  the  usual  classical  circuit  of  a  shot.  By 
virtue  of  the  technological  and  statistical  reasons  powder  elements  always  differ 
from  each  other  and  will  form  a  poly  dispersive  mix.  Moreover  in  regular  shots  are 
as  often  used  the  charges,  consisting  from  different  grades  of  gunpowder,  and 
also  there  are  other  charges  elements:  ignitor,  burning  down  cartridge  and  etc.. 
Thus,  practically  always  at  account  of  a  artillery  shot  it  is  necessary  to  collide 
with  a  problem  of  polydispersive  gas-powder  mix.  In  known  models  of  internal 
ballistics  it  while  is  not  completely  solved.  Is  as  a  rule  used  biphase,  two-high¬ 
speed  (gas  -  gunpowder)  approximation.  The  line  of  the  authors  allocates  as  the 
third  phase  the  ignitor,  but  it  is  considered  to  move  with  the  gas  speed.  In  the 
given  work  new  mathematical  model  of  a  shot,  taking  into  account  the 
polydispersive  mix  and  opportunity  to  disperse  of  the  monoblocks  is  offered. 
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The  model  is  based  on  the  following  assumptions.  The  powder  mix 
includes  N  factions,  each  of  which  moves  with  the  own  speed  and  burns  under 
the  own  laws.  The  products  of  combustion  of  all  factions  are  identical  on 
thermodynamic  properties  and  will  form  uniform  carrying  environment.  The  gas  is 
considered  unviscous  and  unheat-transfered,  and  the  influence  of  viscosity  is  taken 
into  account  only  in  interaction  of  phases.  Particles  are  not  compressable,  and  are 
burning  under  the  georhetrical  law  /!/  with  stationary  speed.  The  transition  of 
particles  from  one  faction  in  other  is  away.  The  initial  period  of  a  shot  is  not 
taken  into  account,  heat-transfer  to  the  walls  of  the  channel  we  neglect. 

The  basic  system  of  the  equations  includes  the  laws  of  preservation  of 
weight,  pulse  and  energy  for  a  gas  phase,  laws  of  preservation  of  weight  and  pulse 
for  each  faction  of  particles,  equation  of  burning  and  change  of  number  of 
particles.  At  the  description  of  transformation  of  monoblocks  it  was  supposed,  that 
the  porous  fuel  (as  MK,  and  PZ)  partially  bums  down  in  a  plane  mode,  and 
partially  disperses  on  particles,  burning  down  in  depth  of  the  volume.  Dispersive 
products  of  PZ  will  form  only  one,  for  example,  m-th_  a  faction,  and  at  MK  — 
1-th  faction  of  a  general  mix  of  n-factions.  The  front  of  burning  is  replaced  by  a 
surface  of  gasdynamics  break,  through  which  with  identical  speed  flow  gas 
products  and  dispersive  particles. 

At  realization  of  accounts  a  grid  was  used  mobile  similar  differencive  net, 
one  end  of  which  fastened  at  the  front  of  PZ  burning,  and  second  —  at  the  front 
of  MK  burning.  At  absence  of  monoblocks  borders  of  settlement  area  served 
bottom  of  the  chamber  of  a  trunk  and  mobile  shell.  The  two-step-by-step  obvious 
circuit  was  applied.  At  a  first  stage  the  Lax-circuit  with  a  step  on  the  top 
temporary  layer,  and  at  second  stage  —  the  integrated  laws  of  preservation  was 
used.  Peculiarity  of  numerical  algorithm  was  occurrence  of  an  additional  condition 
of  stability  by  the  right  part,  connected  with  presence  of  interphases  force 
resistance. 

For  testing  offered  model  of  a  shot  and  developed  numerical  technique  two 
series  of  experiments  on  a  model  ballistic  installation  were  spent.  In  first  was  used 
combined  granulated  charge,  consisting  from  fine  single-channel  and  large  seven- 
channel  gunpowder,  and  in  second  —  plane  and  single-channel  gunpowder. 
Comparison  of  settlement  and  experimental  results  has  shown,  that  a  divergence 
on  the  maximum  pressure  of  a  shot  does  not  surpass  2  %,  and  on  initial  speed  of 
a  shell  -  1.2  %.  The  gasdynamics  analysis  of  a  picture  flow  in  the  burning  space 
has  shown,  that  the  larger  faction  of  seven-channel  gunpowder  is  burning  longer 
single-channel  and  takes  a  great  interest  less  for  a  shell. 

On  the  basis  of  the  advanced  approach  are  investigated  gasdynamics 
peculiarity  of  the  combined  shot  circuit  with  use  of  dispersive  PZ.  With  this 
purpose  a  series  of  experiments  is  carried  out  and  comparison  settlement  and 
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experimental  dates  is  spent.  The  parameters  of  the  coordination  were  served  by 
dispersive  speed  and  characteristic  of  dispersive  phase.  Dynamics  of  dispersive 
products  area  is  investigated  and  is  shown,  that  the  flow  of  powder  gases  from 
PZ  carries  dispersive  faction  to  the  middle  of  volume,  where  it  mixs  up  with 
particles  of  granulated  gunpowder.  The  zone  of  mixture  does  not  reach  bottom  of 
the  chamber  as  fast  dispersal  PZ  and  shell,  connected  with  increase  of  pressure 
and  reduction  of  PZ  weight  begins.  The  back  front  of  particles  takes  a  great 
interest  for  the  PZ  and  up  to  a  moment  of  complete  dispersion  PZ  "is 
accompanied"  by  a  dispersive  zone  in  width  0.2-0.25m.  After  PZ  complete 
dispersion  between  bottom  of  a  shell  and  the  dispersive  zone  will  form  area  of 
pure  gas,  which  in  due  course  extends.  Thus  forward  front  of  m-th  faction 
particles  is  designated.  To  a  moment  of  an  exit  of  a  shell  from  a  trunk  it 
practically  completely  burns  down. 

Parametrical  researches  on  the  basis  of  offered  model  have  allowed  to 
establish,  that  the  circuit  of  a  shot  with  MK  use  can  be  applied  to  modernization 
hunting  and  sporting  shot  flL  Thus  is  shown,  that  at  presence  at  MK  fuel  of 
necessary  delay  of  ignition,  it  is  possible  not  only  in  common  to  use  with 
granulated  gunpowder,  but  also  to  keep  thus  weight  regular  charge.  The  given 
effect  has  completely  proved  to  be  true  at  realization  of  experimental  tests,  during 
which  the  speed  of  fraction  is  increased  up  to  20  %  without  increase  of  the 
maximum  pressure  of  a  shot. 

1,  CepeOpaKOB  M.E.  BHyipenHaa  riaiumcTHKa  CTBOJibHwx  CHCieM  h 
nopoxoBtix  paxeT.  M.:  06opoHrn3,  1962,  703c. 

2.  Ischenko  A.N.,  Starchenko  A.V.,  Khomenko  J.P.  The  peculiarities  of 
combustion  of  granulated  and  combined  powder  charges  in  a  sporting  gun.  Abst. 
Int.  Worksh.  On  Chem.  Gasdyn.  And  Comb.  Energ.  Mater.  Tomsk,  1995,  p.85-86. 


22 


METHOD  and  RESULTS  of  the  DECISION  of  a  RETURN 
PROBLEM  of  INTERNAL  BALLISTICS 


EROFEEV  VK.,  GRIGORIEV  V.  V. 

BSTU,  St. -Petersburg 

Non-stationary  problem  on  definition  of  the  necessary  laws  of  change  gas-and 
thermodynamic  parameters  in  combustion  chamber  of  in  volume  of  any  structures. 
The  problem  discusses  with  the  account  of  physics  and  chemical  conversions, 
burning  surface  of  a  solid  fuel,  assuring  given  laws  of  changes  of  parameters  in  vol¬ 
ume,  kinematics  and  dynamic  laws  of  movement  of  the  object. 

There  are  known  the  decisions  of  a  number  of  return  problems  of  internal  balli- 
tics.  As  a  rule,  there  are  the  problems  of  PAD  designing,  which  have  to  guarantee 
the  necessary  ejection.  As  the  other  turn  there  are  the  problems  of  volumes  desiging 
in  structures  for  creation  of  tlie  required  laws  of  parameters  changing.  The  works 
devoted  of  deciding  the  complex  return  problem  and  for  the  volume,  where  com¬ 
bustion  products  from  the  PAD  are  flow  in  are  not  known  by  the  authors.  The  same 
words  one  can  say  about  the  PAD’S  combustion  chamber  in  view  of  the  possibihty 
to  account  the  physics  and  chemical  reactions  of  gas  mix  for  the  different  object’s 
assignment  laws. 

The  mathematical  model  of  such  complex  return  problem  consists  of  two  equa¬ 
tion  systems,  describing  processes,  occurring  in  given  volume  (1-9)  and  in  a  com¬ 
bustion  chamber  (10-14). 
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At  development  of  mathematical  models  the  following  assumptions  were  accepted: 

-  the  products  of  combustion  represent  an  ideal  gas; 

-  the  concentration  changing  of  gas’s  mix  substances,  occurs  at  the  expense  of  a 
solid  fuel  burning  and  at  the  physics  and  chemical  reactions; 

-  gas-  and  thermodynamic  parameters  in  a  PAD’s  combustion  chamber  are 
rounded. 

The  equation  system  is  realized  as  a  calculation  program  for  IBM  computers. 
The  (1-9)  equation  system  is  decided  by  Runge-Kutt  method.  The  (10-14)  system  is 
rigid  and  for  its  decision  the  Girs  method,  enabling  for  easy  integrating  such  systems 
is  used. 


The  (10-14)  system  for  approximate  valuation  of  burning  surface  can  be  lead  to 
the  next  algebraic  equation: 


wh«e 


_Gk  r  _  _  ^*^10 


A 


The  numerical  researches  results  permit  to  generalize  about  serviceability  of  an 
offered  technique  and  opportunity  of  using  with  considerably  smaller  costs,  than  at 
the  decision  of  a  direct  problem,  in  design  work. 
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THE  ACCOUNT  OF  PERCUSSION  CAP  PRESSURE  WITH 
BALLISTIC  CALCULATIONS  OF  THE  SMALL  ARMS 
AMMUNITIONS  WITH  LOW  BALLISTICS 

MAKAVEYEV  A.  T„  BALTIC  STATE  TECHNICAL  UNIVERSITY, 

198005,  S.-  PETERSBURG,  RUSSIA 

The  initial  pressure  in  chamber  renders  big  influence  on  a  course  of 
ballistic  process  in  the  channel  of  a  barrel.  Particulary,  the  quantity  of  the  powder 
which  has  burned  during  a  shot  depends  from  one  for  the  weapon  with  low 
density  of  loading.  With  realization  of  ballistic  calculations  of  artillery  systems  the 
prime  device  pv  is  accepted  as  initial  pressure  which  will  be  formed  as  a  result  of 
smoky  gunpowder  combustion  of  the  prime  cartridge  but  the  pressure  of  shock 
structure  of  percussion  cap  is  not  taken  into  calculations  as  it  is  a  very  small  in 
comparison  with  pressure  from  smoky  gunpowder  of  prime  device  and  main 
charge  of  smokeless  powder.  However  there  is  a  necessity  of  the  account  of 
percussion  cap  pressure  with  ballistic's  calculations  of  stopping  action 
ammunitions  and  hunting  ammunitions.  It  is  caused  by  the  several  reasons: 
throwing  bodies  of  this  ammunitions  have  the  small  meaning  of  parameter  Cq,  but 
starting  pressure  po  and  shock  structure  of  percussion  cap  pressure  pvk  are 
comparable  to  pressure  from  burning  of  powder  p.  Today  two  kinds  percussion 
caps  applies  for  ignition  of  powder  charge  in  hunting  ammunition-  opened  and 
closed.  Cap  of  "  Tsentroboy"  type  concerns  to  first,  and  second  type  is 
"Zhevelo".  Cap  "Tsentroboy"  is  less  powerful  on  own  action  and  is  usually  used 
with  equipment  of  smoky  gunpowder  ammunition  as  this  powder  is  easily  ignited. 
It  is  more  preferable  to  use  cap  "Zhevelo"  with  equipment  ammunition  by 
smokeless  pbwder.  Cap  "Zhevelo"  is  applied  also  in  ammunition  of  stopping 
action  of  23-mm  rifle  shot  gun  KS-23,  which  was  accepted  on  arms  of  interior 
security  troops  in  1991. 

Percussion  cap  "Tsentroboy"  (State  Standard  N  7574  -  71)  exists  two 
versions:  TsBO  -  with  shock  structure  on  a  basis  of  Hg(ONC)2  and  TsBO-N  - 
with  anticorrosive  structure  on  a  basis  of  C6H(N02)30Pb‘H20.  Prescriptions  of 
shock  structures  of  these  caps  are  following  (in  %  on  weight): 


TsBO 

% 

TsBO-N 

% 

Hg(ONC)2 

16,5  ±1,5 

CfiHCNO^l^OPb-H^O 

35  ±  3 

KCIO, 

55,0  r  3,0 

CfiHsON.o 

3  -t  1 

Sb,0, 

28,5  +  1,5 

Ba(N03)2 

42  ±  3 

PbO, 

5  ±  2 

Granule  mix  of  Sb^O^  and  A1 

15+2 
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Weight  of  shock  structure  in  "Tsentroboy"  type  percussion  caps  is  0,03  ± 
0,002  g.  Percussion  cap  "Zhevelo"  (State  Standard  N  24579  -  80)  also  exists  two 
versions:  "Zhevelo-N"  (anticorrosive)  -  with  structure  on  a  basis 
C6H8(N02)30Pb'H20  and  "Zhevelo-M"  (powerful)  -  with  structure  on  a  basis 
Hg(ONC)2.  Prescription  of  shock  structure  of  cap  "Zhevelo-N"  is  the  same  as  at 
cap  TsBO-N,  and  prescription  of  structure  of  "Zhevelo-M"  is  following: 

Hg(ONC)2  -  50  %; 

KCIO3  -  33  %; 

SbiOs  -  17  %. 

Weight  of  shock  structure  in  percussion  caps  of  this  type  is  0,05  ±  0,003  g. 
Percussion  caps  with  shock  structure  on  a  basis  Hg(ONC)2  have  got  the  greater 
practise  than  anticorrosive  structures  because  of  their  high  energy,  in  particular 
for  ammunition  of  pomp  gun  KS-23  where  "  Zhevelo-M"  cap  is  only  applied. 
Prescription  of  shock  structure  renders  big  influence  on  the  physical  characteristics 
of  products  of  combustion.  The  physical  data  for  two  shock  structures  on  a 
Hg(ONC)2  are  resulted  below. 


_ ^ _ Ca2 _ 

TsBO 

Zhevelo  -  M 

Specific  volume  of  gases  with  normal  conditions,  10'^  /  kg 

107,2 

210 

Flame  temperature  of  gases,  °  C  ' 

4600 

4500 

Adiabatic  coefficient 

1,11 

1,27 

Quantity  of  gases  in  products  of  combustion,  % 

34.3 

68,62 

Quantity  of  the  solid  rest,  % 

65,7 

32,28 

Calculation  of  pressure,  created  percussion  cap,  is  possible  to  make  under 
the  formulas  of  classical  internal  ballistics.  The  pressure,  which  is  created  by 
products  of  combustion  of  shock  structure  in  volume  of  cap,  is  determined  on  the 
formula; 


Pk  =  fk  Wk yg/  (Wk -  Wk yt/  dt) ,  (1) 

where  w^  -  weight  of  shock  structure;  f^  -  force  of  shock  structure;  Wk  -  internal 
volume  of  cap;  dt  -  density  of  the  solid  rests  of  products  of  combustion  shock 
structure;  yg  -  relative  quantity  gas  products  of  shock  structure  combustion; 
yt  -  relative  quantity  of  the  solid  rests  of  shock  structure  combustion. 

The  force  of  shock  structure  of  percussion  cap  can  be  calculated  as  ; 


fk  =  Rk  Tik  , 


(2) 
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where  R^-  specific  gas  constant  for  a  mix  of  gas  products;  Tik  -  explosive 
temperature  of  shock  structure. 

Gas  constant  for  a  mix  of  combustion  gas  products  is  defined  from  the 
energy  equation  of  Clapeyron-Mendeleev: 

Rk  =  Pa  Wik  /  To ,  (3) 

where  pa  =0,1013  MPa  -  normal  atmospheric  pressure;  To  =  273  K;  Wik  - 
specific  volume  of  gases  with  normal  pressure  and  temperature  To  =  273  K. 

Pressure,  which  is  produced  by  products  of  combustion  of  shock  structure 
in  free  volume  of  a  chamber,  is  possible  to  define  from  the  equation  of  adiabatic 
process: 

Pvk  =  Pk  ((  Wk  -  Wkyt/dt)/(Wo  -  w/d  +  Wk  ■  Wkyt/dt) ) ‘‘j  (4) 

where  Wo  -  initial  volume  of  chamber;  w  -  weight  of  powder;  d  -  density  of  a 
grain  powder  ,  k  -  adiabatic  coefficient. 

The  calculations  under  the  formulas  (1),  (2),  (3)  and  (4)  for  two  caps  have 
given  the  following  results: 


Cap 

TsBO 

Zhevelo-M 

Gas  constant,  Rj;,  J/(kg  K) 

39,8 

77,9 

Force  of  shock  stmcture,  fk,  MJ/kg 

0,194 

0,372 

Pressure  of  explosion  products  in  cap  volume,  Pk,  MPa 

30,7 

170,2 

Pressure  of  explosion  products  in  free  volume  of  a  chamber, 
MPa _ 

1,5 

6,3 

The  ballistic  calculations  of  23-mm  standard  ammunition  of  shot  gun 
KS-23  have  shown  that  if  not  to  take  into  calculation  pressure  from  gases  of  shock 
structure  of  percussion  cap  then  the  mistake  with  calculation  of  muzzle  velosity 
makes  from  15  up  to  30  percents  depending  on  density  of  loading. 


REFERENSES 

1.  nX^T  24579-81.  KancrojiH-BociuiaMertHrejiH  "^eBeJio".  -  M.,1981. 

2.  rOCT  7574-71.  KanciojiH-BocnjiaMeHHTejiH  itempajibHoro  6oh  juisi  naipoHOB 
oxoTHHHbHx  pyxeH.  -  M.,  1993. 

3.  KapnoB  n.  fl,  KanciojTH-BocnJTaMeHHTejiH  h  xanciojiH-^^eTOHaTopbi.  -  Jl.,  1935. 


Section  2 

INTRA-CHAMBER  PROCESSES 


30 

MATHEMATICAL  MODELING  OF  THE  IGNITION 
ABNORMAL  MODE  DEVELOPMENT 
IN  THE  SOLID  PROPELLANT  ROCKET  MOTOR, 
CAUSED  BY  THE  NONESTIMATED  ACTIONS 

Alexey  M.  LIPANOV,  Alexander N.  LUKIN  and  All  V.  ALIYEV 

Institute  of Applied  Mechanics,  Ural  Branch  of  the  Russian  Academy  of  Sciences 
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The  tendency  of  development  of  the  modem  large-sized  solid  propellant  rocket 
motors  (SPRM)  is  characterized  by  increase  of  the  combustion  chamber  filling 
degree  by  propellant  and  operation  loads  level.  Improvement  of  the  modem  large- 
sized  SPRM  characteristics,  reduction  of  terms  and  material  expenses  on  their 
improvement  and  manufacturing  depends  largely  on  reached  results  in  understanding 
the  regularities  of  a  course  and  possible  accuracy  of  the  anomalous  physico-chemical 
processes  prediction,  which  arise  during  the  operational  period.  There  are  a  class  of 
the  large-sized  SPRM  with  a  high  intrachamber  loading  coefficient  (  >  0,95  ).  The 
main  peculiarity  of  the  stmctural  diagrams  of  such  propulsion  systems  is  that  the 
firmly  fastened  internal  channel  charge  has  a  partially  non-fastening  and  unarmoured 
end  face  surfaces  in  the  vicinity  of  the  front  bottom.  In  these  conditions  the  end  face 
combustion  surface  can  make  up  to  50  %  and  more  from  the  whole  total  combustion 
surface.  The  charge  end  face  is  usually  fastened  with  the  engine  front  bottom 
with  using  of  the  boot  (the  sealing  ring).  On  the  fig.  1  the  typical  stmctural  diagram 
of  such  SPRM  is  shown  ( 1  -  is  the  fiber  glass  case;  2  -  is  the 


Figure  1. 


Figure  2. 
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pyrotechnic  ignition  system;  3  -  is  the  nozzle  bottom;  4  -  is  the  boot  (the  sealing 
ring) ),  ITie  other  peculiarity  of  this  class  of  SPRM  is  that,  in  normal  conditions,  at 
which  the  propulsion  systems  storage  is  provided,  the  clearances  width  between  the 
engine  case  front  bottom  internal  surface  and  the  charge  end  face  surface  A  <  0.001  m. 
In  the  process  of  the  intrachamber  volume  filling  by  the  high-temperature  combustion 
products  (CP),  coming  from  the  ignition  system  (IS)  and  from  the  propellant  already 
ignited  part,  there  are  joint  elastic  deformations  of  the  propellant  charge  and  engine 
case.  The  engine  case  has  a  “cocoon”  type  design  and  is  executed  from  the  organic 
plastic  materials.  In  these  conditions,  the  thickness  of  the  not-flowing  clearance 
between  the  SPRM  case  front  bottom  internal  surface  and  the  charge  end  face  surface 
can  increase  in  20  -  100  times  at  the  operating  pressures  in  the  combustion  chamber 
5-10  MPa.  On  fig.  2  is  submitted  the  dependence  of  the  front  face  clearance  cross 
section  area  in  a  vicinity  of  the  engine  front  bottom  from  the  longitudinal  coordinate  in 
various  time  moments;  1  -  t  =  0  sec.,  2  -  0.036  sec.,  3  -  0.048  sec.,  4  -  0.060  sec. 
Here  the  significance  of  coordinate  x  =  0  corresponds  to  the  boundary  between  the 
end  face  cavity  and  the  engine  front  volume.  In  a  number  of  cases,  listed  peculiarity 
of  the  considered  SPRM  can  cause  the  increased  sensitivity  of  the  SPRM  operation 
initial  stage  to  a  various  kinds  of  the  nonestimated  actions.  In  general  case,  under 
the  "nonestimated  actions"  are  implied  not  only  deviations  of  working  parameters  of 
the  separate  units  and  elements  of  the  SPRM,  but  also  supposed  design  errors, 
violations  of  technological  process  modes  and  aging  of  the  materials  of  the  elements, 
both  of  the  design  of  SPRM,  and  the  charge. 

In  a  considered  situation,  the  reason  of  an  anomalous  mode  of  ignition  and 
combustion  in  the  SPRM  combustion  chamber  beginning  can  serve  the 
nonestimated  action  from  the  pyrotechnic  IS.  Such  noncalculated  effect  result  in 
nonmonotonicity  of  the  working  processes  course  during  the  propellant  charge 
ignition. 

In  a  number  of  cases,  realization  of  the  situation  is  possible,  when  from  the 
moment  of  beginning  of  the  high-temperature  CP  outflow  from  the  IS  and  up  to  the 
moment  of  the  flame  occurrence  on  the  propellant  charge  surface,  mentioned 
nonflowing  end  face  clearance,  having  significant  elongation,  has  no  time  to  expand 
up  to  the  necessary  sizes.  Beginning  of  a  similar  situation,  in  the  most  responsible 
period  of  activity  of  this  class  of  the  SPRM,  can  initiate  a  number  of  anomalous 
physico-chemical  processes,  which  capable  to  affect  essentially  on  the  engine  further 
operation.  The  most  probable  sequence  of  the  ignition  anomalous  mode  development 
can  look  like  as  follows.  The  clearance,  representing  a  poorly  warmed  up  stagnant 
zone,  will  hinder  to  supply  of  enough  heat  quantity,  necessary  for  the  charge  front  end 
face  ignition.  And  in  case  of  absence  of  beforehand  organized  distribution  of  the  CP 
flows  from  the  IS  on  directions  can  result  in  much  greater  irregularity  of  the  heat 
supply  to  the  charge  various  surfaces.  It  will  lead  to  increase  of  the  front  end  face 
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ignition  time  delay  and  to  the  nonmonotonicity  of  all  process  of  the  charge  ignition. 
Besides,  in  the  beginning  of  IS  functioning,  the  CP  flowing  into  as  from  the  IS,  and 
from  charge  sections  already  joined  to  combustion,  will  affect  on  the  initial  gas, 
placed  in  the  end  face  clearance,  as  a  piston.  Thus,  into  the  depths  of  the  end  face 
clearance  will  begin  propagate  the  compression  waves,  which,  as  a  rule,  are 
transformed  into  an  intensive  shock  wave.  The  last,  when  reaches  the  area  of  interface 
of  the  engine  bottom  with  the  boot  and  the  charge  end  face,  is  reflected  at 
simultaneous  build  up  pressure  in  this  place.  Then,  propagation  of  a  reflected  shock 
wave  begins,  and  the  initial  gas,  which  contiguous  to  the  area  of  interface  of  the 
engine  bottom  with  the  boot  and  the  charge  end  face,  will  contract  up  to  the  high 
pressure.  The  considered  phenomenon  of  gradually  strengthening  compression  waves 
propagation  was  for  the  first  time  investigated  in  the  ordnance  science,  in  1876  -  1878 
by  Major-General  of  the  Russian  Army  Nickolay  V.  Kalakoutsky,  the  outstanding 
Russian  scientist  in  the  fields  of  metallurgy,  physical  metallurgy  and  ordnance.  As  a 
result  of  sequential  development  of  an  indicated  phenomenon,  in  the  area  of  interface 
of  the  end  face  clearance  cavity  with  the  boot,  which  fastening  the  charge  end  face 
with  the  engine  bottom,  can  arise  a  strong  local  dynamic  loads,  causing  either  to  the 
boot  breaking-off,  or  to  the  boot  breaking-off  from  the  engine  case  and  from  the 
charge.  The  numerical  analysis  shows  [  1  ],  that  at  the  clearance  value  equal  to  A  = 
0.001  m.  the  pressure  maximum  level  in  the  end  face  clearance  cavity  is  surpasses 
more  than  in  13  times  the  pressure  level,  corresponding  to  the  clearance  value  equal 
to  A  =  0,1m. 

The  boot  breaking-off  (or  rupture)  from  the  engine  case  is  accompanied  by 
opening  of  a  new  cavity,  in  which  the  incandescent  CP  begins  inflow.  The  last 
circumstance  change  over  the  anomalous  process  development  on  a  new  qualitative 
level,  which  characterized  already  more  large  scale  action  on  the  intrachamber 
processes  in  the  SPRM.  The  further  development  of  the  non-stationaiy  physico¬ 
chemical  processes  can  represent  hazard  for  the  normal  functioning  of  the  SPRM  in 
aU.  In  case,  if  the  surface  on  which  there  is  the  boot  breaking-off  is  great,  the 
infraballistics  processes  development  in  the  engine  chamber  can  happen  on  a 
following  way.  The  CP  pressure  in  the  vicinity  of  the  boot  breaking-off  place  on  a 
some  time  interval  can  decrease,  and  on  the  SPRM  chamber,  as  the  corollary,  will  be 
propagated  the  rarefaction  wave.  In  case  of  the  strong  rarefaction  wave  appearance, 
can  take  place  either  extinction  of  the  propellant  charge  (complete  or  temporary)  or 
transition  to  the  propellant  charge  unsteady  burning  mode  (the  SPRM  “chuffing”)  [  2 
].  Propagation  of  the  compression  wave  along  the  cavity,  formed  at  the  boot  breaking- 
off,  can  result  in  ignition  of  the  propellant  opened  surface  and  to  create  dangerously 
large  pressures  in  the  vicinity  of  the  upper  boundary  of  an  indicated  cavity.  And  it,  in 
turn,  can  result  in  the  propellant  charge  destruction.  The  investigations  of  the  flame 
propagation  processes  along  the  propellant  surface,  which  placed  in  the  narrow  deaf 
and  deformable  channel,  through  the  opened  boundary  of  which,  from  the  IS  the  CP 
flow  into  inwards,  were  conducted  in  the  papers  [  1  ],  [  3  ].  The  difference  of  the 
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present  research  from  indicated  papers  is,  that  here  is  considered  the  situation,  when 
on  the  boundary  between  the  narrow  deaf  end  face  clearance  with  boot,  during  boot 
breaking-off  a  new  cavity  is  opened.  This  cavity  also  representing  a  deaf  channel. 

The  numerical  modelling  of  anomalous  physico-chemical  processes  development  in 
the  SPRM  chamber  after  the  boot  breaking-off  is  carried  out  in  the  following 
statement.  The  CP  flow  in  the  propellant  charge  central  channel  and  in  the  end  face 
clearances  is  considered  in  one-dimentional  gasdynamics  statement.  Thus  it  is 
supposed,  that  the  CP  are  the  chemical  nonreacting  mixture  of  several  components. 
The  heat  flows  from  CP  in  the  propellant  charge  and  the  SPRM  case  are  determined 
with  use  of  criteria  relationship  [  1  ].  Warm-up  of  the  propellant  surface  and  its 
ignition  is  determined  on  the  models  of  Professor  R.  Ye.  Sorkin  of  Russia.  It  is 
supposed,  that  the  boot  breaking-off  happens  instantly,  and  the  parameters  in  the 
vicinity  of  breaking-off  are  described  in  the  one-dimentional  statement.  The 
gasdynamics  equations  are  solved  numerically  by  the  “large  particles”  method  [  4  ], 
which  developed  by  academician  Yu.M.  Davydov  of  Russia,  with  utilization  of  a 
number  of  modifications  [  1  ]. 

For  the  complex  check-up  of  the  developed  mathematical  model  of  the  physico¬ 
chemical  processes,  the  numerical  calculations  results  have  been  compared  with  the 
fire  stand  tests  results.  The  pressure  gauge  was  installed  in  the  engine  front  bottom 
area,  namely  in  the  IS  flange  plug  ( fig.  1,  (2) ). 
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Figure  3.  Figure  4. 

On  the  fig.  3  presented  the  results,  which  corresponds  to  the  processes  development 
without  boot  breaking-off:  1  -  is  the  calculated  CP  pressure  on  the  boundary  between 
the  front  end  face  clearance  and  the  boot;  2,  3  -  are  the  computational  and 
experimental  CP  pressure  in  the  SPRM  front  volume.  On  the  fig.  4  presented  the 
results,  which  corresponds  to  the  processes  development  at  boot  breaking-off:  1  -  is 
the  calculated  CP  pressure  on  the  boundary  between  the  front  end  face  clearance  and 
again  opened  cavity;  2,  3  -  are  the  computational  and  experimental  CP  pressure  in  the 
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SPRM  front  volume.  During  the  calculations,  the  volume  of  the  opened  cavitny 
relied  equal  to  Wo*^''^~0.13m^.  Thus,  the  boot  breaking-off  calculated  pressure  has 
appeared  equal  to  the  value  P*‘'®’^s3.225  MPa,  and  the  breaking-off  pressure  value, 
obtained  by  an  experimental  way  has  made  . 1 53  MPa. 

The  numerical  investigation  of  the  non-stationary  physico-chemical  processes 
development  regularities  after  the  boot  breaking-off,  was  executed  taking  into  account 
the  influence  of  a  number  of  the  factors;  the  boot  breaking-off  pressure  value,  the 
cavity  volume  which  opened  after  breaking-off,  the  length  on  which  there  is  the  boot 
breaking-off  and  the  duration  in  time  the  boot  breaking-off  period.  Established,  that 
the  non-stationary  processes,  caused  by  the  boot  breaking-off,  are  the  most  dangerous 
to  the  SPRM  functioning  at  instant  breaking-off  Conducted  numerical  investigations 
permit  to  make  a  conclusion  that  the  velocity  of  the  flame  propagation  along  the 
propellant  surface,  which  opened  after  the  boot  has  breaking-off,  weakly  depends  on 
the  boot  breaking-off  pressure  and  on  the  length,  on  which  the  breaking-off  occurs. 
However,  the  increase  of  the  above  mentioned  pressure  and  the  length  results  to 
increase  of  the  propellant  ignition  time  delay  in  the  opened  cavity. 

Realization  of  the  described  above,  or  similar,  ignition  anomalous  mode  in  the 
large-sized  solid-propellant  boosters  of  the  Reusable  Transport  Space  System  (RTSS) 
(fig.  5)  is  extremely  undesirable.  The  asymmetry  of  two  solid-propellant  boosters 
distribution  in  the  lay-out  diagram  put  in  the  essential  limitation  on  instability  level  of 
the  ignition  processes  and  subsequent  combustion  in  the  boosters  combustion 
chambers.  Realization  of  the  described  above  an  ignition  anomalous  mode  in  one  of 
the  boosters  can  lead  to  infringement  of  the  symmetry  of  the  RTSS  total  thrust 
distribution  and  appearance  of  the  destabilizing  moment.  And  it,  depending  on  the 
boosters  thrusts  misalignment  level,  can  provide  essential  influence  on  the  RTSS 
motion  stability  on  a  planned  trajectory,  to  create  nonexpected  actions  on  the  payload 
and  even  lead  to  loss  of  the  RTSS  controllability. 

Executed  calculations  demonstrate,  that  the  described  above  ignition 


Figures.  The  RTSS  lay-out  diagram. 

abnormal  mode  appearance  in  the  SPRM  it  is  possible  to  prevent  through  the 
monotony  of  the  engine  operation  during  the  start  period  providing,  by  means  of  the 
IS  CP  flows  selection  both  on  the  direction,  and  On  the  time.  It  can  be  realized  by 
means  of  using  of  the  IS  special  design  scheme,  which  can  provides  at  first  the  CP 
supply  only  in  the  SPRM  charge  channel,  In  this  case,  joint  elastic  deformation  both 
the  propellant  charge  and  the  SPRM  case  will  be  proceed.  And  it  lead  to  the  face 
clearances  volumes  increase  still  up  to  the  moment  of  the  flame  appearance  on  the 
propellant  charge  surface.  Then,  the  CP  supply  in  the  deformed  face  cavity  begins, 
that  will  be  the  favorable  factor  for  the  initial  stage  of  the  SPRM  operation.  For 
example,  in  this  conditions  the  9P/9t  maximum  values  are  reduced,  the  flame 
propagation  velocity  along  the  face  clearance  increase  in  some  times,  the  time 
straggling  of  the  propellant  surface,  located  in  the  face  cavity  connecting  to 
combustion  have  decreased,  and  the  boot  breaking-off  probability  is  essentially 
reduced. 
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THE  NUMERICAL  STUDY  OF  THE  POROUS  COOLING 
EFFECT  ON  THE  CALCULATION  EXAMPLE  OF  HEAT 
TRANSFER  IN  ZONE  OF  THE^^  SHOCK-WALL 
INTERACTIONS 

ISAEV  SA.,  SADOVNIKOVG.S.  (CIVIL  AVIATION  ACADEMY) 

LEONTIEV  AX,  NOSATOVVV  (MSmbyhimNEBawnami) 

I’he  numerical  study  of  injection  effect  on  thermal  output  at  the  hypersonic  flow 
of  viscous  gas  at  the  presence  of  heat  transfer  district  in  the  place  of  the  shock 
drop  has  target  the  substantiation  of  the  concept  of  porous  cooling  in  the  heat 
protection  systems  for  power  plants  of  modern  aircraft.  Elaborated  on  of  the 
solution  by  the  determining  method  of  discretizied  by  means  of  explicit-implicit 
MacCormack  scheme  Navier-Stokes  equations  for  the  laminar  regime  of  flow 
and  Reynolds  equations,  closed  by  means  of  two-parameter  dissipative  model  of 
turbulence,  program  package  [1  -  3  ]  is  applied  for  the  task  solution  about 
laminar  hypersonic  prewall  flow  of  viscous  gas  with  a  drop  shock.  Whereat  for 
the  constriction  of  the  calculating  region  sizes  it’s  accommodated  symmetry  of 
flow  concerning  top  controlling  plane,  parallelled  low  wall  and  coming  out  ^m 
sharp  of  the  shock  generator.  Thus,  it’s  considerred  the  development  of  the 
h>'personic  flow  of  viscous  gas  in  the  channel  of  variable  transversal  section  with 
convergent  middle  partly.  At  the  task  solution  it’s  neglected  possible  at  such 
velocities  of  flow  by  physics  and  chemical  conversions.  At  the  assignment  of 
initial  conditions  the  parameters  of  gas  into  all  calculated  grid  points  are 
admitted  to  the  equal  parameters  of  gas  in  nondisturbed  flow.  On  entrance  the 
parameters  of  flow  are  supported  continual  constant!  On  exit  ones  are  set  of  soft 
boundary  condition  or  the  conditions  of  the  solution  continuation.  On  lower 
CTeHKe  in  the  calculation  process  traditional  slip  conditions  are  fulfilled,  while 
the  wall  temperature  is  supported  constant.  On  tile  overhead  surfece  of  the 
shock  generator  conditions  corresponding  to  thermally  insulated  wall  are 
admitted. 

In  carried  out  analyses  calculated  region  broke  by  grid  with  the  bodying  of 
points  at  the  walls,  especially  at  lower  wall  and  in  the  field  of  shock  drop,  where 
recircuting  zone  is  formed.  Grid  contained  81x61  calculated  points.  Minimal 
step  at  lower  wall  was  chosen  approximately  lO"^.  It’s  necessery  to  notice  ,  that 
the  height  of  channel  section  on  region  exit  was  considered  as  an  characteristic 
linear  size,  and  velocity  and  the  density  of  gas  stream  are  chosen  as  an  scales 
variable.  Comparative  analysis  of  the  calculated  results  for  hypersonic  viscous 
flow  near  impenetrable  wall  and  wall  with  the  district  of  porous  cooling  at  the 
presence  of  falling  the  shock  is  carried  out  in  the  approach  of  the  laminar 
regime  of  flow  (Reynolds  number  is  chosen  equal  to  5000).  Mach  number  was 
assigned  8,2,  and  the  value  of  wall  temperature  and  input  through  the  porous 
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insertion  of  gas  was  admitted  equal  to  0.27.  The  district  of  porous  cooling  was 
situated  in  the  most  heat  transfer  place  (with  x=14  till  x=16)  corresponding  to 
place  of  the  shock  dropping  on  the  wall.  The  impulse  of  gas  injection  was  fixed 
and  was  chosen  equal  to  5%. 


On  the  picture  are  presented  the  results  of  the  comparative  analysis  of 
distributions  along  the  wall  relative  pressure  P/Pinf  (a)  and  relative  heat  current 
Q/Qinf  (b)  in  the  various  moments  of  time  (1  -  T==3;  2  - 10;  3  -  15;  4  -  20;  5  -  20) 
without  injection  (1-4)  and  with  (5).  Here  index  «inf>>  characteri2es  the  gas 
parameter  for  hypersonic  flow  along  flat  wall.  The  analysis  of  the  characteristics 
field  is  carried  out  at  the  achievement  of  the  solution  determining,  i.e.  at  the 
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achievement  of  the  steady  regime  of  the  gas  flow.  This  state  is  achived  at  times 
over  15.  The  pattern  of  the  stationary  viscous  prewall  flow  discovers  the 
decelerating  zone  in  the  most  heating  place  -  the  place  of  the  shock  dropping. 
I'he  recirculating  zone  is  created  in  pointed  place,  although  in  the  case  without 
injection  it’s  e?q)ressed  rather  weakly.  The  presence  of  cooled  insertion  does  not 
accord,  just  like  followed  to  expect,  of  essential  influence  on  the  flow  field.  This 
is  acknowledges  comparative  analysis  of  the  graphic  distributions  of  pressure 
and  temperatures  in  the  prewall  fields  at  presence  and  absence  injection. 

Transformation  with  the  time  of  the  superficial  distributions  of  the  relative 
values  of  static  pressure  and  heat  current  demonstrates  the  process  determining 
at  prewall  flow  without  injection  with  formation  of  thermal  stress  zones 
corresponding  to  place  the  shock-boundary  layer  interactions.  In  this  case 
interestingly  to  markj,  that  flow  pattern  and  distributions  of  local  force  and  heat 
load  generally  are  formed  to  the  moment  of  time  T=10.  Injection  does  practically 
not  accord  sizeable  influence  upon  the  profile  of  relative  pressure  on  the  wall. 
I’here  is  noticeably  only  some  increasing  of  pressure  in  the  place  of  the 
beginning  of  gas  injection.  The  analysis  of  calculated  heat  current  on  the  wall 
displays  that  realization  of  injection  through  porous  insertion  lets  essentially  to 
take  offbeat  loading  on  the  wall,  what  stimulates  the  effectiveness  of  this  type  of 
cooling.  This  conclusion  is  confirmed  in  physical  experiments  [  4  ].  Thus,  in 
carried  out  numerical  analyses  is  substantiated  the  concept  of  the  porous  cooling 
of  thermal  stress  places  of  the  flow  around  elements  of  the  objects  of  various 
destination.  Whereat  is  marked  local  character  of  the  actions  of  this  type  of 
cooling,  not  bringing  to  essential  influence  on  the  pattern  of  flow  around  the 
object. 

This  work  is  carried  at  financial  support  of  RFFI  under  the  project  N96-01- 
00298. 
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MAGNETIC  FIELD  EFFECT  ON  HEAT  EXCHANGE 
AT  THE  BURNING  SURFACE 

/.  G.  BOROVSKOY,  A.B.  VOROZHTSOV 

RESEARCH  INSTITUTE  OF  APPLIED  MATHEMATICS  AND 
MECHANICS,  TOMSK,  RUSSIA 

A  laminar  flow  of  a  noncompressible  conducting  medium  in  a  plane 
channel  of  height  2b  exposed  to  a  transverse  uniform  magnetic  field  is 
considered.  The  bottom  and  top  walls  of  the  channel  are  the  burning  surf  aces  of 
a  solid  propellant,  the  conducting  lateral  walls  being  shorted.  The  variation  of 
channel  height  due  to 
propellant  consumption  is 
ignored.  The  magnetic  flux 
density  vector  in  the 
Cartesian  coordinates  is 
B  ==  (0,  0)  and  has  one 

component  in  the  flow  zone, 

By -const.  'I'hc  value  of 
electric  current  density  j, 
according  to  generalized  Ohm's  law,  also  has  one  component  with 

electric  currents  shorting  to  the  lateral  walls  of  the  channel  and  the 
ponderomotive  force  vector  opposing  the  flow  f  =  0, 0),  -  -j^By. 

U  follows  from  the  solution  of  Hartman's  problem  for  the  considered  flow 
in  the  absence  of  heat  and  mass  exchange  with  the  channel  walls  [1]  that  with 
other  conditions  being  equal  the  increased  magnetic  induction  of  the  external 
field  leads  to  increased  friction  stress  near  the  wall.  Accepting  the  Reynolds 
assumption  that  the  dynamic  and  thermal  fields  are  similar  for  the  given 
problem  [2],  one  should  expect  certain  intensification  of  the  heat  transfer  to  the 
channel  walls.  This  in  turn  can  increase  the  surface  temperature  and  burning  rate 
since  usually  for  solid  propellants  valid  is  the  exponential  dependence  of  the 
mass  burning  rate  m  on  the  siuface  temperature  0^,  at  a  constant  pressure  [3] 

m(e^)  =  ~  (1) 

where  =  7^/7^  is  the  ratio  of  the  propellant  surface  temperature  to  the  gas 
temperature  7^  in  the  flow  core,  and  are  constants. 

The  approximate  value  of  magnetogasodynamic  effect  on  burning  surface 
temperature  is  estimated  from  the  heat  balance 


(2) 
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Dividing  both  sides  of  Eq.  (2)  by  Tg  and  taking  into  account  the  Stanton 
number  definition 


where  Wg  is  the  mass  gas  velocity  in  the  flow  core,  after  transformations  we  have 

Si  ^  exp(-ei,  I Q,).  0) 

e,-e„  A.  w 

Equation  (3)  applies  both  to  the  case  with  applied  magnetic  field  (in  this 
case  the  parameters  are  marked  with  the  index  H)  and  to  the  case  with  the 
magnetic  field  absent  (index  o).  To  estimate  the  Stanton  number  in  the  presence 
of  magnetic  field  we  make  use  of  the  solution  of  Hartman’s  problem 


MiHa)  = 


W  thilia) 

3  Ha-ikHaY 


where  th(Ha)  is  the  hyperbolic  tangent  of  the  Hartman  number  Ha, 

’I'hen,  substituting  the  Reynolds  analogy,  which  in  terms  of  Eq.  (4)  has  the 
form  SIj^ISIq  =  into  Bq.  (3)  we  have  the  final  equation 


St',  w{Ha)  =  exp(-e„IQ,), 


where  the  Stanton  number  St^  constitutes  the  heat  exchange  coefficient  in  the 
absence  ol'  magnetic  field,  the  constant  rn^,  and  the  specific  heat  coefficient  of 
solid  propellant  This  substitution  is  possible  if  the  magnetic-field-induced 
changes  in  the  mass  velocity  of  combustion  products  in  the  flow  core  (unlike  the 
mass  binning  rate)  are  neglected. 

Figures  1  show'  the  numerical  solution  of  nonlinear  equation  (5)  at  the 
dimensionless  parameter  ranging  within  e  [0,  3].  The  ratios  0^.  and 
are  independent  of  the  Stanton  number  SVq  and  weakly  depend  on  the 
temperature  The  results  obtained  suggest  that  with  increasing  Hartman’s 

number  the  propellant  surface  temperature  increment  is  15  %  (at  Ha  <  20),  while 
the  burning  rate  shows  more  than  two-fold  increase,  which  is  accounted  for  by 
the  exponential  dependence  (I ). 
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The  results  presented  are  qualitative  and  are  indicative  of  possible  essential 
influence  of  external  magnetic  field  on  the  combustion  of  solid  propellants.  The 
authors  did  not  extend  the 
range  of  Ha  variation  at  4 

rn^lniQ »  1  since  the 

intense  efflux  from  channel 
walls  would  violate  the  3 

assumptions  accepted  in  the 
approximation  of  both  ^ 

laminar  and  turbulent 
boundary  layers  [2].  Hence,  ^ 

Eq.  (4)  used  as  the  basis  for 
the  present  analysis,  could 
also  be  violated.  Thus,  a  1 

detailed  investigation  of  the  ^  5  10  15  20 

problem  can  be  performed  Ha 

and  quantitative  results  can  Fig.  1.  Plots  for  dimensionless  mass  burning  rate 

be  obtained  only  by  the  versus //a  and  l-c^=0;  3-1.0;  5-2.0. 

solution  of  the  complete 

Navier-Stokes  equations  for  a  compressible  gas,  supplemented  by  equations  for 
diffusion  and  magnetic  induction  as  well  as  by  the  right  sides  responsible  for  the 
chemical  reactions  proceeding  in  the  gas  flow  and  at  the  burning  surface. 
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THE  EFFECT  OF  EXTERNAL  ELECTRIC  FIELD  ON  HEAT  AND  MASS  TRANSFER 
CHARACTERISTICS  IN 

THE  POWER  FACILITY  COMBUSTION  CHAMBER 

YAGODNIKOV  D.A.,  VORONETSKII  A.V.,  PELEVIN  F.  V. 

Baiiman  Moscow  State  Technical  University 

Experimental  study  results  of  the  external  electrical  field  (EEF)  effect  on  operating  process 
characteristics  and  heat  and  mass  transfer  peculiarities  in  a  model  porous-wall  combustion  chamber  (CC) 
are  given.  It  is  found,  that  when  a  900  V  constant  negative  potential  is  applied  to  the  CC  the  combustion 
products  heat  flow  to  CC  w'all  is  increased,  and  in  case  of  a  positive  potential  the  heat  flow  is  decreased. 
In  both  case,  it  was  been  established  that  ignition  and  combustion  processes  is  intensified  and  blow-off 
velocity  grows.  The  data  obtained  is  explained  by  means  of  ionic  wind  coinciding  in  direction  with 
positive  ions  movement.  With  negative  voltage  CC  the  ionic  wind  directed  toward  to  chamber  wall 
interacts  with  the  wall  boundary  layer  and  determines  a  modification  of  its  geometrical  and 
thermophysical  parameters. 
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MATHEMATICAL  MODELLING  OF  TRANSITIONAL 
REGIMES  FOR  CONVERSION  OF  HIGH  ASH 
BITUMINOUS  COAL  IN  THE  PRESSURISED 
CIRCULATING  FLUIDIZED  BED  FACILITY 

MAY  ST  REN  KO  A.  Yu.,  PATSKOV  V.P.,  TOPALA.I. 
Coal  Energy  Technology  Center  (CETC) 

National  Academy  of  Sciences  of  Ukraine 
MINENERGO  of  Ukraine 


The  conversion  (pyrolysis,  combustion  and  gasification)  of  low 
grade  high  ash  coals  in  the  combined  cycle  units  using  pressurised 
circulating  fluidized  bed  combustion  (PCFBC)  is  one  of  the  promising 
ways  to  increase  the  efficiency  of  combustion  while  decreasing  NO^  and 
SOx  emissions  [  1  ] . 

As  a  part  of  activity  aimed  at  implementation  of  new  clean  coal 
technologies,  the  Coal  Energy  Technology  Center  (CETC)  is  developing 
pilot- scale  facility  with  coal  consumption  of  1  ton/ hr  to  carry  out 
investigation  on  conversion  of  bituminous  coals  at  the  PCFBC  under 
pressure  of  up  to  2.5  MPa  and  with  high  number  of  recirculation  (more 
than  100)  [IJ.  The  main  components  of  the  facility  being  developed  are: 
i)  pyrolyser  with  slowly  movable  bed  of  fuel  at  the  return  circuit  of 
coke  residue,  2)  reactor,  hydrodynaniically  separated  at  fluidized  bed 
zone  and  freeboard  space,  to  produce  syn-gas,  3)  cyclone-separator. 

The  mathematical  models,  algorithms  and  program  for  numerical 
analysis  of  transitional  regimes  of  conversion  process  for  polyfractional 
high  ash  coal  in  the  above  PCFBC  facility  are  considered. 

The  mathematical  model  of  pyrolyser  is  made  using  the  fixed  bed 
assumption.  The  interaction  of  two  phases  are  considered,  namely:  solid 
phase,  including  porous  coal  particles,  coke  residue  and  ash,  which  are 
entered  from  above  and  gas  phase,  containing  components  of  fuel 
pyrolysis,  which  are  passed  through  bed  from  below.  The  fuel  pyrolysis 
components  consist  of  light  gaseous  products  (CO2,  CO,  H2O,  H2,  N2, 
CH4,  C2H4,  C2H6)  and  tar  vapour  as  well  as  hydrocarbons  of  phase 
(evaporation-condensation)  and  chemical  reactions  inside  the  porous 
structure  of  coal  particles  [2].  The  tar  composition  is  modelled  by  seven 
key  components:  benzol,  toluol,  phenol,  phenatren,  undekan, 
dimetilnaphthaline  and  metiltetralin  [3].  The  mass  and  energy  transfer 
over  the  bed  is  due  to  effective  diffusion  and  heat  conductivity,  filtration 
and  convective  heat  / mass  transfer  between  fuel  particles  and  ash  inside 
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the  free  space  of  the  bed.  The  effective  diffusion  and  heat  conductivity. 
cou^^ecti\-e  and  filtration  transfer,  primary  extraction  of  vapour  and 
volatile,  tar  and  hydrocarbons  phase  transition  (evaporation- 
condensation),  secondary  heterogeneous  gasification  reactions  for  coke 
residue  with  pyrolysis  components  are  considered  for  porous  coal 
particles.  It  is  also  taken  into  account  the  heating  up  of  inert  ash 
pa ] tides  of  bed. 

The  main  assumptions  made  for  the  fluidized  bed  zone  of  the 
reactor  producing  syn-gas  are  as  follows:  non-stationary  diffusive  model 
that  accounts  for  longitudinal  effective  diffusion  and  heat  conductivity 
for  solid  phase  transfer,  gas  component  filtration,  pha.se-to-phase 
heat.,  mass  transfer,  convective  and  radiative  heat,,  mass  transfer 
between  fuel  and  inert  ash  particles  and  gas  fiov%  near-  wall  heat 
transfer,  surface  heterogeneous  reactions  of  combustion /'gasification  of 
coke  residue,  homogeneous  reaction  tor  CO  oxidation. 

The  transfer  processes  for  bubble  phase  are  described  using  non- 
stationary  model  of  ideal  forcing  out  that  takes  into  account  phase-to- 
phase  heat,/  mass  transfer  and  homogeneous  reaction  of  CO  oxidation. 

The  surface  concentrations  of  gas  components,  temperature  of  coke 
and  inert  ash  particles  are  calculated  taking  non-stationaiy  model  of 
ideal  mixing  that  takes  into  account  convective  heat,  'mass  transfer 
between  particles  and  flow  in  gas  phase,  radiative  heat  transfer  between 
particles  and  gas,  surface  reactions  of  fuel  combirstion. /gasification. 

The  size  distribution  functions  for  solid  particles  are  determined 
basing  on  non-stationaiy  diffusive  model  that  accounts  for  effective 
longitudinal  diffusive  and  circulation  of  particles  as  well  their  change 
due  to  surface  chemical  conversion  and  carry  over  [4]. 

In  the  freeboard  space  of  the  reactor  producing  syn-gas,  the 
upstream  movement  of  polyfractionai  fuel  coke  surrounded  b}^  kigh 
temperature  gas  flow  in  the  entrained  flow  regime  is  considered.  The 
piocess  models  are  ba.sed  on  theory  ot  multi-rate  penetrating  phases  and 
turbulent  heterogeneous  streams  [5]  in  one  and  two  dimension  non- 
stationan,)'  statements.  The  surface  heterogeneous  reactions  of  combustion 
and  gasification  of  coke  residue,  homogeneous  reaction  of  CO  oxidation, 
convective  heat  /  mass  transfer  between  phases,  convective  and  radiative 
heat  transfer  between  particles,  flow  and  facility  walls,  turbulent 
transfer,  gravitation  forces,  pressure  drop,  phase-to-phase  friction, 
heating  up  of  reactor  walls  are  taken  into  account. 

The  models  of  transfer  processes  in  the  cyclone  are  based  on 
Lagiange  trajectory  theory  of  particle  mo^^emeiit  inside  the  twisted  gas 
flow  and  Eiler  assumption  on  gas. 
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The  balance  boundary  conditions  which  bind  proposed  assumptions 
into  the  whoJe  mathematical  description  of  the  process  in  the  facility  aie 
formulated. 

The  iterative  algorithms  of  process  calculation  (like  Zeidel  method) 
for  each  components  of  the  facility  are  developed. 

The  solution  of  the  equations  describing  the  pyrolyser  process  is 
provided  by  double-sweep  method.  An  special  feature  of  the  algorithm  is 
an  additional  assumption  for  the  system  about  thermodynamic  balance 
between  gas  and  condensed  phase  inside  the  pores  of  coal  particles,  that 
is  necessary  for  iterative  estimation  of  phase  transition  rates  of  liquid 
tars  and  hydrocarbons  [2]. 

The  systems  of  equations  for  solid  phase  of  fluidized  bed  and  for 
size  distribution  functions  are  solved  by  double-sweep  methods,  as  to 
bubble  phase  and  single  particles  -  like  relaxation  method  for  solving  of 
systems  of  hyperbolic  differential  equations. 

The  systems  of  equations  for  the  freeboard  space  and  the  cyclone, 
when  one  dimension  statement  is  used,  are  solved  the  same  manner.  To 
solve  two  dimension  tasks,  the  longitudinal-cross  double-sweep  method  is 
used  [5], 

The  method  of  automatic  choice  of  time  step  provided  b\' 
Richardson  extrapolation  is  used  for  each  algorithms. 

'I'he  ])rograi.n  to  calculate  process  parameters  for  the  faciiitv  is  made 
using  Fortran -77  (v.1.0  for  Wiindows) 
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GAS  DYNAMIC  EXPERIMENTAL  DATA  BASE 
OF  SURFACE  GROUND  STRENGTHS 

Vasin  V.A.,  Sannikav  K.A. 

198005,  Russia,  Si.-Peiersburg,  Bailie  Slate  Technical  Umversiiy, 

The  data  accrued  in  d^a  base  (DB)  are  used  as  training  and  examimr^ 
sampling  at  execution  prognosticate  of  accounts  on  erosion  of  surface  ground 
strengths.The  necessary  volumes  of  experimental  data  are  accumulated  by 
statement  of  special  laboratory  experiments  on  physical  models.The  unique  data 
on  erosion  of  surface  ground  strengths  in  conditions  are  used  as  the  basic  points 
for  checks  of  accuracy  of  statistical  models,  constructed  on  results  laboratory 
model  of  tests.  The  opportunity  of  direct  valuation  of  parameters  crater  by  direct 
carry  given  model  of  experiment  on  nature  the  object  is  provided. 

The  gas  dynamic  experimental  d^a  base  of  surface  ground  strengths 
represents  the  managing  programme  and  information  environment,  ensuring,  in 
regime  of  interactive  dialogue  with  personal  computer  ( PC  ): 

-  management  and  accumulation  of  base  of  experimental  data: 

-  fast  search  of  necessary  information; 

-  updating  of  structure  DB; 

-  initial  formation  for  ccn'elation,  regression  and  variance  of  analysis; 

-  prognosis  the  accounts  of  sizes  crater; 

-  documenting  of  results. 

The  data  composition. 

Accrued  in  supervision  bank  for  erosion  by  processes  in  surface  ground 
strengths  the  information  is  directed  on  use  of  her  in  quality  of  training  sample 
of  statistical  models  given  at  set-up  for  forecasting  of  sizes  forming  in.  ground 
crater.  The  statistical  uniformity  is  reached  criterion  by  foim  of  representation  of 
data.  Quality  of  determining  criteria  of  similarity,  forming  structure  of  data, 
nominates  the  criteria,  pursuant  to  which  the  gears  of  formation  crater  act.  The 
most  representative  by  criteria  of  similarity  for  processes  considered  here  the 
numbers  Maha  and  Fruda  are  considered.  Because  of  inconsistency  of 
requirements  to  structures  size  invariant  of  similarity,  propound  by  these  criteria, 
kind  of  data,  accommodate  in  bank,  is  defined  alternatively  on  some  of  one  of 
mentioned  criteria. 

At  use  in  quality  of  determining  criterion  of  similarity  of  number  Maha, 
table  of  supervision  is  made  following  imsiring  parameters  and  rizing  invariant  of 
similarity; 
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Characteristic  gas  dynamc  level  ofIoa± 

-  impulse,  on  unit  of  working  volume  nozzie  J'  -J  / 

H  s/m^\  -  average  traction,  on  unit  of  working  area  of  target  section  traction 
P^Pfd^,NAn\  -  traction  armament  of  product  T-gm,^p\  -  the  relative  distance 
from  null  levd  traction  to  ground>''=y/rf;  -  unsizing.  the  time  of  effect  of  jet  on 
ground  ( number  Struha!  )t'~t  -  v/y; 

The  physical-mechamcal  properties  of  ground; 

-  indication  dynamic  densitometer  Cy^;,  quantity  of  impacts  /  the  depth  of 
immersing;  -  humidity  w,  %;  -  pore  e;  -  density  p;  -  pressure  a,  -  module  £, 

nmS 

The  sizes  crater: 

~  relative  diameter  D  /  d;  -  relative  depth  H  /  d. 

At  use  in  quality  of  determining  criterion  of  similarity  of  number  Fruda  in 
table  of  supervision  unsizing  the  parameters  remain  the  same,  that  and  in 
previous  item.  Sizing,  invariant  of  similarity  present  to  base  object  on  formulas; 

/  =  E'  =  E‘Ki 

Where  scale  factor  (  factor  of  reduction  );  - 

diameters  of  target  section  of  nozzle  of  base  and  obser^^able  objects;  f,P\cr\E'  - 
adduced  pecific  pulse,  specific  traction,  pressure  of  ground  and  its  module. 

The  base  object  nominates  the  most  the  characteristic  object  from 
considered  class.  At  placing  in  bank  the  data  informations,  received  by 
laboratory  tests  of  physical  models,  as  such  object  is  used  prototype  of  model. 

The"  data  are  adduced  in  natural  kind.  Transformed  to  mentioned  above 
by  way  in  criterion  form,  these  data  are  united  with  other  data,  received  in  other 
operational  conditions.  Thus  the  statistical  uniformity  of  incorporated  data  is 
preserved,  and  training  sample  becomes  more  representative  one. 

The  increase  of  reliability  prognosis  valuations  of  sizes  crater  at 
simultaneous  decrease  of  operational  restrictions  on  possible  area  their 
distributions  can  be  achieved  by^  updating  DB  by  new^  data,  received  on  ground 
platforms  with  other  properties,  at  other  constructive  -  composite  and  traction 
characteristics  booster. 

The  structure  of  DB. 

DB  is  executed  in  kind  of  user  interface,  is  constructed  on  system  of 
multilevel  menu.  DB  is  written  on  language  Clippers.  1.  The  data  are  placed  on 
groups: 

-  supervision  for  model  by  stait-up,  classified  on  types  of  ground,  criteria, 
scales  of  simulation  and  prototypiein  natural  measurement  of  parameters; 
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-  supeivision  foi  natural  by  start-up,  classified  on  objects  in  natural 
tneasurement  of  parameters; 

-  supervision  for  stan-up  in  criterion  kind. 
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FigJ  Fig, .2 

On  fig.  1  and  2  the  results  of  statistical  processing  of  tests  on  gas  t^asnic 
base  of  surface  gritty  strengths  are  submitted.  The  volume  of  tests  has  made  about 
40  experience.  The  analysis  of  accuracy  and  adequacy  of  model  is  conducted. 
The  schedule  on  fig.  1  is  constructed  on  “paper  probability  value”,  ropriate  to 
normal  law  of  distribution.  The  remainders  arrange  on  direct  line,  that  peaks 
about  weak  interrelation  between  remamders,  i.e.  about  their  casual  character. 
Fig.  2  illustrates  the  examples  of  arrangement  settlement  craters  on  attitude  to 
uiure  to  objects  for  gritty  ground  of  various  humidity,  pore  and  density.  The 
analysis  of  accuracy  and  adequacy  of  statistical  model  is  conducted.  On  fig.  2  by 
vertical  strokes  the  95  %  the  confidence  intervals  for  settlement  significances  of 
diameters  a'^ers  are  shown. 


Section  3 

TWO-PHASE  FLOWS 
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NUMERICAL  INVESTIGATION 
OF  TURBULENT  TWO-PHASE  FLOW 
IN  THE  CHANNEL  WITH  POROUS  WALLS 

EMELYANOV  V.N.,  VOLKOV  K.N, 

Baltic  State  Technical  University 

The  sectionf?  of  gaseous  highways  with  penetrable  walls  are  one  of  most 
encountered  hydrodynamic  compound  units  of  thermal  and  power  devices. 
The  flow  in  the  channel  with  distributed  gas  blowing  is  a  model  of  flow  solid 
propellant  combustion  products  in  the  gap  maked  up  wall  of  combustion 
chamber  and  burning  surface  cartridge. 

Turbulent  flow  with  the  suspended  particles  has  been  simulated  on 
the  basis  of  averaged  equation  systems  of  interpenetrating  continuum 
mechanics,  turbulent  viscosity  conception  and  equations  of  two-parameter 
dissipative  turbulence  model.  Correlative  moments  and  additional  source 
terms  £p  and  associated  with  dispersed  phase  have  been  calculated 
out  of  spatial  and  temporal  averaging  method  [1]  in  local  and  uniform 
approximation.  Iterphase  iteraction  have  been  defined  with  the  drag  force 
according  to  Stokes  formula.  Admixture  particles  are  rigid  undeformed 
spheres  of  equal  diameters.  Thermal  and  physical  properties  of  phases 
have  been  assumed  permanent  values.  The  frame  of  reference  has  selected 
so  that  axis  x  coinsides  with  the  channel  centerline  and  axis  y  has 
directed  perpendicularly  it.  Positive  reference  of  lateral  component  has  been 
conducted  to  the  side  contrary  of  blowing  direction. 

The  simulation  of  two-phase  turbulent  flow  in  the  channel  with  porous 
walls  has  been  developed  in  the  framework  of  calculation  general  conception 
of  accelerating  gas-particle  flows  proposed  by  the  authors  [1.  2].  Let  us 
introduce  dimensionless  coordinates  ^  x  /  Vyj  ,  rj  ■=  y  f  and  shall  seek 
the  solution  in  the  form 

PgUgy^  =  PwVy,  ^  Ug{ri)  r^,  pgVgy^  =  -p,oVy,  Vg{r])  r^,  dg  = 

Pp'^pV^  ~  Piu'^w^'Kpiji)  Pp'^’py^  ~  ~  Pw'^wVp{Tj)  'dp  —  'dy^Tpi^rj^, 


k  =  vlK{r]),  £  =  vlE{r})/ry,, 

where  r^,  —  half  width  of  channel,  Vy,  —  blowing  rate,  x  =  0  foi"  planar 
flow  and  X  =  1  for  axially  symmetric  flow. 
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Then,  neglecting  with  the  terms  of  multiplicity  of  a  ,  the  turbulent 
two-phase  flow  calculation  in  the  duct  with  distributed  gas  blowing  has 
been  reduced  to  the  solution  of  boundary  problem  for  nonlinear  ordinary 
differential  equations 


U,  ~  V!  =  0; 


{U,U,)' -  {Vfi-',)' ^ 


{U,T,)  -  iV,T,)'  = 


1  Cn  K  1 

- 1 — - T 

Pr  E  '  “ 


+  f2i^p'0  {Tg  —  Tp) ; 


v;=o-, 


{U,U,)  -  (V,U,)'  =  /a  {U,  -  U,)  +  ; 

(V,T,)  -  iV,T,y  =  h  (T,  -  T,)  +  . 


{U,K,)  -  (V,K)'  = 


Ok  E  I 


+  P-E-E,: 


(U,B)  -  iV,E}'  = 


(± 

+  Sl^i 

\Re 

+  ^  {CeiP  -  caE)  - 


where 


/i  =  l/Ski,.  /2  =  1/Skr.  E,  =  re,lvl,  , 

Mathematical  formulation  of  the  problem  has  been  completed  with  the 
representation  of  boundary  conditions  at  the  centerline  and  channel  walls. 
The  boundary  conditions  for  averaged  flow  properties  under  the  rj  =  1 
have  been  determined  by  means  of  assignment  of  wall’s  temperature,  no-slip 
condition  for  longitudinal  component  velosity  and  normal  blowing  predicate 
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to  the  channel  walls  for  lateral  velosity  component  of  the  carry  phase.  The 
symmetry  flow  conditions  have  been  supposed  under  the  rj  =  0  . 

The  obtained  distributions  of  turbulence  kinetic  energy  and  energy 
dissipation  rate  have  been  used  for  statistical  simulation  of  particles  motion 
trajectories  in  the  stochastic  and  inhomogeneous  medium. 


Distributions  of  turbulence  kinetic  energy  (1)  and  energy  dissipation  rate  (2)  at  the 
various  cross-section  of  the  channel  with  porous  walls.  The  lines  correspond  to  the 
following  values  ^  =  5  (1);  10  (2);  15  (3);  20  (4);  25  (5). 


Turbulent  structure  formation  and  particle-phase  flow  patterns  in 
channels  are  studied  depending  on  size  and  volume  fraction  of  particles. 
The  obtained  results  are  most  efficiently  can  be  used  for  solution  of  applied 
problems  and  for  recommendation  establishment  according  to  practical 
implementation  of  optimal  operating  conditions  of  intra-chamber  processes 
in  the  rocket  solid  propellant  motors  and  power  devices.  Developed  models 
are  the  base  for  investigations  of  turbulent  transfer  processes  in  the  domains 
with  more  complicated  geometry. 
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TWO-PHASE  GAS-SOLID  FLOW  AROUND  A  CYLINDER: 
EFFECTS  OF  PARTICLE  SIZE  DISTRIBUTION 
AND  PARTICLE-PARTICLE  COLLISIONS 

A.K  VOLKOV,  YU.M.  TSIRKUNOV 
Baltic  State  Technical  University,  198005  St.  Petersburg,  Russia 


In  computational  models  of  gas-solid  particle  flows  particles  are  usually 
assumed  to  have  an  equal  size.  At  the  same  time  there  is  always  a  particle 
size  distribution,  and  a  particle  radius  can  vary  in  a  rather  wide  range. 

The  kinetic  model  and  the  DSMC  method  were  proposed  in  papers  [1,2] 
for  the  simulation  of  mono-sized  particle  phase  motion  in  transient  flow  regime 
on  the  Knudsen  number  in  the  pseudogas  of  particles.  The  DSMC  method 
is  similar  to  the  one  in  rarefied  gas  dynamics.  It  can  be  extended  rather 
simply  to  the  case  of  different  particle  sizes  in  the  dispersed  medium.  The 
aim  of  the  present  work  is  to  study  how  the  particle  size  distribution  influence 
the  dispersed  phase  flow  field  and  the  action  of  particles  on  the  streamlined 
surface.  The  cross-flow  of  a  dilute  gas-particle  mixture  around  a  circular 
cylinder  is  considered  as  an  example. 


The  particle  size  distribution  can  be  described  by  the  particle  concentra¬ 
tion  distribution  function  in  a  particle  radius  fp.  It  is  often  convenient  to 
deal  with  the  particle  volume  concentration  (Q;p5°(rp)drp  is  the  volume 
concentration  of  particles  with  a  radius  from  to  rp-f  drp)  [3].  For  the  DSMC 
simulation  the  distribution  function  at  the  entrance  boundary  of  the 

calculational  domain  must  be  given  (the  entrance  boundary  is  located  in  an 
undisturbed  flow).  A  random  radius  rp  of  each  new  particle  flying  into  the 
calculational  domain  is  sampled  from  the  distribution  function  /^(rp)  of  the 


number  particle  density  =  r  ^gLi^p)/  J  r^^gl^{rp)drp). 

particle  phase  macroparameters,  such  as  the  volume  concentration  a-^ 


The 
in  an 


arbitrary  point  in  the  flow  and  the  energy  flux  from  particle  phase  to  a 


body  surface  are  determined  from  the  integrals  taken  over  translational  and 


rotational  velocities  and  over  radiuses  rp  of  particles 


J  m^f{r,Zp)dZp,  E^{ry,)  =  j  AEp\v^  ■n,y\f{r^,z^)dz^.  (1) 


Here  r  is  a  radius-vector  of  a  point  in  the  physical  space,  rw  a  radius-vector  of 
a  point  on  the  body  surface,  p®  and  mp  the  material  density  and  the  mass  of  a 
particle;  Zp  =  (?/p,rp),  =  (up,ajp);  f{r,Zp)  is  the  distribution  of  particles 

in  velocities,  radiuses  and  locations  in  the  physical  space  (J  f{r,Zp)dzp  — 
np{r)  is  the  number  particle  density  at  the  point  with  the  radius-vector  r, 
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//(^!-Sp)dyp  =  np(r)/°(r,rp)),  nw  is  the  outer  surface  unit  normal  at  the 
point  with  the  radius-vector  r^;  =  E^[z~)  -  Ep  is  the  kinetic 

energy  of  a  particle,  z~  vl  z'^  =  2+  {z~ ,  n^)  are  parameters  of  a  particle  before 
and  after  its  rebound  from  a  surface  at  the  point  with  the  radius-vector 
Ew  in  (1)  is  one  of  the  key  parameters  in  the  process  of  abrasive  erosion  of  a 
body  in  the  gas-particle  flow. 

In  the  present  study  the  log-normal  law  is  accepted  for  the  particle 
distribution  in  an  undisturbed  flow  [3] 

Inrp  -Inr^^yl 

^/21na  ;  J  ■ 

When  calculating  the  motion  of  particles,  the  aerodynamic  force,  the 
Magnus  force  and  the  damping  torque  are  taken  into  account  as  it  was 
done  in  [1,  2].  The  influence  of  the  particle  phase  on  the  carrier  gas  flow 
is  assumed  to  be  negligible.  The  semiempirical  impact  model  proposed 
in  [4]  is  used  to  determine  parameters  of  a  particle  just  after  its  rebound 
from  a  body  surface.  Hydrodynamic  interaction  between  particles  are  not 
taken  into  account.  Velocities  of  two  particles  after  their  binary  colhsion  are 
calculated  from  the  classical  impact  theory  for  undeformable  rotating  spheres 
with  the  use  of  two  restitution  ratios  and  art  for  the  normal  and  tangential 
components  of  the  contact  point  velocity  [5].  Kinetic  energy  losses  of  colliding 
particles  depend  on  the  values  of  am  and  art- 

Motion  of  particles  of  synthetic  corundum  in  the  potential  incompressible 
flow  of  the  carrier  gas  (air)  around  a  cylinder  of  mild  steel  is  studied 
numerically.  The  governing  parameters  have  the  following  values:  am  =  0.5, 
art  -  0.9,  pI  =  3950  kg/m^,  =  100  m/s,  =  1.197  kg/m^  = 
1.71  •  10“®  kg/(m-s),  radius  of  a  cylinder  has  been  taken  equal  to  1  m.  The 
particle  volume  concentration  apoo,  the  most  probable  particle  radius  in  an 
undisturbed  flow  fp  =  exp  [—  In^  a]  and  the  parameter  a  in  the  formula 
(2)  are  varied.  The  last  parameter  describes  the  particle  size  dispersion  in  an 
undisturbed  flow  D?-p  =  exp[31n^  a]  (exp[ln^  cr]  -  l)- 

In  the  numerical  simulations  the  truncated  distribution  function  p^(rp) 
approximating  9^{rp)  with  high  accuracy  is  used  (^^(rp)  = 
where  C  =  £?^(rp)drp,  rp  G  [rmin,  rmax])-  Boundary  values  and 

Tmin  liave  been  determined  from  the  conditions:  p^(rniax)  and  < 

O-Olp^(fp),  C  >  0.99.  In  all  calculations  r^ax  and  Train  satisfy  the  following 
inequalities:  r^in  >  1  pm  h  r^ax  <  250  pm. 

The  scheme  of  the  algorithm  and  the  parameters  of  the  numerical  method 
(size  of  a  mesh  cell,  a  number  of  modelling  particles,  etc.)  are  given  in  [2]. 

The  particle  volume  concentration  distribution  along  the  axis  of  symmet¬ 
ry  depending  on  cr  is  shown  in  figure  !,«.  These  results  refer  to  the  case  of 


5^(^p)  = 


1 


ln<T 


exp 
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infinitely  small  particle  concentration  in  the  undisturbed  flow  when  particle'- 
particle  collisions  are  negligible.  One  can  see  that  the  distribution  changes 
strongly  when  particle  size  dispersion  increases  {a  changes  from  1  to  1.2^). 
Peaks  usually  observed  in  the  case  of  mono-sized  particles,  first  become 
smooth  and  then  disappear  at  all.  Conversly,  the  particle  size,  dispersion 
plays  no  role  in  distribution  of  macroparameters  (1)  along  the  body  shape 
(curves  1-3  in  the  figure  2  practically  coincide). 


Figure  1.  Distribution  of  the  relative  particle  concentration  Qp/op,:,©  along  stagnation 
streamline  ^  =  0  in  the  flow  around  a  cylinder,  fp  =  100  fim;  (a)  collisionless  particle 
phase  motion,  (6)  particle  phase  motion  with  particle-particle  collisions.  Curve  1,  cr  =  1 
(monodisperse  particles);  2,  <7  =  1.2;  3,  cr  =  1.2^;  4,  cr  =  1.2®. 


Opposite  situation  takes  place  when  particle-particle  collisions  are 
essential.  In  this  case  the  particle  flow  field  structure  does 'not  change 
qualitatively  when  the  particle  size  dispersion  is  varied  in  the  wide  range 
(see  fig.  1,6).  But  the  particle  concentration  distribution  along  the  body 
shape  and  the  distribution  of  the  dispersed  phase  energy  flux  to  the  surface 
change  noticeably  (see  curves  4-6  in  the  figure  2).  It  should  be  noted  that  the 
role  of  particle-particle  collisions  in  forming  of  particle-phase  flow  structure 
decreases  with  increasing  of  cr.  In  particular,  the  distributions  of  both,  the 
particle  concentration  and  the  particle  energy  flux  to  the  surface  along  the 
cylinder  contour,  calculated  at  cr  =  1.2^  are  closer  than  at  cr  =  1.2  to  those 
obtained  for  the  collisionless  particle  phase  model. 

The  same  results  are  obtained  for  the  other  values  of  fp  (for  7‘p  =  25/jm 
and  50/im). 

Thus,  even  rather  large  dispersion  of  a  particle  size  Drp  calculated  for 
the  truncated  distribution  5^(rp),  results  in  the  change  of  macroparameters 
(1)  at  the  front  surface  of  a  cylinder  only  in  10  -4-  20%  compared  with  those 
in  the  flow  of  mono-sized  particles  of  radius  fp.  Therefore,  force  and  erosive 


effects  of  the  particle  phase  on  a  blunt  body  are  expected  to  depend  on  Drp 
slightly  enough. 
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Figure  2.  Distribution  of  the  relative  particle  concentration  cvp/apoo  and  the  relative  particle 
energy  flux  to  the  surface  Ev,/Eoo  along  the  cylinder  contour,  fp  =  100  pm;  curves  1- 
3,  ctpoo  0  (collisionless  particle  phase  motion);  4-6,  <ipoo  =  3  •  10“^  (particle  phase 
motion  with  particle-particle  collisions);  1,  4,  cr  =  1  (monodisperse  particles);  2,  5,  cr  =  1.2; 
3,  6,  cr  =  1.2^.  Drp  is  the  dispersion  of  the  truncated  distribution  function  5^(rp). 
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NUMERICAL  SIMULATION  OF  REFLECTION 
OF  SOLID  PARTICLES  FROM  THE  ROUGH  SURFACE 
IN  TWO-PHASE  FLOW 

S.V.  PANFILOV,  YUM.  TSIRKUNOV 

Baltic  State  Technical  University,  198005  St.  Petersburg,  Russia 

In  two-phase  gas-solid  flow  around  the  bounding  rigid  surface  coarse¬ 
grained  particles  collide  with  the  surface  and  rebound  from  it.  Reflected 
particles  can  essentially  influence  the  two-phase  flow  structure  near  the 
surface.  The  motion  of  reflected  individual  particles  is  largely  determined  by 
their  rebound  parameters  which  depend,  first  of  all,  on  the  particle  collision 
velocity  and  the  particle  collision  angle  given  by  particle  trajectory  inclination, 
as  well  as  on  the  physical  and  mechanical  properties  of  the  particle  and  surface 
materials  [1,2].  At  the  same  time  the  appreciable  dispersion  of  parameters 
of  reflected  particles  is  always  observed  in  experimental  studies  [1].  The 
reasons  of  this  phenomenon  are  random  in  nature.  Among  them  are  particle 
size  distribution  in  real  gas-particle  niixtures,  possible  rotation  of  particles, 
deviation  of  a  particle  shape  from  a  sphere,  as  well  as  surface  roughness. 

Up  to  now  there  are  no  sujSicient  detailed  experimental  data  to  develop 
mathematical  models  of  particle-wall  collisions  taking  account  of  all  factors 
named  above.  In  this  connection  the  analysis  of  every  separate  factor  is  of 
significant  interest,  because  it  allows  to  formulate  and  investigate  assumptions 
which  then  can  be  used  in  more  complete  physical  and  numerical  models, 
taking  into  account  all  important  features  of  the  collision  process. 

At  present  the  two  sufficiently  developed  semiempirical  models  of  a 
spherical  particle  collision  with  the  smooth  wall  are  known  [3,4].  They  are 
based  on  the  momentum  and  the  angular  momentum  equations  for  a  particle 
and  contain  two  experimental  coefficients. 

In  the  paper  [3]  the  model  of  irregular  particle  reflection  from  the  surface 
with  stochastic  roughness  profile  was  also  proposed.  A  local  collision  angle 
was  represented  as  a  sum  of  the  regular  collision  angle  ai  between  the  particle 
collision  velocity  vector  and  the  smooth  surface,  and  the  inclination  7  arising 
due  to  surface  roughness.  The  roughness  profile  was  not  considered  by  itself 
and  7  was  assumed  to  be  a  random  value  with  a  Gaussian  distribution  with  a 
mean  value  of  0°  and  a  given  standard  deviation  of  A7.  Under  a  suitable  value 
of  A7  this  model  gave  very  good  agreement  between  numerical  simulations 
and  experiments  in  two  examples,  namely  on  distributions  of  the  normal 
velocity  restitution  ratio  for  particle  reflection  from  a  target,  and  on  profiles 
of  particles’  velocity  fluctuations  along  the  flow  direction  in  the  plane  channel 
two-phase  turbulent  flow. 
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Nevertheless  the  model  of  irregular  reflection  [3]  calls  for  further  analysis 
and  justification.  First,  the  point  is  that  the  distribution  function  of  random 
angle  7  depends  in  fact  on  ai  for  the  given  roughness  profile,  and  therefore 
it  must  be  related  with  both  the  roughness  profile  and  the  regular  collision 
angle  .  However,  the  question  of  how  to  determine  this  relation  has  been 
only  touched  on  in  [3]  and  remains  open.  Besides,  the  particle  rebound  angle 
distribution  depends  not  only  on  the  distribution  of  local  collision  angles  Ofi  + 
7,  but  also  on  the  kinematics  of  the  reflection  process  (for  example,  on  the 
possibility  of  multiple  collisions  within  one  hollow)  that  can  not  be  taken  into 
account  in  the  model  [3]. 

The  purpose  of  the  present  work  is  the  numerical  simulation  of  the 
particle  reflection  angle  distribution  functions  for  two  types  of  the  roughness 
profile,  periodic  (sine-shaped)  and  stochastic,  and  for  different  values  of  the 
regular  collision  angle  ai . 

The  roughness  profile  and  the  particle  motion  were  supposed  to  be  two- 
dimensional  and  the  particle  shape  to  be  spherical.  The  characteristic  scales  of 
roughness  profiles  (the  amplitude  and  the  period  of  sinusoid  or  the  roughness 
depth  and  the  mean  cycle  of  roughness)  were  given  from  the  ranges  mentioned 
in  [3].  The  stochastic  roughness  profile  y{x)  was  simulated  as  a  cubic  spline 
drawn  through  N  points  with  coordinates  {xi^yi),  that  were  determined  as 
follows:  xi  =  0,  Xi  =  Xi-i  h/2  +  {i  =  2,3,...,^');  =  y,  {i  = 

1,2,...,  N)^  where  the  step  h  was  equal  to  the  mean  cycle  of  roughness,  and 
random  values  ^  and  rj  were  sampled  from  the  Gaussian  distribution  with 
a  mean  value  of  zero.  Standard  deviations  A<f  and  At]  can  be  determined 
from  measurements  of  real  roughness  profiles.  Typical  values  of  the  governing 
parameters  were  taken  as  follows;  iV  =  50,  h  =  100/;/m,  =  bfim  and 

Ay  —  3/im.  In  the  computational  simulation  of  the  distribution  of  particle 
reflection  angles  a2  for  each  considered  value  of  ai  a  number  of  particle- wall 
collisions  was  taken  equal  to  7500  for  the  sinusoidal  roughness  and  about  10000 
for  the  stochastic  one.  In  calculations  the  presence  of  ’’shady”  zones  on  the 
profile  which  particles  can  not  fall  into  at  the  sufficiently  small  value  of  ati  and 
the  possibility  of  multiple  particle-surface  collisions  within  a  roughness  hollow 
were  taken  into  account.  The  semiempirical  particle-wall  colhsion  model  [4] 
was  used  to  determine  the  particle  parameters  after  every  single  reboimd.  The 
particle  radius  was  varied  within  the  range  rp  =  5  -f  50/im  and  the  particle 
collision  velocity  was  taken  equal  to  lOOm/s. 

Typical  patterns  of  the  considered  roughness  profiles  and  the  correspond¬ 
ing  particle  scattering  indicatrixes,  i.e.  the  reflection  angle  distributions  in 
polar  coordinates,  are  shown  (for  rp  =  5/im)  in  the  figure. 

For  the  sinusoidal  profile  there  are  two  dominant  directions  of  particle 
reflection.  The  similar  profile  was  considered  earlier  in  [5].  It  should  be 
noted,  however,  that  the  model  [5]  for  particle  reflections  from  a  sine¬ 
shaped  surface  in  the  case  of  sufficiently  small  ai  is  invalid  because  the 
local  surface  inclination  angle  was  determined  in  [5]  from  the  profile  equation 


Scattering  indicatrixes  of  particles  reflected  from  the  surface  with 
sinusoidal  {a,c,e,g)  and  stochastic  {b,d,f,h)  roughness; 
- ,  regular  reflection  of  particles  from  the  smooth  surface. 
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y  =  a  sm{2nx/ L-\-(p),  where  the  random  phase  (p  was  considered  as  distributed 
uniformly  in  the  range  [0,27r].  Such  distribution  of  (p  is  accurate  only  at 
ai  =90°  but  the  distribution  of  (p  differs  significantly  from  the  uniform  one 
when  ai  is  small. 

It  is  seen  from  the  figure  that  the  configuration  of  particle  scattering 
indicatrixes  for  the  stochastic  roughness  differs  qualitatively  from  the  case  of 
the  periodic  one.  In  spite  of  multiple  outliers,  the  distribution  tends  to  become 
monomodal  which,  possibly,  can  be  obtained  with  significant  increasing  of 
the  surface  length  (with  increasing  of  N).  But  a  mean  value  of  the  particle 
reflection  angle  distribution  will  differ  from  the  regular  angle  0^2,  at  least  for 
small  (Ui  (see  the  figure). 

This  study  was  carried  out  due  to  the  financial  support  of  the  Russian 
Foundation  for  Basic  Research  (grant  JVs  96-01-01467)  and  the  St.  Petersburg 
University  Competitive  Centre  on  Basic  Natural  Sciences  (grant  Ks  95-0- 
4.2-5). 
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MODEL  OF  PROBE  PARTICLES  IN  APPLYING  TO 
GAS-PARTICLE  FLOW  COMPUTATIONS. 

D.  K.  ZAITSEV 

St.Petersburg  State  Technical  University 

Gas-particle  flows  are  found  in  many  industrial  applications  such  as  cyclone 
separators  and  classifiers,  pneimiatic  transport  of  powder,  droplet  and  coal  combustion 
systems,  spray  diying  and  cooling,  sandblasting  and  so  on.  Such  flows  are  usually 
characterised  by  considerable  coupling  between  phases.  For  instance,  when  a  spray  is 
injected  into  a  hot  gas  stream,  one  observes  thermal  coupling  through  heat  transfer  to 
the  droplets,  momentum  coupling  through  aerodynamic  drag  responsible  for  droplet 
motion,  and  mass  coupling  through  evaporation. 

Two  basic  approaches  are  successfully  used  to  model  behaviour  of  discreet  phase 
in  such  flows:  Lagrangian  and  Eulerian.  In  the  Eulerian  approach,  a  "gas"  of  particles 
is  introduced  and  both  phases  are  treated  as  separate  interpenetrating  continua,  each 
one  with  its  specific  density,  velocity,  temperature,  etc.  The  governing  equations  with 
phase  coupling  terms  are  solved  mutually  for  all  phases  (components).  Probably  the 
weakest  point  of  this  approach  is  the  necessity  of  specifying  effective  transport 
coefficients  for  the  "gas"  of  particles  in  turbulent  carrier  flow;  all  the  developed 
closure  models  need  some  experimental  fitting  for  the  particular  flow  type  considered. 
Many  artificial  assumptions  have  to  be  introduced  in  respect  of  the  effect  of  particles 
on  the  carrier  flow  turbulence.  One  more  disadvantage  of  the  Eulerian  treatment  of 
dispersed  phase  concerns  modelling  the  particles-boundaiy  interaction;  in  fact,  the 
boundary  conditions  required  by  the  governing  equations  do  not  correspond  to  the 
physics  of  the  phenomenon. 

The  Lagrangian  approach  is  known  to  be  less  efficient  while  computations  but 
more  physically  correct  and  simple  in  implementation  than  the  Eulerian  one.  In  this 
approach,  a  lot  of  probe  particles  are  tracked  using  a  previously  computed  carrier  flow 
field,  and  an  averaging  procedure  is  employed  to  obtain  local  parameters  of  the 
particle  cloud  and  generate  phase  coupling  terms  needed  for  the  carrier  flow  governing 
equations.  Two-way  coupling  between  the  phases  can  be  achieved  through  global 
iterations.  Effect  of  the  earner  flow  turbulence  is  typically  modelled  by  introducing 
random  fluctuations  of  the  flow  field  while  tracking  the  particles.  Since  all  the  probe 
particles  are  tracked  independently  of  each  other,  the  Lagrangian  approach  is  only 
applicable  to  dilute  mixtures  with  no  or  negli^ble  particle  collisions.  The  conventional 
upper  limit  of  the  volume  fi-action  of  particles  in  gas-particle  flows  is  0,001  that 
corresponds  to  the  mass  loading  of  order  unity.  Below  this  limit  the  model  of  probe 
particles  yields  well  predictions  of  turbulent  gas-particle  flows  while  introducing 
relatively  simple  and  physically  ejear  closure  assumptions. 

Several  modifications  of  the  model  of  probe  particles  are  reported  in  the 
literature;  the  formulations  4j®5r  from  each  other  by  some  details  of  generating  random 
quasi-turbulent  flow  field  as  well  as  by  the  physical  effects  taken  into  account  while 
particle  tracking.  A  simple  model  presented  below  is  aimed  at  simulation  of  a  steady- 


state  dilute  gas-particle  flow  with  solid  spherical  particles;  heat  and  mass  transfer 
between  the  phases  are  not  considered. 

Motion  of  a  probe  particle  is  governed  by  the  equation  (1).  The  aerodynamic 
force  on  the  particle  is  represented  by  the  drag  one  (2);  other  forces  are  negligible 
because  of  hi^  particle-to-gas  density  ratio. 

dv  _  =  df  ^ 

F  =  37in8(u-v)(l  +  0.179^/Re;;  + 0.013 Rep),  Rep<6-10^  (2) 

To  simulate  particle  dispersion  due  to  turbulence,  the  acmal  velocity  u  of  the 
carrier  gas  along  the  particle  trajectory  is  represented  as  the  sum  of  the  local  mean 
velocity  (available  from  the  computed  flow  field)  and  some  fluctuation  (3).  The 
fluctuation  is  randomly  sampled  from  the  normal  (Gaussian)  distribution  with  mean 
zero  and  standard  deviation  according  to  the  turbulent  energy  k;  the  turbulence  is 
assumed  isotropic.  The  sampled  fluctuation  is  associated  with  an  eddy  whose  size  U 
and  life-time  tg  are  evaluated  through  the  turbulent  energy  k  and  the  rate  of  its 
dissipation  e;  the  constant  C-r=0.28  was  chosen  by  comparing  numerical  results  with 
fundamental  derivation  by  Hinze  (1975)  for  turbulent  dispersion  of  light  particles. 


u  =  u  +  u' 


=  ^k; 
3 


It  is  accepted  that,  the  generated  fluctuation  remains  unchangeable  within  the 
eddy  during  its  life-time.  As  soon  as  the  life-time  has  expired  or  the  particle  has  left  the 
eddy,  a  new  fluctuation  is  chosen  assuming  no  correlation  with  the  previous  sample. 

Neglecting  fluctuations  of  gas  density  and  the  volume  fraction  of  particles  in  the 
mixture,  the  carrier  flow  field  is  governed  by  the  equations  (4).  They  differ  from 
corresponding  single-phase  equations  by  particle-induced  source  terms  (5);  the  comer 
brackets  denote  time  averaging  along  a  particle  trajectory  and  summation  over  all  the 
probe  particles  within  a  control  volume.  Except  the  last  one,  the  source  terms  are 
available  from  particle  tracking  and  do  not  need  artificial  modelling. 


V .  (pu)  =  0  ;  V .  (pu  ®  Q  +  pT  -  =  Fp 


V .  puk  -  -^Vk  =P-p8  +  P 


V .  ^pu8  -  ^  Vej  =  -(C,,P  -  C,2PS  +  C,3Pp) 


Fp=-F  ;  P,=-(u'.F);  C3-I.0 


(5) 
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The  model  has  been  implemented  using 
arbitrary  structured  multiblock  3D  grids. 
Several  tests  were  made  to  validate  the 
developed  code.  Particularly  in  the  diagram 
to  the  left,  numerical  results  on  dispersion  of 
separate  particles  in  a  grid-generated 
turbulence  are  compared  with  measurements 
by  Snyder  &  Luinley  (1971).  The  carrier 
flow  field  was  specified  according  to  the 
experimental  data.  Injection  of  10000  probe 
0  100  200  300  400  particles  provided  statistically  stationary 

Time  from  station  one  (msec)  solution  with  the  true  Gaussian  distribution 

of  particle  deflection  from  the  flow  axis. 

As  another  test,  computations  were  made  for  a  planar  mixing  layer  with 
polydispersed  (3-lOOmkm)  water  drops  according  to  flow  conditions  reported  by 
Chang  et  al.  (1993).  Unlike  the  former  test,  both  the  carrier  flow  and  particle 
trajectories  were  simulated;  measured  flow  parameters  were  used  as  inlet  boundary 
conditions  while  at  computations.  Following  the  reported  data,  8  fractions  of  probe 
particles  were  introduced.  To  obtain  a  statistically  stationary  solution,  19200 
ti-ajectories  were  computed  (2400  particles  in  each  of  the  eight  fi-actions;  800  particles 
per  cell  at  the  inlet).  The  prediction  was  encouragingly  accurate  except  the  lateral 
component  of  the  carrier  gas  velocity;  as  a  result  the  lateral  velocity  of  particles  was 
slightly  underestimated  too. 

The  developing  region  (x,0<12)  of  a  particle  laden  round  jet  was  investigated  by 
Mostafa  et  al.  (1989);  the  measurements  provided  a  basis  for  testing  the  model  at 
relatively  high  mass  loading  of  the  flow  (up  to  unity).  13200  probe  particles  were  used 
in  computations.  A  satisfactory  agreement  between  measured  and  computed  flow 
fields  was  obtained  despite  rather  unreliable  data  on  the  inlet  conditions. 

The  work  has  been  earned  out  under  support  of  “Research  in  Brussels  95”  grant, 
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LASER  DOPPLER  SYSTEM  FOR  DIAGNOSTICS  OF  HIGH¬ 
SPEED  KINETIC  FLOWS 

ANISIMOV  YU.l,  AGAPOV A.S,  LASHKOVVA.,  MASHEKI.CK, 
MASHEKA.CH. 

St.-Petersburg  State  university 

High-speed  kinetic  flows  of  dispersed  particles  (flows  of  particles  by  the 
sizes  5-300  microns  witli  speeds  2-15  Km/s  without  gas),  are  investigated  now 
rather  heavily. 

Such  flows,  created  in  laboratory  conditions,  can  simulate  micrometeoric 
influence  on  elements  of  space  veliicles,  to  serve  for  researches  in  the  field  of  hi^- 
speed  impact,  to  be  used  in  technological  purposes.  Diagnostics  of  particle 
parameters  over  body,  reflected  and  scattering  flows  of  matter  at  siuface  of  target  is 
tlie  special  interest  at  research  of  processes  interaction  of  such  kinetic  flows  with 
target. 

The  laser  system,  serving  for  diagnostics  of  the  specified  flows,  consists  of 
tlie  single  fi'equency  tuning  and  stabilized  Argon-ion  CW  laser  with  output  power  up 
to  100  mW,  the  system  of  illumination  of  researched  area  of  a  flow,  the  gathering 
system  of  scattering  radiation,  the  absorptive  frequent  detector  on  vapor  molecular 
iodine- 127,  two  channels  of  measurement  of  scattering  light:  measuring  (past 
tlirough  frequent  detector)  and  basic  (falling  on  detector),  two-channel  fast  ADC 
and  the  computer.  The  main  principles  of  operation  of  Doppler  systems  with  tlie 
fi'equent  molecular  detector  and  their  application  to  diagnostics  of  stationary  flows 
are  described  in  the  papers  /I,  2/.  Peculiarity  of  the  system  described  in  the  present 
work  is  its  high  temporary  and  space  resolution,  received  at  the  expense  of 
application  of  fast  light  detectors  with  threshold  frequency  more  than  100  MHz,  and 
fast  ADC  witli  operation  time  25  ns.  The  specified  parameters  of  channels  of 
measurement  of  scattering  light  have  allowed  to  receive  measuring  volume  with 
characteristic  sizes  0,5x0, 5x1,0  mm  at  diagnostics  of  flows  with  speeds  up  to  15 
Km/s.  A  particle,  passing  through  area  of  measurement  scatters  laser  light,  forming 
pulse  of  radiation  by  duration  lO'”^  -  10'^  s.  The  system  measures  scattered  light  in 
two  channels  simultaneously  -  basic  and  measuring.  Their  attitude  does  not  depend 
on  ^plitude  of  light  pulse  and  is  determined  by  Doppler  shift  frequency  of 
scattering  radiation  and  kind  of  the  gear  characteristic  of  the  absorption  frequency 
detector.  The  subsequent  numerical  processing  of  signals,  made  in  view  of 
temporaiy/  width  of  each  pulse  spectrum,  permits  to  receive  distribution  of  particle 
speeds.  The  quantity  of  registered  particles  for  one  cycle  of  measurement  is 
determined  by  volume  of  mtemal  memory  used  ADC  (four  F4226  in  standard 
CAMAC,  two  ADC  on  channel,  started  synchronously,  with  sliift  25  ns).  Maximum 
time  of  measuring  process  for  one  cycle  is  50  -  100  mks.  The  described  system  at 
first  stage  permits  to  measure  only  one  projection  of  speed  of  researched  flows, 
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however  potential  opportunities  of  a  described  method  permit  to  create  the  similar 
analyzer  as  for  two,  and  for  three  components  of  a  vector  of  speed. 
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COMPARATIVE  NUMERICAL  ANALYSIS  OF  VISCOUS  AND 
INVISCID  FLOWS  IN  A  SQUARE  CAVITY 


Lobov  S.L.,  Romanov  0.  Y. 

The  Baltic  State  Technical  University 

The  main  features  of  the  inner  flows  in  SRM  chamber  are  the 
recirculation  zones.  These  zones  arise  due  to  main  body  of  the  gas  flow 
separating  from  a  rigid  wall,  which  caused  by  joint  action  of  positive  pressure 
gradient  and  viscosity. 

Since  the  gaseous  products  of  solid  propellant  burning  have  a  very  small 
viscosity  (large  Re),  the  numerical  computation  of  the  Navier-Stokes  equations 
present  great  difficultes.  The  first  is  that  the  conventional  schemes  with  small 
numerical  visconsity  have  very  strong  stability  limitation  on  the  grid  size  and 
time  step  [  1  ],The  second,  in  each  gridpoint  finite-difference  solution  is 
qualitatively  viscous  only  when  local  grid  Reynolds  number  not  more  than  2  and, 
consequently,  the  local  grid  size  is  comparatable  with  kinematic  viscosity 
coefficient  [  2  ],  That  why,  a  method  that  enables  to  obtain  a  stable  and  accurate 
solution  is  limited  by  the  computing  time  and  storage  available. 

One  of  the  possible  ways  to  solve  this  problem  is  to  use  an  inviscid 
vortical  flow  model.  The  above  assumption  is  undoubtedly  true  in  cases  where 
the  streamlines  are  open-ended.  There  are  some  cases,  in  which  the  main  body 
separating  from  the  ridig  wall  is  determinated  by  boundary  conditions  without 
apparent  viscosity  influence.  For  instance,  this  takes  place  in  a  case  of  flow  in  a 
region  bounded  extermaly  by  ridig  wall  with  rectangular  elementes,  or  in  a  point 
of  joint  of  mass  injection  part  of  boundary  and  solid  wall.  Batchelor  G.  K.  has 
theoretically  proved  that  for  small  but  non-zero  values  of  kinematic  viscosity  in 
two-dimensional  closed  flow  the  central  part  asymptotically  approached  to 
'inviscid  core'  with  uniform  vorticity.  However,  the  r^ation  between  value  of 
vorticity  inside  the  recirculation  zone  (WO)  and  the  same  on  boundary  line  of 
main  stream  body  (Wl)  remains  to  be  determinated  by  the  conditions  that  the 
viscous  boundaiy  layer  surrouding  this  region  must  also  be  in  steady  motion  [3]. 

The  parametric  numerical  investigations  based  on  viscous  and  on  inviscid 
models  have  been  fulfiled  to  definite  above  mentioned  relation  and  to  find  a 
lowest  Reynolds  value,  since  which  the  vorticity  in  the  centre  of  cavity  became 
quasi  uniform.  The  classical  cavity  flow  problem  has  been  modified  by  the 
entering  of  boundary  part  of  main  flow  body  to  the  computational  domain.  The 
scale  size  equal  to  the  bight  of  cavity  (H)  .  The  computional  grid  is  rectangular 
and  imiform.  The  equation  system  is  written  in  term  of  a  vorticy-stream  fimction. 
The  task  is  resolved  in  axisimmetric  variant  but  at  large  values  of  radial 
coordinates  (R>8H).  The  govering  equations  are  used  here  in  implicit  finite- 
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difference  formulation.  All  metric  derivatives  are  evaluated  by  central  difference 
formulas. 

The  calculations  demonstrated,  that  at  the  Re=800  in  vorticity  distribution 
in  central  part  of  cavity  the  plate  appeared.  When  the  Reynolds  number  further 
grows,  this  plate  increases  and  occupies  the  most  of  the  cavity  space.  In  the  main 
flow  zone  vorticity  distribution  apparently  trends  to  consideration  of  vorticy 
along  streamlines  since  Re>200.  The  relation  (WOAVl)  is  near  5-6. 

The  accurate  viscous  flow  solutions  obtained  for  the  peculiar  geometric 
elements  then  applied  to  solve  the  global  flow  problem  in  whole  domain  of  SRM 
inner  volume  using  approximate  inviscid  model.  Thus,  the  presented  technique 
ensures  the  required  accuracy  at  resonable  computing  time  and  storage 
requirements. 
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COMPUTATION  OF  TWO-PHASE  FLOW  AND  SLAG 
MOTION  IN  SOLID  ROCKET  MOTOR 

Lobov S.L.,  Romanov  O.Y. 

The  Baltic  State  Technical  University 

The  combustion  of  an  alumized  solid  propellant  produces  not  only  gases 
but  also  liquid  particles.  These  particles  are  the  source  of  slag  material  that  may 
remain  in  the  motor  during  firing  and  thereby  cause  excessive  heating  and 
subsequent  insulation  erosion  where  they  collect. 

A  time-dependent,  axisymmetric  numerical  solution  technique  has  been 
developed  to  analyse  the  coupled  gas-particle  flow  and  the  motion  liquid  slag  in  a 
solid  rocket  motor. 

In  this  work  for  description  of  two-phase  flow  numerical  model  based  on 
a  combined  Eulerian-Lagrangian  analysis  was  used  [1].  The  time-dependent  flow 
description  was  ensures  by  simulating  the  quasi  steady-state  flowfield  at 
sequence  of  bum  times. 

The  previously-computed  single-phase  flow  results  are  taken  as  the  initial 
quess  for  the  two-phase  flow.  This  quasi  gas  phase  solution  is  solved  in  Eulerian 
mesh.  Initialy  the  gas  is  inviscid,  vortical  and  incompressible,  without  chemical 
reactions.  The  mathematical  formulation  for  gas  is  used  here  in  term  of  vorticity- 
stream  function. 

This  model  treats  the  particles  as  a  discrete  phase.  The  discrete  particle 
phase  solution  is  obtained  by  integrating  their  dynamical  equation  of  motion 
along  their  trajectories.  The  combustion  of  an  aluminized  solid  propellant 
produces  AL/AL2O3  agglomerates  and  small  particles  on  the  propellant  burning 
surface,  which  are  ejected  into  the  gaseous  surrouding.  Along  agglomerate 
trajectoiy'  remained  in  agglomerate  metal  is  oxidized.  Some  of  the  new  AL2O3  is 
shed  as  small  particles,  while  the  rest  remains  attached  to  the  AL/AL2O3 
agglomerate.  This  events  lead  to  bimodal  size  distribution  of  k-phase,  consisting 
of  large  AL/AL2O3  composite  particles  and  small  'smoke'  particles  [2],  The 
'smoke*  particles  have  neglectible  values  of  velocity  lag  from  gas,  that  why  the 
particle  phase  analysis  is  only  by  agglomerate  solution  restricted.  In  this  study 
agglomerate  burning  and  breacup  are  accounted  for.  In  a  SRM,  the  particle-to¬ 
gas  density  ratio  is  usually  lage,  ^erefore,  the  only  important  forces  on  a  particle 
are  the  inertia,  drag  and  external  acceleration. 

In  model  the  relationship  between  the  metal  burning  rate  and  spread 
gaseous  oxidizer  density  distribution  in  inner  volume  of  SRM  is  also  account  for. 

The  time-dependent  slag  motion  solution  is  solved  in  a  Eulerian  mesh. 
The  slag  material  is  considered  as  chemically  inertial  film  flow  along 
axisimmetric  surface.  The  important  forces  on  a  film  are  the  gas  pressure 
gradientes,  the  viscous  friction  between  gas  and  film  and  external  acceleration. 
The  local  momentum  and  mass  flux  density  also  influences  on  a  slag  motion. 
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The  particle  flow  and  slag  motion  information  has  been  used  to  estimate 
the  amount  of  slag  in  a  apogee  SRM,  which  agrees  favourably  with  availible  test 
data.  The  results  of  numerical  investigation  were: 

1) the  distribution  of  mass  flux  density  for  oxide  and  nonbumed  metal 
along  the  countour  of  aft  chamber  bottom  and  submerged  nozzle; 

2)  space  distribution  of  combustion  effeciency  for  solid  propellant  charge 

body; 

3)  the  time-dependent  slag  form,  determinated  by  slag  material  motion. 
The  prezented  numerical  solution  technique  gives  the  possibility  to 

determinate  the  behavior  of  the  particle  laden  gas  as  it  moves  in  the  motor,  and 
the  slag,  which  remains  in  the  chamber.  This  knowledge  is  neccesary  for 
optimum  design  of  SRM. 
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COMPUTATION  OF  CHARACTERISTICS  OF  LARGE-SIZE 
OXIDE  PARTICLES  IN  CHAMBER  OF  ALUMINIZED  SOLID 
PROPELLANT  MOTOR 
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Abstract 

Combustion  of  aluminized  solid  propellants  (ASP)  is 
characterized  by  formation  of  particles  of  two  different  sizes.  Large 
size  particles  may  reach  up  to  hundreds  and  thousands  micrometers. 
It  is  these  particles  that  cause  slag  deposition  in  rocket  chambers  and 
some  accompanying  phenomena.  Large-size  oxide  particles  are 
formed  as  the  result  of  evolution  of  agglomerates  in  the  combustion 
product  flow  after  they  leave  the  burning  propellant  surface.  A  model 
for  large-size  particle  formation  is  developed  as  a  result  of 
experimental  and  theoretical  investigation.  Parametric  analysis  of  the 
model  as  applied  to  a  large  dimension  rocket  engine  had  revealed  that 
characteristics  of  large-size  oxide  particles  depend  on  propellant 
properties,  engine  geometry,  time  of  their  presence  in  a  combustion 
chamber.  The  model  allows  to  predict  large-size  particles  properties, 
to  improve  the  technique  for  calculation  of  two  phase  flows  in  the 
chamber  and,  in  turn,  to  obtain  more  precise  characteristics  of  a 
"propellant-engine"  system. 


Introduction 

Burning  of  aliminized  solid  propellant  (ASP)  produces  condensed 
phase  combustion  products  (CPCP)  which  consist  of  two  types  of 
particles:  large-size  and  small-size.  The  former  may  reach  hundreds 
and  even  thousands  micrometers,  the  latter  is  about  1  micrometer  in 
size.  In  the  rocket  chamber  the  small-size  particles  are  in  equilibrium 
with  gas  phase.  The  large-size  particles  are  characterized  by  their 
velocity  lag  as  compared  with  gas.  It  is  precisely  these  particles  that 
cause  slag  deposition  on  the  inner  surfaces  of  the  chamber  and 
accompanying  phenomena  connected  with  additional  transfer  of 
substance  and  energy,  i.e.,  more  intensive  destruction  of  heat  isolation 
and  slag  deposition. 
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The  CPCP  are  formed  in  the  course  of  two  processes:  metal 
agglomeration  on  the  burning  propellant  surface  and  CPCP  evolution 
in  the  combustion  product  flow.  Agglomeration  leads  to  the  fact  that 
the  CPCP  leaving  the  burning  propellant  surface  consist  of 
agglomerates  and  highly  dispersed  oxide  particles  (HDOP)  of  about  1 
micrometer  in  size.  Quantitative  proportion  of  agglomerates  and 
HDOP  and  their  dispersity  are  determined  by  peculiarities  of  the 
agglomeration  process.  The  evolution  of  CPCP  includes  a  number  of 
physical  and  chemical  transformations  of  the  particles  in  combustion 
product  flow.  This  paper  is  devoted  to  modeling  of  large  size  particle 
formation  in  combustion  chamber  with  agglomeration  characteristics 
assigned.  To  solve  this  problem  it  is  necessary  to  model  the  evolution 
of  a  single  agglomerate  and  the  evolution  of  two  phase  flow. 

Modeling  of  a  Single  Agglomerate  Evolution 

The  model  was  developed  on  the  base  of  various  experimental 
data  on  agglomerate  characteristics  at  different  stages  of  evolution  [1]. 
The  model  takes  into  account  the  following  phenomena: 

-  burning  of  agglomerate  metal  in  gaseous  phase  mode, 

-  chemical  interaction  of  metal  and  oxide  in  agglomerate  with 
formation  of  gaseous  products; 

-  structure  of  agglomerate  composed  of  metal,  oxide,  and  gaseous 
products  in  different  proportion; 

-  deposition  on  agglomerates  of  HDOP  from  the  flow. 

The  developed  model  combines  a  number  of  individual  models 
for  specific  phenomena  occurring  in  the  evolution  [1,2].  As  these 
models  are  rather  complex  and  require  much  computer  time  when 
used  then  fast-operating  formal  models  were  developed  for  some 
specific  phenomena.  Reliability  of  single  agglomerate  evolution 
model  was  confirmed  by  comparison  of  the  model  and  experimental 
data  on  evolution  of  agglomerates  in  linear  flow. 

Modeling  of  Two  Phase  Flow  Evolution 

Two  phase  flow  consists  of  two  interacting  objects:  1.  gas  phase 
with  HDOP  and  2.  large-size  particles  of  CPCP.  Chemical,  termal, 
and  dynamic  interactions  between  this  objects  take  place  in  the  flow 
that  causes  changes  in  velocity,  chemical  contents,  and  temperatures 
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of  the  objects.  In  the  particle  burning  zone  HDOP  are  formed  and 
then  transported  to  gas  phase.  In  turn,  HDOP  of  the  gas  phase  may 
deposit  on  the  large-  size  particles. 

Modeling  of  the  flow  evolution  in  rocket  chamber  is  based  on  the 
following  approach.  Velocity  field  of  the  gas  phase  is  calculated  on 
the  base  of  a  model  for  vortex  motion  of  non  viscous  and 
incompressible  gas  developed  by  V.N.Emelyanov  at  Baltic  State 
Technical  University  (Saint  Petersburg,  Russia).  Then  on  the  base  of 
the  gas  velocity  field,  motion  of  large-size  particles  is  calculated.  In 
this  procedure  the  model  of  single  agglomerate  evolution  is  used  that 
allows  to  take  into  account  all  the  phenomena  occurring  in  the 
evolution.  Calculations  are  made  for  different  fractions  of  the 
particles  and  for  different  places  on  propellant  burning  surface  which 
produce  the  particles. 

Chamber  cavity  is  divided  into  sells.  Dimensions  of  a  sell  depend 
on  precision  of  calculations.  Temperature,  oxidizer  and  HDOP 
concentrations  in  each  sell  are  calculated  as  a  result  of  evolution  of 
agglomerates  crossing  a  gas  stream  which  passes  through  the  sell.  This 
parameters  are  determined  with  the  help  of  equilibrium 
thermodynamic  approach. 

A  problem  of  modeling  HDOP  formation  in  solid  propellant 
combustion  has  not  been  solved  completely  yet.  However,  an 
experimental  data  are  obtained  which  allow  to  think  that  HDOP 
dispersity  depends  on  sizes  of  the  agglomerates.  At  the  burning 
propellant  surface,  mass-medium  diameter  843  of  HDOP  is  about  0.8  - 
0.9  pm.  These  HDOP  are  formed  as  a  result  of  burning  process  of 
metal  particles  which  do  not  take  part  in  the  agglomeration  [3]. 
Agglomerate  combustion  may  change  this  value  up  to  4  -  5  pm  [4,5].  It 
may  be  assumed  also  that  agglomerate  lag  influences  dispersity  of 
HDOP  formed  in  agglomerate  burning  process.  These  facts  allow  to 
believe  that  HDOP  dispersity  may  change  significantly  during  two 
phase  flow  evolution  in  rocket  chamber.  In  this  paper,  HDOP 
dispersity  was  accepted  to  be  constant  during  evolution  but  dependent 

on  the  portion  of  metal  involved  in  agglomeration  (  Z^);  HDOP 

diameters  increase  if  parameter  increases. 

In  computer  modeling  it  was  accepted  that  agglomerate 
deposition  takes  place  if  agglomerate  trajectory  crosses  inner  surface 
of  combustion  chamber.  However,  behavior  of  deposited 
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agglomerates  is  not  examined.  The  model  allows  to  calculate 
parameters  of  large-  size  fraction  of  CPCP  in  any  section  of  the 
chamber.  Input  information  for  the  calculation  is  the  following: 

-  geometry  of  propellant  charge  and  combustion  chamber; 

-  propellant  burning  rate  and  chamber  pressure; 

-  parameters  of  CPCP  at  propellant  burning  surface  (proportion  of 
agglomerates  and  HDOP,  agglomerate  size  distribution); 

-HDOP  size; 

-  thermodynamic  characteristics  of  the  propellant. 

Computer  Modeling 

Computer  program  was  developed  for  above  model.  The 
program  was  used  for  calculating  behavior  of  large  -  size  CPCP 
particles  in  the  last  stage  of  their  evolution  in  rocket  chamber.  CPCP 
evolution  in  combustion  chamber  is  considered  to  end  with  either 
agglomerate  deposition  on  chamber  inner  surface  or  agglomerate 
entry  into  motor  nozzle  where  their  evolution  continues  in  supersonic 
two-  phase  flow. 

Results  of  the  agglomerate  deposition  are  represented  in  Fig.l  and 
2.  Fig.  1  shows  section  of  rocket  chamber  and  propellant  burning 
surface.  Computations  were  made  for  the  following  conditions: 
pressure  P  =  6MPa,  acceleration  1.5g,  oxide  part  in  agglomerate  by 

mass  ri=0.35,  parameter  =0.72,  HDOP  diameter  5=2|um, 

agglomerate  mass-medium  diameter  D43  =327 pm.  Fig.  2  shows 
distribution  of  CPCP  deposition  along  the  inner  surface  ABCDEFGH 
of  the  chamber.  The  deposited  CPCP  consist  of  80%  oxide  by  mass. 

The  study  allows  to  find  out  characteristics  of  propellant  and 
combustion  chamber  which  strongly  influence  deposition  intensity. 
Obvious  ways  to  deposition  reduction  are  decrease  of  gas  velocity, 
curvature  of  flow  lines,  agglomerate  size  and  their  quantity.  The  same 
task  may  be  solved  with  unusual  approaches:  by  temperature  increase 
and/or  organization  of  local  zones  with  oxidizing  potential  increased 
in  gas  phase.  These  factors  may  cause  intensification  of  reaction 
between  aluminum  and  alumina  in  the  agglomerate,  and,  in  a  limiting 
case,  fragmentation  of  agglomerates. 

Propellant  features  influence  parameters  of  large-  size  CPCP 
changing  characteristics  of  agglomeration  process.  In  turn, 
parameters  of  two-  phase  flow  at  propellant  burning  surface  have 
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strong  influence  on  the  flow  evolution  process.  Decrease  in  size  of 
agglomerates  and  oxide  content  in  them  leads  to  intensification  of 
reaction  between  aluminum  and  oxide  in  agglomerates.  Reductions  in 

Zm  parameter,  HDOP  size,  and  agglomerate  velocity  lags  cause 
decrease  of  HDOP  deposition  on  the  agglomerates. 


Figure  1.  Section  of  solid  propellant  motor  chamber. 
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CPCP  evolution  in  one-dimensional  flow  was  calculated  for 
different  two  phase  flow  conditions  provided  by  using  two  different 
propellants:  Prop.N  1  and  Prop.N2,  Figures  3  and  4  show  results  of 
the  calculation:  agglomerate  size  distribution  at  propellant  burning 
surface  and  in  the  last  stage  of  the  evolution  process  when  combustion 
efficiency  reaches  up  to  95% 


Figure  3.  Comparison  of  particle  size  distribution  for  Prop.Nl. 

Prop.N  1  is  characterized  by  the  following  initial  parameters  of 

CPCP  (near  propellant  surface):  D43=l30nm,  Z^=:0.26,  ti=0.34. 

These  values  provide  predominance  of  oxide  removal  from 
agglomerates  in  comparison  with  HDOP  deposition.  Therefore  in  the 
case  of  Prop.Nl,  CPCP  evolution  is  accompanied  by  decreasing 
agglomerate  sizes  in  consequence  of  metal  burning  and  oxide  removal. 
As  a  result,  dispersity  of  large-  size  particles  increased  and  the 
distribution  function  changed  as  it  is  shown  in  Fig.3.  Initial 
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parameters  of  CPCP  for  Prop.N2:  D43=327|um,  =0.72,  r)=0.35 

provide  much  lower  rate  of  oxide  removal  and  more  intensive 
deposition  of  HDOP  on  agglomerates.  These  facts  reduce  dispersity 
of  large-  size  particles  (Fig.  4). 


Figure  4.  Comparison  of  particle  size  distribution  for  Prop.N2. 


Structure  of  large-  size  CPCP  particles  is  changing  during 
evolution  (Fig.5).  When  leaving  the  burning  propellant  surface, 
agglomerates  are  aluminum  drops  with  oxide  ’’caps"  [6].  Agglomerate 
may  contain  gas  cavities  (bubbles)  which  are  results  of  gaseous 
products  formation  in  the  chemical  reaction  between  metal  and  oxide 
of  agglomerates.  The  agglomerate  structure  is  near  to  equilibrium  and 
it  is  determined  by  surface  properties  and  amounts  of  the  substances 
in  the  agglomerate.  Intensification  of  the  reaction  and  decrease  in 
"cap"  thickness  may  lead  to  out-break  of  bubbles  through  the  "cap" 
and,  consequently,  to  oxide  removal  from  the  agglomerate.  At  the 
last  stage  of  the  evolution,  the  above  reaction  may  continue  but  oxide 
removal  stops  in  consequence  of  agglomerate  structure 
transformation  (oxide  "cap"  thickness  increasing)  because  of  metal 
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burning  out.  Thus,  conditions  for  bubble  existence  in  agglomerates 
remain  in  final  stage  of  the  evolution.  The  bubbles  merge  and  after 
metal  burned  out,  agglomerates  turn  into  oxide  drop  with  one  bubble, 
as  a  rule.  Size  of  the  final  bubble  depends  on  pressure,  particle  size 
and  temperature.  The  bubble  may  amount  to  25  -  40  %  of  the  oxide 
drop  volume  and  be  not  located  symmetrically  in  the  drop. 


motion  direction 


Figure  5.  Scheme  of  agglomerate  transformation  to  final  oxide 

particles. 


Thus,  at  the  last  stage  of  the  evolution,  large-size  particles  of 
CPCP  are  oxide  drops  containing  bubbles  displaced  from  the  drop 
center  in  direction  of  the  drop  motion.  This  drop  structure  is  rather 
steady  for  a  wide  range  of  conditions  of  two  phase  flow  of  combustion 
products  in  rocket  chamber. 

Conclusion 

The  work  is  based  on  contemporary  knowledge  about  large-  size 
particle  formation  in  rocket  solid  propellant  combustion.  The  work 
resulted  in  developing  methods  for  prediction  of  large-  size  particles 
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characteristics  in  rocket  chamber.  The  result  of  the  work  would  be 
useful  for  optimization  of  the  "propellant  -  motor"  system. 
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ON  REASONS  FOR  STABILITY  LOSS  OF  MAGNESIUM  VAPOR- 
PHASE  COMBUSTION 

U.I.  GOLDSHLEGER,  E.Ya.  SHAFIROVICH 
Institute  of  Stmctural  Macrokinetics,  Russian  Academy  of  Sciences 
Chemogoiovka,  Moscow  142432,  Russia 

Combustion  of  metals,  in  spite  of  the  often  observed  imsteady  phenomena  (a 
decrease  in  light  emission  intensity,  jets,  explosions,  fragmentation),  is  usually 
considered  a  quasi-steady  vapor-phase  or  heterogeneous  process.  It  is  generally 
agreed  that  for  volatile  metals  ignition  always  leads  to  stable  vapor-phase  burning. 
However,  in  experimental  investigation  of  magnesium  interaction  with  CO,  a  pulsating 
combustion  mode  has  been  unexpectedly  found  [1].  Later,  in  thermogravimetric 
studies  of  Mg  oxidation  in  carbon  oxides,  a  "slow  heterogeneous  combustion"  mode 
has  been  found,  with  kinetic  self-acceleration  of  heterogeneous  reaction  being  inherent 
[2,  3].  In  CO2  atmosphere  the  temperature  existence  region  of  such  a  mode  separates 
slow  oxidation  and  fast  vapor-phase  burning,  the  low  boundary-  of  tltis  region 
coinciding  with  the  melting  point  of  Mg.  In  CO  atmosphere  a  high-temperature  mode 
of  "slow  heterogeneous  combustion",  apart  from  the  low-temperature  one,  is  observed 
that  separates  regions  of  pulsating  and  vapor-phase  combustion. 

In  the  present  paper  a  possible  mechanism  of  interaction  between  macro- 
kinetically  characteristics  and  physical/chemical  processes  in  a  surface  film  and  gas 
phase  is  discussed  that  makes  it  possible  to  explain,  in  the  framework  of  a  unified 
concept,  the  existence  of  the  unusual  combustion  modes  of  Mg  mentioned  above. 

Analysis  of  experimental  results  has  shown  that  the  slow  heterogeneous 
combustion  process  is  mainly  detennined  by  heat  and  mass  transfer  of  liquid 
magnesium  through  a  growing  porous  oxide  film  and  by  its  evaporation  from  smface 
of  menisci  in  pores.  Good  wetting  by  liquid  Mg  of  its  own  oxide  results  in  continuous 
pulling  of  metal  into  the  porous  film.  The  pressure  difference  needed  for  tlie  liquid 
movement  is  caused  by  capillary  forces  on  the  liquid/gas  interface  in  the  pores  and  on 
surface  of  a  bubble  that  forms  inside  the  metal  drop  and  grows  during  the  process. 
Due  to  high  velocities  of  capillary  flow,  the  vapor-gas  reaction  occurs  close  to  the 
film/oxidizer  interface.  The  growth  of  the  film  leads  to  increasing  the  effective  surface 
of  evaporation  and,  therefore,  to  macro-kinetic  acceleration  of  chemical  reaction. 
Moreover,  with  increasing  the  bubble  the  curvature  of  concave  liquid  menisci  in  the 
pores  decreases.  Pressure  of  vapor  above  the  metal  surface  increases  (Thomson 
effect).  This  is  an  additional  reason  for  the  kinetic  self-acceleration  of  reaction  during 
slow  heterogeneous  combustion. 

Let  us  consider  features  of  metal  vapor-phase  combustion  in  strong  (CO2)  and 
weak  (CO)  oxidizers. 

During  combustion  of  Mg  in  CO2  the  gas-phase  reaction  occurs  in  a  thin  zone 
detached  from  the  particle  surface.  The  great  heat  flux  from  high-temperature  flame 
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results  in  the  high  rate  of  metal  vaporization.  The  vaporization  occurs  in  Stefan  mode 
with  noticeable  gas  velocities  in  space  between  the  particle  surface  and  the  flame 
zone.  For  this  reason,  the  oxidizer  concentration  close  to  the  particle  surface  is  low 
and  the  oxide  fihn  staicture,  fonned  at  initial  phases  of  the  process,  does  not  change 
during  vapor-phase  combustion.  The  burning  phase,  in  spite  of  possible  formation  of 
an  "outer"  oxide  shell  in  the  flame,  may  be  considered  a  stable  quasi-steady  process. 

During  combustion  of  Mg  in  CO,  when  the  adiabatic  flame  temperature  is  low, 
tlie  gas-phase  reaction  occurs,  on  the  contrar>^,  occurs  in  an  extended  zone  adjusting 
the  particle  surface.  In  this  case  the  temperature  difference  between  the  flame  zone 
and  the  particle  surface  is  noticeably  less.  Since  the  heat  flux  from  the  flame  is  small, 
the  evaporation  mode  is  close  to  the  diffusion-controlled  one.  Therefore,  the  counter¬ 
flow  of  Mg  vapor  (that  prevents  ingress  of  CO  to  the  particle  surface)  is  significantly 
weaker  than  in  the  combustion  in  CO2.  For  this  reason,  in  the  mode  of  Mg  vapor- 
phase  combustion  in  CO,  together  with  the  gas-phase  reaction  in  volume,  the 
heterogeneous  reaction  may  always  occur  on  the  paiticle  surface  at  the  temperature 
close  to  the  boiling  point  of  Mg.  The  increase  in  the  particle  surface  during  vapor- 
phase  inflammation  promotes  great  thennal  acceleration  of  the  heterogeneous  reaction. 
On  the  other  hand,  the  characteristic  grain  size  of  the  suiface  film  falls  exponentially 
with  increasing  the  temperature,  resulting  in  drastic  changes  in  the  film  properties. 
This  may  entail  tiie  great  kinetic  inhibition  of  the  heterogeneous  reaction  that  is  not 
compensated  by  the  higher  rate  of  heat  release. 

The  changes  in  the  film  properties  are  an  inhibiting  cause  with  respect  to  both 
heterogeneous  and  gas-phase  reactions.  With  the  capillary  mechanism  of  metal 
transport  to  the  outer  boundary  of  the  film,  the  decrease  in  the  characteristic  size  of 
pores  results  also  in  the  decrease  in  the  effective  surface  of  evaporation. 
Simultaneously  the  saturated  vapor  pressure  of  Mg  decreases  due  to  the  increase  in  the 
meniscus  ciuvnture  in  the  pores.  Both  these  causes  decrease  flow  of  Mg  vapor  to  the 
flame  zone,  resulting  in  the  decrease  in  the  gas-phase  reaction  intensity,  and  therefore, 
to  stability  loss  of  vapor-phase  burning. 

Therefore,  existence  of  feedback  between  the  rate  of  heat  release  and  the 
changes  in  structure  and  properties  of  the  film  at  the  vapor-phase  burning  phase  may 
lead  to  instability,  particularly  to  the  pulsating  mode  observed  experimentally.  The 
feasibility  for  realization  of  this  mode  depends  on  the  state  and  properties  of  the  film 
formed  at  the  pre-ignition  phase.  If  the  oxide  film  fonned  at  the  initial  phases  of  the 
process  is  loose,  only  partial  healing  of  the  pores  may  take  place,  insufficient  for 
cessation  of  vapor-phase  burning.  With  increasing  the  ambient  temperature  the  density 
of  the  original  film  increases.  For  the  reasonably  dense  film  fast  healing  of  the  pores 
and  extinction  of  the  flame  may  occur  after  vapor-phase  ignition.  The  particle  returns  to 
the  slow  heterogeneous  combustion  mode  when  the  characteristic  size  of  the  film  pores 
is  sufficiently  great  and,  therefore,  the  growth  of  the  film  does  not  prevent  from  thermal 
self-acceleration  of  reaction,  resulting  in  a  repeat  flash.  As  this  take  place,  the  higher 
the  particle  temperature  at  the  phase  of  slow  heterogeneous  combustion,  the  shorter  the 
time  before  the  repeat  flash. 
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To  realize  the  pulsating  combustion  mode,  it  is  necessary  tliat  the  activation 
energy  of  the  healing  process  is  higher  than  that  of  chemical  reaction.  In  this  case,  the 
value  of  temperature  always  exists  theoretically  when  the  characteristic  size  of  pores  at 
the  phase  of  slow  heterogeneous  combustion  is  small  to  an  extent  that  thermal  self- 
acceleration  of  reaction  is  not  sufficient  for  equilibrium  breakdown  in  the  system.  As  a 
result,  at  high  ambient  temperatures  the  combustion  mode  may  become  heterogeneous 
again,  as .  we  obseived  in  the  experiments . 

This  work  was  supported  by  the  Russian  Foimdation  for  Basic  Research  (Grant 
95-03-09829). 
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STUDY  OF  DISPERSE  PHASE 
OF  HMX-BASED  ALUMINIZED  PROPELLANTS. 

1.  STEADY-STATE  COMBUSTION  CHARACTERISTICS 

O.G.  GLOTOV,  V  E.  2ARK0,  V.  V.  KARASEV 

Institute  of  Chemical  Kinetics  and  Combustion 
Russian  Academy  of  Sciences,  Novosibirsk  630090,  Russia 

The  efficiency  of  use  of  aluminized  propellants  in  solid  rocket  motors  is  deter¬ 
mined  not  only  by  conversion  of  chemical  energy  in  the  heat  of  combustion  but  also  by 
formation  of  condensed  combustion  products  that  causes  two-phase  loss  of  specific 
impulse,  slagging  in  the  combustion  chamber  and  erosion  of  the  nozzles. 

The  objective  of  this  study  was  to  get  experimental  information  about  generation 
of  condensed  combustion  products  by  HMX-based  model  metalized  propellants  with 
inert  and  energetic  binder  both  in  steady-state  and  nonsteady  regimes  of  combustion 
(the  last  one  is  the  topic  of  companion  paper).  The  data  obtained  by  sampling  tech¬ 
nique  [1-2]  on  size  distribution  in  the  range  from  0.5  pm  to  and  on  content  of 
non-consumed  Al  in  condensed  combustion  products  (CCP)  formed  in  course  of  sta¬ 
tionary  combustion  are  reported  in  present  paper. 

Two  metalized  propellants  were  designated  as  I  and  E  in  accordance  with  the  na¬ 
ture  of  binder.  The  propellants  characteristics  are  presented  in  Table  1 .  The  inert 
binder  consists  of  butadiene  rubber  20%,  transformer  oil  76%,  and  additives  4%  while 
tlie  energetic  binder  consists  of  nitiyl  rubber  plasticized  with  DEGDN  98%  and  addi¬ 
tives  2%.  Propellant  E  contained  numerous  macro  pores  and  thus  had  relatively  low 
density. 

Table  1.  Propellant  formulations  (%  mass.)  and  characteristics 


Pro-.  Energetic 

Inert  Al 

AP 

AP 

HMX 

Density, 

05 

li 

s.. 

pel-  binder 

binder  1)30= 14 

200-315 

5=6000 

>315 

g/cm^ 

r  in  mm/s, 

lant 

-  pm 

pm 

cmVg 

pm 

6<P<64  atm 

I 

15  16.8 

13.1 

31.8 

23.3 

1.7 

r=1.13-P°^^ 

E  20 

20 

18 

7 

35 

1.6 

r=1.88-P°52 

Coarse  size  CCP  particles,  propellant  I 

The  combustion  behavior  of  propellant  I  is  characterized  with  a  weak  agglomera¬ 
tion  ofAl  particles.  At  all  pressures  under  study  (6,  22,  64  atm)  no  more  than  10%  of 
mass  of  sampled  CCP  belong  to  particles  with  size  exceeding  100  pm.  In  fact,  the 
coarse  particles  give  negligibly  small  contribution  mto  total  mass  of  unbumt  Al.  The 
mass  size  distribution  data  are  presented  in  Table  2  in  a  form  of  mean  sizes  cal¬ 
culated  in  the  size  range  120-1000  pm.  The  mean  sizes  decrease  with  pressure.  A 
similar  trend  (Table  2,  2nd  column)  is  revealed  for  incompleteness  of  ^/  combustion 
which  is  determined  in  present  work  as  the  ratio  of  mass  of  unbumt  Al  in  CCI^  parti¬ 
cles  and  mass  of  original  metal  in  the  propellant  sample  rf=  mf  .  Such 
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observations  are  typical  for  aluminized  propellants  combustion  [3],  A  trivial  explana¬ 
tion  of  this  trend  is  in  the  fact  that  the  smaller  time  is  available  for  Al  particles  to  melt 
and  agglomerate  on  the  burning  surface  with  enhanced  burning  rate  under  elevated 
pressure.  A  possible  cause  of  relatively  weak  agglomeration  of  Al  in  combustion  of 
propellant  I  could  specific  properties  of  binder  which  prevents  adhesion  between  metal 
particles  in  preheat  layer.  According  to  the  observations  made  with  high  speed  movie 
and  morphological  analysis  the  aggregates  of  Al  particles  leaving  the  surface  are  only 
partially  sintered  clusters  which  ignite  mostly  at  a  finite  distance  above  the  surface.  It 
should  be  underlined,  however,  that  total  completeness  oiAl  combustion  for  propellant 
I  is  rather  high  due  to  relatively  small  particle  sizes  and  dieir  fast  bum  out  in  flame. 

Coarse  size  CCP  particles,  propellant  E 

An  agglomeration  of  Al  particles  in  combustion  of  propellant  with  energetic  binder 
is  rather  pronounced.  Table  2  presents  characteristic  sizes  Dmn  of  agglomerates  calcu¬ 
lated  in  size  range  120  pm  -  Mass  size  distribution  fimctions  for  CCP  and  for 
unbumt  Al  in  CCP  do  not  essentially  change  in  tlie  pressure  range  1-60  atm,  and  ob¬ 
served  rise  in  D43  is  mainly  due  to  increase  of  number  and  mass  of  particles  larger  than 
1000  pm.  The  porous  structure  of  propellant  E  probably  causes  enhancement  of  ag¬ 
glomeration  due  to  formation  of  Al  particle  clusters  from  neighbor  "pockets"  in  the 
bulk  of  propellant.  Completeness  of  Al  combustion  (Table  2,  2nd  column)  at  low  pres¬ 
sure  is  higher  than  at  moderate  pressures.  One  may  suppose  that  presence  (only  in  case 
of  propellant  E)  of  "matrix"  agglomerates  [4],  the  number  of  which  diminishes  at  ele¬ 
vated  pressures,  is  the  cause  of  more  effective  combustion  of  metal  at  low  pressures 
[1].  The  structure  of  matrix  agglomerate  provides  better  opportunities  for  metal  oxida¬ 
tion  due  to  direct  contact  of  oxidizer  with  neat  metal  at  the  surface  of  agglomerate.  The 
obser\’ed  trend  in  Al  consumption  completeness  can  be  treated  as  indirect  confirmation 
of  hypothesis  that  variation  of  tlie  structure  of  agglomerates  can  make  an  effect  on 
completeness  of  metal  combustion.  Slight  increase  of  metal  combustion  completeness 
at  elevated  pressures  (22  and  62  atm)  can  be  simple  consequence  of  temperature  in¬ 
crease  at  the  burning  surface  and  prolonged  residence  time  for  agglomerates  in  gase¬ 
ous  combustion  products  (due  to  reduction  of  velocity  of  gases  issued  from  the  burning 
surface  with  pressure). 

Fine  CCP  particles,  propellants  I  and  E 

The  typical  feature  of  mass  size  distribution  of  fine  CCP  for  both  propellants  under 
study  is  the  presence  of  local  maximum  (M-peak)  in  size  range  2.4-5  pm  at  all  tested 
combustion  conditions.  We  calculated  mass  size  distribution  of  oxide  particles  assum¬ 
ing  that  all  virgin  Al  particles  were  converted  into  oxide  ones  with  depositing  17%  of 
theoretical  available  oxide  back  on  the  surface  of  Al  particles.  The  results  of  these  cal¬ 
culations  show  that  the  calculated  peak  for  oxide  particles  practically  coincides  with 
main  experimental  peak  in  fine  CCP  distribution  but  its  amplitude  is  much  less  than 
that  of  M-peak  for  CCP  particles.  Thus,  it  should  be  suggested  some  another  (for  ex¬ 
ample,  condensation)  mechanism  for  delivering  large  niunber  of  oxide  particles  into  the 
M-peak. 
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At  elevated  pressures,  an  L-peak  appears  in  the  size  range  0.5-1. 2  pm  with  am- 
phtude  strengthening  with  pressure.  The  estimated  residence  time  for  fine  particles  in¬ 
creases  with  pressure  and  equals  6  ms  at  22  atm  and  10-12  ms  at  62-64  atm  for  both 
propellants.  At  the  same  time,  it  is  known  that  the  longer  residence  time  for  smoke 
(including  submicron)  particles,  the  larger  tlieir  size  [5],  If  one  takes  into  account  that 
sampling  technique  used  permits  to  analyze  particles  exceeding  0.5  pm  in  size,  the  ap¬ 
pearance  of  L-peak  can  be  explained  by  enlargement  of  real  submicron  particles  due  to 
increase  of  their  residence  time  in  flame  witli  pressure. 

To  conclude,  the  results  of  tliis  study  showed  significant  difference  in  metal  com¬ 
bustion  behavior  for  studied  HMX-based  propellants.  Propellant  with  inert  binder  is 
characterized  with  relatively  weak  agglomeration  and  bad  ignition  of  Al  particles. 
Therefore  the  Al  combustion  incompleteness  is  rather  high  at  low  pressures  but  it  is 
very  small  at  elevated  pressures  (^^O.Ol  at  64  atm).  In  case  of  propellant  with  ener¬ 
getic  binder  the  Al  combustion  incompleteness  is  relatively  high  (77=0.25-0.32)  due  to 
large  size  of  agglomerates  formed  and  ignited  at  the  burning  smface.  Some  better 
completeness  of  Al  combustion  at  low  pressures  is  probably  caused  by  presence  of 
“matrix”  agglomerates  with  better  access  of  oxidizer  to  the  surface  of  neat  Al  The 
causes  of  different  combustion  behavior  of  metalized  propellants  are  still  unknown. 

Table  2. 


Incompleteness  ofAl  combustion  rj  and  mean  diameters  (pm)  of  coarse  CCP 
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STUDY  OF  DISPERSE  PHASE 
OF  HMX-BASED  ALUMINIZED  PROPELLANTS*. 

2.  TRANSIENT  COMBUSTION  CHARACTERISTICS 

V.  V  KARASEV,  O.  G.  GLOTOV,  A.  B.  KISKIN,  V  N.  SIMONENKO,  A.  G.  SVIT, 

V.  E.  ZARKO 

Institute  of  Chemical  Kinetics  and  Combustion 
Russian  Academy  of  Sciences,  Novosibirsk  630090,  Russia 

The  objective  of  this  study  was  receiving  information  on  characteristics  of  tran¬ 
sient  combustion  of  HMX-based  model  aluminized  propellants  with  inert  (I)  and  ener¬ 
getic  (E)  binder.  The  experiments  were  perfonned  on  the  measurement  of  recoil  force 
and  visualization  of  the  flame  zone  via  high  speed  photography. 

The  propellant  characteristics  and  experimental  data  on  particles  sampled  in 
course  of  stationary  combustion  are  presented  by  Glotov  0.  G.  and  co-authors  in  com¬ 
panion  paper  in  tliis  proceedings. 

In  order  to  obtain  data  on  transient  burning  rate  response,  the  experiments  were 
conducted  for  both  propellants  under  action  of  sinusoidal  radiant  flux  at  atmospheric 
pressure.  The  set  up  [1]  that  included  computer  controlled  laser  (YAG,  X=1.06  pm, 
120W)  and  microforce  transducer  for  measuring  recoil  force  of  burning  propellant 
sample  was  used.  It  was  found  that  transient  combustion  behavior  of  studied  propel¬ 
lants  is  quite  different  despite  similar  stationary  burning  rate  laws  under  constant  radi¬ 
ant  flux.  The  following  dependencies  of  burning  rate  r  (mm/s)  on  radiant  flux  q 
(cal/(cra^-s))  were  obtained  for  radiation  driven  combustion: 

r=1.44+3.77-10-V'7.03-i0"’V  for  propellant  I  and  r=1.93+4.23-10-V3-65-10V  for 
propellant  E.  Propellant  I  demonstrates  remarkable  response  to  perturbations  of  radiant 
flux  with  resonance  frequency  being  equal  about  2-3  Hz.  For  the  propellant  E,  the  re¬ 
coil  force  signal  has  significant  level  of  noise  (intrinsic  oscillations)  that  makes  it  diffi¬ 
cult  studying  transient  combustion  response.  It  is  assumed  that  in  combustion  behavior 
of  propellant  E  the  important  role  belongs  to  the  presence  of  large  agglomerate  parti¬ 
cles  ignited  at  or  close  to  the  burning  surface.  Detachment  of  agglomerates  from  the 
surface  may  cause  large  scale  perturbations  of  tlie  signal  of  recoil  force. 

The  experiments  on  visualization  were  performed  only  with  propellant  I  at  atmos¬ 
pheric  pressure  using  shadow  image  high  speed  movie  technique  [2].  Combustion  tests 
were  made  under  action  of  constant  and  step-wise  radiant  flux  with  frequency  close  to 
the  resonant  for  tlie  burning  rate  response.  The  cylindrical  propellant  sample  with  di¬ 
ameter  7.5  and  height  7-10  mm  was  tightly  inserted  into  15  mm  long  quartz  tube.  The 
propellant  sample  was  mounted  horizontally  and  visualized  area  located  flush  the  edge 
of  quartz  tube  and  had  tlie  dimensions  6x4  mm^.  The  flash  illumination  frequency  was 
equal  to  2  and  3.4  kHz,  image  magnification  on  the  film  amounted  2.5“^  and  S'*.  In  or- 
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der  to  obtain  data  on  correlated  with  the  phase  of  irradiation  temporal  particle  charac¬ 
teristics  special  design  time-marking  system  has  been  developed.  The  marks  on  the 
fihn  made  by  emitting  diode  allow  to  bring  in  correspondence  the  sequence  of  the  film 
fi*ames  and  temporal  phases  of  radiant  intensity  variation. 

It  should  be  noted  that  the  possibilities  of  studying  the  characteristics  of  two  phase 
flow  in  flame  of  propellant  I  are  veiy  restricted.  The  formed  on  the  surface  particles 
are  the  coral  (grape)  shape  aggregates.  Most  of  recorded  on  the  films  aggregates  has 
characteristic  size  less  than  100  jam  and  about  30%  of  their  number  have  size  less  than 
10-20  jam  which  is  below  resolving  limit  of  high  speed  movie  registration.  The  aggre¬ 
gates  leave  the  burning  surface  non-ignited  and  they  are  ignited  mainly  at  the  distance 
1-5  cm  downstream  in  the  flame.  When  constructing  histograms,  a  diameter  of  aggre¬ 
gates  was  detennined  as  diameter  of  the  circle  witli  the  area  equal  to  particle  image 
area. 

On  tlie  basis  of  film  treatment  the  following  parameters  were  determined: 

Hf  -  averaged  (mean)  number  of  particles  in  a  frame  obtained  during  given  time 
interval  of  histogram  accumulation  (40  ms  in  the  case  of  constant  flux;  25  ms  or  1/8  of 
total  period  of  pulsed  irradiation  in  the  case  of  step-wise  radiant  flux); 

-Oio,  ^>20,  and  D30  -  mean  particle  diameters  in  size  range  0-2 1 0  jmi. 

The  value  of  relative  error  Anc  Irif  was  estimated  as  <0.06  for  experiments  with 
constant  radiant  flux  and  as  <0.02  for  experiments  witli  step-wise  radiant  flux.  The  er¬ 
ror  zlDmn  of  particle  size  determination  was  estimated  as  a  half  of  the  width  of  histo¬ 
gram  size  sub  range  and  equals  (30/2)  pm. 

Tlie  number  of  particles  in  the  individual  frame  on  the  film  is  relatively  high  (about 
15-25  on  each  frame)  and  the  velocities  of  their  motion  are  high  enough  (up  to  7  m/s). 
Therefore  there  is  impossible,  as  a  rule,  to  identity  just  tlie  same  particle  on  neighbor 
frames.  To  ensure  statistic  validi^'  of  experimental  data  we  treated  not  each  frame  but 
every  5th  frame  on  the  film.  Thus,  the  data  on  particles  concentration  and  size  distri¬ 
bution  obtained  via  high  speed  movie  visualization  in  present  work  should  be  treated 
as  local  and  instantaneous^  i.  e.  corresponding  to  the  given  visualized  scene  at  given 
time  interval. 

The  instantaneous  liistograms  can  be  effectively  used  in  evaluation  of  external  ra¬ 
diation  attenuation  by  flame  and  damping  of  acoustic  (pressure)  oscillations  in  rocket 
chamber.  To  construct  liistograms  corresponding  to  volume  flow  rate  of  disperse  phase 
in  flame  one  should  know  the  particle  velocity  field  in  tlie  flame.  The  last  histograms 
can  be  used  in  estimates  for  metal  combustion  completeness  and  for  two  phase  losses 
of  specific  impulse  . 

The  experirnents  with  constant  irradiation  i5r-12  cal/(cm^-s)  showed  substantial 
scatter  in  experimental  data.  However,  only  the  variation  of  “instantaneous”  value  of  «f 
is  larger  tlian  tlie  estimated  experimental  error.  The  variations  of  Wf  and  can  be 
caused  by  heterogeneous  structure  of  propellant  and  relatively  short  time  of  accumula¬ 
tion  of  the  data. 
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In  experiments  with  step-wise  oscillating  radiant  flux  during  first  half  (0,1  s)  of  pe¬ 
riod  the  burning  surface  is  irradiated  with  flux  16  cal/(cm^  s).  During  second  half  of  the 
period  irradiation  intensity  equals  8  cal/(cm^*s).  The  films  under  treatment  contained 
information  about  6  whole  periods.  In  course  of  treatment  each  half  of  the  period  was 
divided  onto  4  equal  portions  and  finally  all  histograms  corresponding  to  appropriate 
(according  to  the  current  phase  of  irradiation)  portions  were  summarized.  This  was 
done  aiming  to  diminish  “statistical”  noise  due  to  the  effect  of  occasional  micro  non¬ 
uniformity  of  composition  and  structure  of  the  propellant  sample.  The  results  for  the 
combustion  under  step-wise  oscillating  irradiation  are  presented  in  Figure.  To  make  it 
obvious  every  8  points  on  the  plots  in  Figure  are  repeated  in  a  cyclic  manner  in  accor¬ 
dance  with  the  process  repetition  frequency  5  Hz.  Each  point  on  the  plots  is  the  result 
of  treatment  of  2100-2600  particle  images  that  makes  it  possible  to  assume  that  the 
effect  of  micro  compositional  heterogeneity  of  the  propellant  sample  can  be  neglected. 

The  results  of  tliis  study  have  shown:  a)  There  is  distinguished  response  of  dis¬ 
perse  phase  parameters  to  variations  of  the  external  radiant  flux  but  temporal  variations 
of  the  parameters  exceed  remarkably  the  experimental  error  only  for  «f;  b)  Comparison 
of  the  Wf  response  witli  recoil  force  response  at  the  same  radiant  flux  oscillations  al¬ 
lows  to  conclude  that  the  value  of  «f  is  proportional  to  the  burning  rate  r  if  one  assume 
tliat  the  recoil  force  c)  The  mean  value  of  under  step-wise  irradiation  by  1.5 

times  higher  then  that  un¬ 
der  constant  irradiation  de¬ 
spite  the  same  mean  energy 
input  to  the  burning  sui- 
face.  Study  of  tliis  effect 
could  be  interesting  and 
important  from  the  point  of 
view  of  suppression  of 
acoustic  pressure  oscilla¬ 
tions  in  rocket  chamber. 


Figure.  Temporal  de¬ 
pendencies  of  particle 
characteristics  in  flame 
zone  of  propellant  I  com¬ 
busted  under  step-wise  ra¬ 
diant  flux,  pressure  1  atm. 
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COMBUSTION  PRODUCTS  OF  GAS  MIXTURES 
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Now  in  different  industry  branches  there  are  technological  processes  which  at 
possible  accidental  regimes  are  followed  by  the  fast  energy  release.  It  causes  the 
overloading  of  the  construction  and  under  certain  conditions  actually  destroys  tiiem 
[1].  Because  these  processes  often  realizes  in  large-scale  or  sometime  even  unique 
installations,  the  scale  model  experiments  are  used  widely  to  estimate  the  strength  of 
these  objects.  And  the  accidental  energy  releasing  is  simulated  usually  by  the 
detonation  of  solid  high  explosives  (HE)  [1,2].  As  shown  [3],  the  accidental  energy 
releasing  simulation  by  solid  HE  detonation  in  the  scale  model  experiments  in  many 
cases  results  some  uncertainties  and  does  not  guarantee  correct  conservative 
evaluation  of  the  adequate  response  of  a  full-scale  objects  to  a  disruptive  accident.  To 
increase  the  credibility  to  a  strength  prediction  ,  the  preliminary  calculated  accidental 
energy  release  must  be  reproduced  experimentally  with  adequate  parameters. 

In  this  report,  we  briefly  describe  the  resiftts  of  experimental  investigations  of 
closed  volume  combustion  of  the  condensed  explosive  materials  (EM)  ignited  by  the 
hot  products  of  hydrogen-air  mixtures  combustion  and  show  the  possibility  of 
developing  the  simulator  which  provides  the  monotone  energy  releasing  with 
maximum  combustion  products  pressure  ranging  to  100  MPa  through  the  time  of  20 
to  400rasec.  These  parameters  exceed  these  ones  of  all  known  load  simulators  both  in 
duration  of  process  and  in  amplitude  value  [1,4]. 

The  experiments  were  carried  out  in  the  cylindrical  steel  combustion  chamber 
with  inner  volume  of  0.012  m^.  The  used  explosive  materials  were  the  pyroxylin  gun 
powder  PI  5/1  and  gelatin  explosive  (GE)  based  on  finely  divided  penthrite.  Powder 
charge  was  an  assembly  of  tubular  powder  elements  every  of  which  was  0.35  m  long. 
The  mass  of  charge  was  varied  from  0.15  kg  to  0.5  kg.  The  GE  charge  was  formed  as 
a  hollow  cylinder,  its  mass  of  the  same  range  was  varied  by  the  thickness  of  the 
explosive  layer  and  the  height  of  cylinder. 

The  EM  charges  were  ignited  by  the  hot  products  of  combustion  of  10% 
hydrogen-airmixture  at  initial  pressure  of  3  MPa.  The  technique  of  mixing  the  gases“ 
was  the  same  which  described  elsewhere  [5].  The  mixture  was  fired  in  the  center  of 
the  chamber  by  explosion  of  the  thin  Ni-Gr  alloy  wire  by  the  capacitor  discharging 
[6]. 

In  experiments  the  pressure  and  the  temperature  of  gas  in  combustion  chamber 
and  the  heating  of  its  wall  were  measured.  From  the  pressure  records  the  maximinn 
combustion  products  pressure  Pk  and  the  time  of  EM  charge  combustion  Tf  (equals  t 
time  of  maximum  pressure  realizing,  see  Fig,  1)  were  determined. 
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Experimental  values  of  peak  pressure  Pk 
and  temperature  Tk  were  compared  with  the 
results  of  calculations  of  equilibrium  parameters 
of  combustion  products  for  used  EM.  The 
possibility  of  21  components  presence  in  the 
combustion  products  has  been  taken  into  account. 
And  it  has  been  assumed  that  the  explosive 
materials  are  in  closed  volume  with  heat  insulated 
wall  and  the  gas  velocity  equals  zero  in  the 
chamber.  The  program  of  calculation  includes  the 
solution  of  simultaneous  equations  set,  namely 
matter  balance  equation,  heat  balance  equation, 
density  constancy  equation  and  equilibrium 
equation. 

Charge  EM  combustion  time  (or  in  other 
words  time  of  pressure  increasing  to  maximum 
value)  is  determined  by  the  type  of  EM  and 
charging  density,  which  defined  as  a  charge  mass 
referred  to  volume  of  combustion  chamber. 
Increasing  of  charge  density  reduces  the  combustion  time  for  both  type  of  EM 
through  the  range  investigated.  Such  reducing  is  caused  by  the  increasing  of  normal 
combustion  rate  with  pressure  increasing  [7]. 

As  shown  in  Fig.  3  in  the  equal  conditions  GE  charge  provides  increasing  of  the 
combustion  time  because  the  normal  combustion  rate  of  gun  powder  is  much  higher 
than  that  of  GE  [7]. 

The  parameter  which  defines  the  duration  of  combustion  is  the  characteristic 
geometrical  size  of  single  powder  element  or  the  thickness  of  GE  layer  /.  Production 
process  of  gelatin  explosives  enables  to  product  the  plates  with  wide  range  of 
thickness.  Thus  the  variations  of  combustion  time  fi-om  40  ms  to  400  ms  are  possible 
through  the  wide  range  of  charging  density  using  the  GE  charge. 

The  value  of  the  peak  combustion  pressure  is  determined  by  the  parameters  of 
hydrogen-air  mixture  and  the  charging  density.  Using  the  10%  hydrogen-air  mixture 
with  initial  pressure  of  3  MPa  the  variations  of  charging  density  up  to  65  kg/m  allow 
to  obtain  the  peak  pressure  in  combustion  chamber  from  15  MPa  to  100  MPa  as  for 
gun  powder  PI  5/1  and  GE  (see  Fig.2). 

For  powder  charges  through  the  whole  parameters  range  investigated  the 
differences  between  experimental  and  pre-calculated  values  of  peak  pressure  are 
within  error  limits. 

Thus  the  state  of  products  of  combustion  of  condensed  EM  ignited  by  the 
hydrogen-air  mixture  combustion  products  may  be  considered  as  the  thermodynamic 
equilibrium  state,  and  one  can  use  the  pre-calculated  data  to  predict  the  impact  load. 

The  level  of  peak  pressure  values  realized  both  for  gun  powder  charge  and  for 
GE  charge  is  greater  than  the  arithmetic  sum  of  pressures  caused  by  separate 
combustion  of  each  components  (see  Fig,2).  Considering  the  analysis  of  combustion 


Fig.  1.  CP  pressure  V(t)  for 
charging  density  variation 

1-  gelatine  explosive,  p=5kg/m^ 

2-  gun  powder,  p=25kg/m^ 

3-  gun  powder,  p=41.7kgW 
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products  one  can  conclude  that  the  noticed 
excess  of  pressure  is  caused  by  the  after¬ 
burning  of  CO  to  CO2  in  the  unusable  oxygen 
in  the  chamber  after  the  hydrogen-air 
combustion  has  been  completed.  For  gun 
powder  charge  the  after-burning  contribution 
to  increasing  of  pressure  is  greater  as  the 
charging  density  is  higher  up  to  20  kg/m  ,  and 
for  GE  charge  -  up  to  12  kg/m  . 

Thus,  the  results  presented  show  that  the 
closed  volume  combustion  of  such  explosive 
materials  as  gun  powder  and  gelatin  explosive 
ignited  by  the  hot  products  of  hydrogen-air 
mixture  combustion  is  the  stationary  process  in 
charging  density  range  up  to  65  kg/m  .  This 
energy  source  allows  to  produce  the  monotone 
increasing  of  load  with  predicted  values  of 
pressure  in  range  from  15  to  100  MPa.  And 
the  time  of  realizing  the  peak  pressure  may  be 
varied  from  40  to  400  ms  by  monitoring  both  choosing  the  type  of  EM  and  the 
thickness  of  EM  charge. 
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Fig.2.Maximum  combustion  products 
pressure  pKVs.charging  density. 

GE:  1-  theory,  A  -  experiment 

powder:  2-theory,  •  -  experiment 

3-  only  powder  combustion 

4-  arithmetic  sum  of  pressures 
caused  by  separate  combustion  of  gas 
and  powder  components. 
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INFLUENCE  OF  SOLID  PROPELLANT  STRUCTURE  ON  METAL 
AGGLOMERATION. 

Babuk,  V.A.,  Vasilyev,  V.A.,  Sviridov,  V.  V 
Baltic  State  Technical  University 
lst.Krasnoarmeyskayastr.,L  Saint  Petersburg,  198005,  RUSSIA 

Meta!  agglomeration  in  solid  propellant  combustion  is  complex  and  multi-stage 
process.  Agglomeration  includes  a  number  of  phenomena.  Necessary  condition  for  the 
agglomeration  is  foiTnation  so-named  “frame  layer”  of  the  propellant  burning  surface 
[1],  Frame  layer  is  porous  structure  formed  with  metal,  metal  oxide,  and  carbon. 

Structure  of  composite  solid  propellant  has  an  sufficient  influence  on  “frame 
layer”  formation.  Change  of  propellant  structure  leads  to  change  in  agglomeration  and, 
consequently,  in  parameters  of  condensed  phase  combustion  products  (CPCP). 

Propellant  structure,  can  be  characterized  with  the  help  of  the  terms  “pocket”  and 
“interpocket  bridge”.  “Pocket”  is  a  part  of  binder  surrounded  by  oxidizer  particles. 
“Interpocket  bridge”  is  a  part  of  binder  located  between  two  neighbouring  oxidizer 
particles  which  take  part  in  “pocket”  formation. 

“Frame  layer”  appears  in  “pocket”  burning.  Binder  of  “interpocket  bridges”  does 
not  take  part  in  “frame  layer”  fonnation.  Hence,  only  metal  of  “pocket”  is  involved  in 
agglomeration.  So,  characteristics  of  the  “pockets”  have  strong  influence  on 
parameters  of  agglomerates  [2].  To  verify  an  influence  of  propellant  structure  on  metal 
agglomeration,  experimental  investigation  was  carried  out. 

Particles  of  CPCP  were  quench-collected  and  analyzed.  Also,  experimental  data 
were  compared  with  results  of  numerical  modeling  of  propellant  structure.  Computer 
calculations  were  made  with  the  help  of  propellant  structure  model  developed  in  [3]. 
For  computer  modeling,  a  method  of  “statistical  tests”  was  used.  Parameters  of 
propellant  structure  were  calculated  on  the  base  of  information  about  distribution 
function  for  oxidizer  particle  centers  in  propellant  and  size  distribution  function  for 
the  oxidizer  particles. 

Propellants  with  isoprene  rubber  as  binder  and  ammonium  perchlorate  (AP)  as 
oxidizer  were  studied.  Propellants  were  differed  with  oxidizer  dispersity  only.  Mass 
content  of  components  in  the  propellants  is  shown  in  Table  1.  Propellants  includes 
aluminum  powder  of  up  to  50  fim  particle  diameter.  Four  fraction  of  AP  particles  were 
used  for  propellant  composition. 

Fraction  #1 ;  particle  size  is  of  1  /im. 

Fraction  #2:  particle  size  is  less  than  50  ^m. 

Fraction  #3:  particle  size  is  within  160-315  pm. 

Fraction  #4;  particle  size  is  within  400-600  pm. 

Proportions  of  tlie  AP  fractions  in  the  propellants  are  shown  in  Table  2. 
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Table  1 


Components 

Aluminum 

A1 

Ammonium 

Perchlorate 

NH4CIO4 

Isoprene 

Rubber 

-CHj-  C  =  CH- 
1 

C 

Additives 

Oil 

C19H35 

Contain,  % 

24 

64 

2.6 

0.4 

9 

Table  2 


Propellant 

#1 

#2 

#3 

#4 

#5 

Proportion  of  AP 
fraction 

Fr.#l-Fr.#2 

50%-50% 

Fr.#2 

100% 

Fr.#2-Fr.#3 

40%-60% 

Fr.#3 

100% 

Fr.#4 

100% 

Experiments  consist  in  combustion  of  small  propellant  samples  in  constant 
volume  bomb  (CVB)  filled  with  inert  gas.  Burning  agglomerates  were  quenched  and 
collected.  Distribution  functions  and  chemical  compositions  for  samples  of  quench- 
collected  agglomerates  were  determined. 

Two  main  parameters  characterize  agglomeration  process; 

1 .  a  part  of  metal  taking  part  in  agglomeration 

1.  degree  of  integration  of  metal  particles  (increase  of  particle  size  due  to 
agglomeration). 

Particle  integration  degree  can  be  characterized  by  mass-medium  diameter  D43  of 
agglomerates.  Influence  of  propellant  structure  on  agglomeration  is  well  seen  at  low 
pressure  (  about  1  MPa).  Increase  of  AP  dispersity  in  propellant  first  (from  Propellant 
#5  to  Propellant  #2)  leads  to  decrease  in  and  D43,  and  then  (from  Propellant  #2  to 
Propellant  #1)  these  parameters  increase.  Nature  of  these  dependencies  is  considered 
to  connected  with  propellant  structures.  Decrease  of  oxidizer  particle  sizes  causes 
increase  in  part  of  “interpocket  bridges”  in  comparison  with  part  of  “pockets”  in 
propellant.  This  lead  to  decrease  of  when  moving  from  Propellant  #5  to 
Propellant  #2  (Fig.l).  In  the  case  of  Propellant  #1,  AP  particles  are  supposed  to  be  not 
able  to  form  “pocket”  due  to  their  small  sizes.  So,  structure  of  Propellant  #1  is  similar 
to  homogeneous  propellant  structure.  (Homogeneous  propellant  can  be  characterized 
as  propellant  consisted  of  one  big  “pocket”).  That  is  why  part  of  metal  taking  part  in 
agglomeration  (Zm**)  for  Propellant  #1  at  low  pressure  is  more  than  the  same  for 
Propellant  #2  (Fig.  1). 

Used  model  of  propellant  structure  allows  determine  parameters  and  quantity  of 
“pockets”  in  propellant.  Pocket  sizes  have  influence  on  agglomeratie  sizes,  but 
“pocket”  model  can  be  used  for  agglomerate  size  distribution  calculations  in  the  case 
of  so-named  “pocket  mode”  of  agglomeration.  Agglomerate  dispersity  depends  also  on 
propellant  burning  rate,  “frame  layer”  features  and  others  factors.  As  a  result, 
agglomeration  can  go  under  so-named  “interpocket”  and  “prepocket”  regimes. 
However,  parameters  of  “pockets”  make  it  possible  to  evaluate  agglomerate  sizes. 


0  100  200  300  400  500 

D^3,  fim 

Figure  1.  Dependence  of  Zn,^  on  mass>medium  diameter  of  AP  particles  in 
propellants. 

Table  3  shows  parts  of  “pockets”  in  binder  of  propellants  (gp)  calculated  with  the 
help  of  propellant  structure  model.  These  data  confirm  dependence  of  7^  on 
propellant  structure.  Table  3  shows  also  mass-medium  diameters  for  agglomerates 
D43®  and  for  “pockets”  in  propellants  043^  calculated  with  the  help  of  propellant 
structure  model.  Both  experimental  investigation  and  calculations  show  that  043^ 
decrease  when  moving  from  Propellant  #5  to  Propellant  #2  and  then  D43®  increase  for 
Propellant  #1  in  comparison  with  Propellant  #2. 

Table  3 


So,  the  result  of  this  work  allow  to  draw  the  conclusion  that  developed  propellant 
structure  model  can  be  used  for  evaluation  of  agglomeration  characteristics. 

Experimental  investigation  was  supported  by  United  States  Air  Force 
European  Office  of  Aerospace  Research  and  Development  (contract  F6 1708-96- 
W0269) 
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MODEL  OF  ALUMINUM  PARTICLE  IGNITION  IN  SOLID  ROCKET 
PROPELLANT  COMBUSTION 

Babuk  V.A.,  Naslednikov P.A. 

Baltic  State  Technical  University 
Ist.Krasnoarmeyskaya  str.,1.  Saint  Petersburg,  198005,  RUSSIA 


Use  of  AI  as  a  metal  fuel  is  characteristic  for  modem  solid  rocket  propel- 
lants(SRP).  Essential  feature  of  aluminised  SRP  combustion  is  agglomeration 
(merging)  of  metal  particles  in  surface  layer  of  burning  propellant  that  leads  to  a  num¬ 
ber  of  negative  consequences.  When  modelling  agglomeration  phenomenon  it  is 
needed  to  consider  various  stages  of  metal  fuel  evolution  in  comWstion  process  in¬ 
cluding  ignition  of  Al  particles  in  surface  layer.  Knowing  domain  of  Al  particles  igni¬ 
tion  it  is  possible  to  describe  conditions  in  which  heterogeneous  combustion  of  parti¬ 
cles  takes  place  and  hence,  one  can  predict  fraction  of  Al  oxide  in  agglomerats  being 
formed  [2].  Model  of  ignition  presented  permits  to  take  into  account  influence  of 
composition  and  stmcture  of  SRP  and  also  of  pressure  in  combustion  chamber 
ofSPRE. 

The  basis  for  model  developed  is  as  follows. 

When  Al  particle  appears  in  surface  layer,  its  oxide  film  gains  some  properties  that 
remain  invariable  up  to  particle  ignition.  Initial  oxide  film  has  crystal  structure  and 
performs  protective  function  (with  respect  to  oxidation  process).  Because  of  high  rate 
of  heating  it  cracks  under  thermal  stress.  Oxidation  of  metal  taking  place  around  the 
cracks  leads  to  formation  of  amorphous  oxide.  Such  an  oxide  does  not  perform  pro¬ 
tective  function  and  crumbles  under  the  action  of  thermal  stress. 

All  particles  are  equidistant  from  upper  bound  of  "frame  layer"  [2].  They  are  subjected 
to  identical  conditions,  i.e.  heterogeneity  of  propellant  structure  and  hence  nonhomo¬ 
geneity  of  surface  layer  hadn't  been  taken  into  account. 

Consideration  of  all  processes  had  been  realized  applying  to  "averaged  cell"  which 
is  burning  correspondingly  to  burning  of  the  bulk  of  propellant  [3]. 

When  "averaged  cell"  burning  variations  of  temperature  and  concentration  of  com¬ 
bustion  products  are  assumed  to  be  one-dimensional. 

Properties  of  gaseous  phase  depend  on  point  position  respective  to  gasification 
surface  of  binder  and  oxidizer.  The  availability  of  relatively  large  heterogeneity’s 
(such  as  dispersed  oxidizer)  in  propellant  structure  leads  to  production  of  highly  tur- 
bulized  jets  from  binder  and  oxidizer.  In  the  context  of  scheme  assumed,  concentra¬ 
tion  is  linear  with  respect  to  coordinate  and  it  can  be  determined  by  calculation  of 
mixing  interval  length  [1,5],  At  the  end  of  this  interval  gas  composition  is  close  to 
equilibrium  which  can  be  determined  applying  methods  of  equilibrium  thermodynam¬ 
ics  to  mixture  consisting  of  oxidizer  and  binder  taken  in  ratio  which  is  characteristic 
for  the  "pocket". 
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It  is  diffusion  stage  that  limits  oxidation  rate  of  AJ  (but  during  model  trials  also  dif¬ 
ferent  variants  were  considered  [6]). 

Analysis  of  ignition  possibility  is  carried  out  having  parameters  of  media  and  ones 
of  particles  fixed  at  levels  corresponding  to  a  certain  position  in  a  surface  layer. 

Temperature  distribution  in  particle's  cross  section  is  absent. 

Change  of  mass  of  metal  owing  to  chemical  interaction  is  not  considered. 

Convective  heat  exchange  is  determined  by  characteristics  of  "fi-ame  layer"  as  a  po¬ 
rous  structure  and  by  A1  particle's  size  [4], 

The  media's  oxydizing  components  are  :  H2O,  CO2. 

The  basis  of  the  problem's  mathematical  description  is  heat  balance  equation. 
When  surface  temperature  is  Ts  and  media  temperature  is  Tg^  ignition  condition  is 
equality  between  feeding  heat  current  (Q+)  and  draining  heat  current  (Q_)  and  also 
excess  of  (Q+)  over  (Q_)  in  a  given  temperature  vicinity. 

For  different  SRP  classes,  ignition  places  in  "frame  layer"  are  different  [2]. 
This  is  because  of  difference  in  ignition  process  conditions.  Analysis  of  the  model 
permits  to  determine  ignition  temperature  values  (in  interval  from  900K  to  HOOK)  for 
a  number  of  SRP  compositions.  These  values  are  in  reasonably  good  agreement  with 
experimental  data. 

The  model  gives  possibility  to  reconstruct  a  number  of  regularities  that  have 
been  obtained  experimentally.  Analysis  of  Tbm  dependence  on  a  number  of  parameters 
was  carried  out  for  one  particular  SRP  composition  belonging  to  class  B  according  to 
terminology  of  article  [2]  (fig.l).  It  can  be  seen  that  the  difference  Tbm(P)-Ts(P) 
does  not  exceed  lOOK  and  in  fact  Tbm  does  not  change  when  P>3MPa.  Like  in  ex¬ 
periments,  the  place  where  ignition  occurs  coincides  with  surface  of  burning  propel¬ 
lant  (i.e.  upper  surface  of  "frame  layer"). 


T,  K 


Fig.l  Temperature  of  particle  Ts  and  that  of  media  Tbm  as  functions  of  combustion 
chamber  pressure  P  during  ignition. 
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Usage  of  highly  dispersed  A1  powder  (d<lmkm)  leads  to  fast  decrease  of  Tbm  (fig.2). 


d,  mkm 

Fig, 2  Temperature  of  media  Tbm  as  function  of  A1  particle's  diameter  during  ignition. 

Thus  the  model  of  ignition  developed  meets  experimental  data  in  quantitative 
and  qualitative  sense. 

The  authors  hope  to  apply  tlie  created  model  to  method  of  prediction  of  ag¬ 
glomeration  process  characteristics  for  different  SRP. 
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STUDY  OF  HIGHLY  DISPERSED  OXIDE  FORMATION  IN 
ALUMINUM  PARTICLE  COMBUSTION 

V.A.Babuk,  M.S.Malakhov 
Baltic  State  Technical  University 

l-stKrasnoarmeyskayastr.,  1,  Saint  Petersburg,  198005,  RUSSIA 

As  we  know,  aluminum  powder  is  an  organic  part  of  many  types  of 
composite  solid  propellants.  It  is  established  that  in  the  process  of  burning 
propellant  agglomeration  has  been  put  into  effect  on  its  surface  layer.  Some 
particles  of  metal  flow  together  into  drops,  exceeding  by  one  order  ie  size  of 
initial  particles.  The  rest  of  the  particles  are  burnt  down  in  the  smface  layer  of 
burning  propellant. 

It  is  accepted  to  divide  the  particles  of  conditioning  phase  of  products  of 
combustion  of  aluminum  into  2  fractions:  the  first  one  -  agglomerates  (particles 
with  sizes  more  than  30-40  pm)  and  the  second  are  the  particles  of  Highly 
Dispersed  Oxide  (HDO)  with  sizes  less  than  30-40  pm. 

Supposedly,  the  particles  of  three  fractions  will  have  to  do  with  HDO, 
which  differ  from  each  other  by  way  of  formation.  The  first  fraction  is 
represented  by  the  particles  formed  as  a  result  of  vapor-phase  combustion  of 
initial  nonagglomerated  particles.  The  second  fraction  -  the  particles  of  oxide  as 
a  result  of  initial  burning  of  nonagglomerated  particles  in  heterogeneous  mode. 
The  third  is  represented  by  particles  which  are  the  result  of  vapor-phase  burning 
of  agglomerates. 

In  this  way,  the  preliminary  analysis  shows,  that  the  flmction  of  size 
distribution  of  HDO  will  have  complex  polimodal  character. 

Authors  of  this  paper  created  the  method  of  research  of  HDO  particles, 
the  description  of  which  is  in  the  work  [1].  The  method  described  serves  for 
direct  determination  of  parameters  of  HDO  particles  and  realizes  with  the  help 
of  contact  research  method.  Due  to  this  method  some  dispersible  characteristics, 
chemical  composition  of  HDO  particles  may  be  determined;  it  is  possible  to 
research  selected  particles  with  the  use  of  a  microscope. 

With  the  help  of  this  method  an  experimental  investigation  on 
determination  of  parameters  of  HDO  particles  near  the  surface  of  burning 
propellant  has  been  held.  The  HDO  particles  of  the  third  fraction  are  represented 
by  less  part  (5  -  20%)  from  all  HDO  particles  selected  because  the  process  of 
burning  of  agglomerates  stopped  near  the  surface  of  propellant  burning 
(agglomerates  were  quenched  with  inert  medium  at  the  distance  6  mm  from 
burmng  surface).  So,  the  results  of  this  work  describe  the  process  of  HDO 
formation  in  the  process  of  combustion  nonagglomerated  metal. 

Several  propellant  compositions  were  studied  which  characterized  by 
different  oxidizer  dispersion  and,  as  a  result,  difference  in  agglomerate  sizes  and 
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agglomeration  degree,  that  should  lead  to  different  correlation  between  HDO 
particles  of  three  different  fractions. 

It  is  important  to  mention  that  HDO  particles  stick  together  actively  after 
selection,  forming  coagulants.  Their  definition  was  realized  by  dispersion  in  a 
liquid  with  a  help  of  ultrasound.  We  can  not  confirm,  that  all  coagulation  were 
destroyed  and  there  were  no  repeated  coagulation.  Thus,  in  that  way  it  is 
necessary  to  count  that  all  got  parameters  of  the  particles  can  be  raised  to  high. 
Figure  1  shows  distribution  ftmctions  of  HDO  particles  for  different  levels 
of  pressure,  built  on  the  base  of  experimentally  gotten  histogrammes. 
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Fig.l.  Functions  of  density  of  HDO  particles  distribution  for  the  different 
compositions  gotten  in  the  conditions  of  different  levels  of  pressure. 

Composition  N1  -  oxidizer  dispersity:  50%  ~1  pm  ,  50%  less  than  50  pm 
K22  -  100%  less  than  50  pm 
:N25-  100%  400-600mkm 

As  follows  from  all  the  results,  there  are  two  main  fractions;  in  the 
intervals  20-200  mpm  and  300-1000  mpm.  Supposedly  the  first  is  represented 
by  the  particles  of  the  lesser  part  fraction,  formed  as  a  result  of  vapor-phase 
burning  of  nonagglomerated  particles.  These  particles  -  10  %  by  mass  from  all 
selected  HDO  particles. 

Graphically,  the  dependence  of  mass-medium  diameter  of  HDO  particles 
from  the  dispersion  of  oxidizer  is  represented  in  Figure  2. 


Fig  2.  Dependence  of  mass-medium  diameter  of  HDO  particles  on  oxidizer 

dispersity. 
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It  is  necessary  to  notice  that  the  particles  of  the  first  fi'action  will  be  the 
lesser  parts  (the  smallest),  because  the  size  of  the  initial  metal  particles  is  much 
lesser  than  of  agglomerates  size.  Supposedly  their  sizes  are  situated  in  interval 
10-200  mpm. 

Using  this  method,  we  can  affirm  that  the  losses  of  particles  with  sizes 
less  than  300  mpm  in  size  are  very  little.  A  number  of  experiments  took  place  on 
determination  of  mass  imposing  of  HDO  particles  selection  for  control 
compositions  and  it  was  established  that  70-85%  particles  have  been  selected. 
In  this  case  the  losses  of  particles  are  more  for  the  composition  with  less  level  of 
the  agglomeration,  thus  is  more  the  share  of  products  of  combustion  of 
nonagglomerated  particles.  So,  we  can  propose,  that  the  part  of  submicron 
particles  of  HDO  in  size  10-200  mpm  is  more  than  it  considered  to  be  earlier 
and  may  consist  of  20-50%  fi-om  the  mass  of  all  HDO  particles,  forming  in  the 
process  of  combustion  of  nonagglomerated  metal.  Further  investigation  of 
conformities  to  natural  laws  foundation  of  HDO  must  be  directed  to  verify  the 
data  proposal. 

This  investigation  was  supported  by  United  States  Air  Force 
European  Office  of  Aerospace  Research  and  Development  (contract 
F61708-96-W0269). 
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AERODYNAMICS  OF  A  CLOUD  OF  HOT  SOLID  AEROSOL 
PARTICLES 

S.l.Igolkin 

Baltic  State  Technical  University 

The  physics  model  and  the  expressions  are  supposed  for  the  parameters 
evaluation  of  anomalous  motion  of  hot  aerosol  particles  cloud  in  atmosphere  and  air 
flows.  The  equations  are  obtained  related  to  the  velocity  of  the  cloud,  geometry, 
concentration  of  tire  particles  and  tlieir  temperature. 

Physics  model  and  equation  of  motion 

Let  suppose  the  hot  aerosol  particles  are  confined  in  the  sphere,  having  cross 
section  So.  The  particles  are  faintly  connected  each  other,  for  example,  by  the 
electrostatic  forces.  They  have  non-spherical  shape  and  their  concentration  is  not  too 
high  to  suppose  the  cloud  to  be  blow  up  by  air  jets  and  flows.  The  air  penetrates 
tlirough  the  front  boundary  of  the  cloud  having  Vo  velocity.  It  is  heated  by  particles 
when  moves  inside  the  cloud,  expands  and  leaves  it  at  Vi. 


Follow  equations  are  hold  true  in  a  fi’ame  of  reference  related  to  the  cloud. 
Flow  continuity  equation: 

Po^o‘^0  ~  ~  G  ,  (1) 

Equation  of  gas  momentum; 

G  =  -  F,  +  G  K,  ,  (2) 

Equation  of  rectilinear  motion  for  particles,  being  expressed  by  the  balance  of 
reactive  driving  force  Fr  and  aerodynamic  drug  force  Fcx: 
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pR  =  Fc,  , 


(3) 


or 

WS,(V,-V,)  =  C,S,S}^,  (4) 

Here 

K=(K, +  rj/2,  (5) 

Cx  -  average  drug  coefficient, 

Sp  -  summary  cross  section  of  all  particles. 

At  constant  pressure  flow  continuity  equation  leads  to  follow  relation: 

VJV,  =  TJT^  .  (6) 

Taking  it  into  accoimt  and  transformed  (4)  we  get; 


45,  t;  -  t; 


(7) 


The  relationship  of  the  particles  parameters  and  heated  air  temperature  can  be 
received  from  the  heat  transfer  equations  as  well. 

Let’s  introduce  tlie  follow  meanings. 

(X  —  heat-transfer  coefficient; 

Spp  -  summary  surface  of  all  particles; 

Tp  —  surface  temperature; 

r  =  (7;  +  r„)/2-  average  temperature  of  heated  air. 

In  the  assumption  of  constant  temperature  of  the  particles  heat  transfer  equation 
is  written  as  follow: 


aS,,{T,-r)  =  (T,-T,)p,V,S,C, 


(8) 


or: 

y  -  ^  ^PP  ~  7" 

^  CpP,  S,  T,-T, 


(9) 
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Excluding  from  ( 7  )  and  ( 9 )  the  unknown  temperatures  we  get  the  follow 
expression  for  aerosol  cloud  motion  through  undisturbed  air: 


2a  Sp, 


pp 


^pPo^x  ^P 


T 

Tp 


Cy 


-)~1 


(10) 


The  exit  temperature  of  air  is  determined  by  (7)  and  the  velocity  of  heated  air 
passed  through  the  cloud  is  given  by  (6).  The  real  parameters  of  the  particles  Cx,  Sp, 
Spp  have  to  be  set  or  determined  in  a  special  investigation.  All  these  parameters  are 
completely  dependent  on  the  way  of  aerosol  formation. 

Thus,  at  some  combinations  of  sizes,  shapes,  concentration  and  temperature  of 
the  dispersed  particles  the  motion  of  aerosol  cloud  in  undisturbed  air,  in  the  air  flows 
and  jets  can  take  place.  This  result  give  us  satisfy  explanation  of  unpredictable 
behaviour  of  active,  hot,  for  example  condensing  and  relaxing,  two-phases  cloud 
formations. 

The  energy  source  for  such  motion  is  the  temperature  drop  caused  to  the  heating 
of  the  gas  by  hot  particles.  The  nature  of  driving  force  affected  on  the  aerosol  cloud  is 
the  same  as  one  to  use  in  a  classic  ramjet  engine. 
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THE  ROLE  OF  CONVECTION  AND  HEAT 
EXCHANGE  BY  RADIATION  IN  PROCESSES  OF 
HETEROGENEOUS  BURNING  OF  A  PARTICLE 

KALINCHAK  V.  V'.,  ORLOVSKAYA  S.  G.,  PRUDNIKOVA  J.  V., 
LEKHNO  D,  A. 

Odessa  State  University  by  L  L  Medmikov 

The  burning  of  dispersible  fuels  in  fire-chamber  of  torch  type  or  with 
a  circulating  boifing  layer  occurs  at  intensive  heat  and  mass  exchange  of 
particles  with  gas  phase.  The  role  of  the  exchange  can  change  with  changes 
of  the  diameter  and  the  velocity  of  particles.  An  influence  of  natural  and 
induced  convection  on  steady  (burning)  and  critical  (ignition  and  extinction) 
hightemperature  states  of  carbon  particle  taking  into  account  its  heat  transfer 
by  radiation  to  cold  reaction  chamber  walls  has  been  studied  in  the  work. 

Taking  into  account  that  the  value  of  the  Bio’s  criterion  C  1  the  change 
of  particle  temperature  in  time  is  expressed  by  equation 

1 

— CiPid  “  QcA  ~  Qm.r.  ~  (Jn  7i(t  =  0)  =  /Ito, 

I'he  density  ol‘  chemical  heat  secretion  is  determined  in  the  main  by 
cheniical  reactions  C  +  O2  — ^  CO2  (I)  and  2C  +  O2  — ^  200  (II) 

(Ich  ~  {kiQ I  +  kiQ^)  , 

Qmx  —  OfiT]  —  T2)  —  is  tile  density  of  the  heat  flow  by  molecular  and 
convective  mechanism, 

qr  ~  ecr  is  the  density  of  the  heat  flow  by  radiation; 

where  <x  X2'  Nu/d,  ft  =  Nu/d,  D  -  X-if  {c,2pi) , 

Nu  =  2A-  b\/We.  •  b  =  0,65,  Re  =  — ,  U2  =  Pr  ■  D, 

1^2 

The  increase  of  the  particle  diameter  depends  on  the  mass  transfer  coefFici- 
ent  /?  and  constants  of  rates  of  chemical  reactions  (I)  and  (II)  in  the  explicit 
form: 

^  + 1)“,  d{t  =  0)  =  di.. 

The  relative  particle  velocity  U  is  obtained  from  the  movement  equation 

Pi  ~  =  (Pi  -  P2)f/  -  U (f  =  0)  =:  Uiny 
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where  the  coefficient  of  the  particle  head  resistance  depends  on  the  Reynolds 
number  and  in  the  interval  1  <  lie  <  400  it  is  determined  by  the  formula 
Klyachko  L.A.: 

,  2d  4 
*  "  ftc  We 

A  change  of  heat,  and  mass  transfer  coefficients  {a  and  P)  in  time  is 
determined  by  dependences  Ti{i)^  d{L)  and  There  exist  critical  values  a 
and  P  characterizing  transition  from  lowtemperature  states  into  highterapera- 
ture  ones  and  vice  versa,  'rhe  steady  high-  and  lowtemperature  states  take 
place  in  the  certain  range  of  coefficients  of  heat  and  mass  transfer.  An 
operation  of  these  coefficients  is  possible  by  giving  the  initial  diameter  and 
particle  velocity.  An  investigation  of  stability  of  stationary  regimes  of  heat 
and  mas.s  transfer  {dT\/dt  ~  0)  gives  a  possibility  to  define  critical  values  of 
particle  velocities  and  diameters. 

The  dependence  of  coefficients  ct  and  P  determining  steady  and  critical 
states  of  the  particle  versus  particle  temperature  is  obtained  from  the  condition 
of  the  stationary  state  ~  Qm-r,  +  (?r  in  the  form: 

P  —  ^  _  (^1^1  T  kiQi)  {a  ±  \/A}  —  B) 

<’2P2  C2P2  {T)  -  72) 

where  A  =  O.h  (1  -  Ai  -  As) ,  B  =  Aj  A2; 

[IciQ]  +  ^2^2)  {kiQi  +  k7Q2)P2nox 

Introducing  the  complex  /  =  'ZdjNu  we  can  write  the  stationary  state 
condition  in  the  form: 

j  ^  ^  ^ _  2A2(71  -72) _ 

Aw  {kiQi  +  k2Q2) P^riox  (A  ±  VA^T. 

How  it  follows  from  this  equality,  for  7\  — >■  72  the  /  — »•  0,  and  /  — 00 
for  7i  — I-  T,(,.  Here  the  chemical  heat  secretion  power  becomes  negligible  as 
compared  with  m  the  first  case  and  with  qr  in  the  second  one. 

The  complex  /  makes  it  possible  to  take  into  account  a  change  both 
the  particle  diameter  and  the  particle  velocity.  Therefore  the  results  of 
investigations  of  dependences  d(T))  and  f/(7’))  we  can  represent  as  one 
generalized  dependence  f{Ti).  It  is  illustrated  by  figure  in  which  the  results 
of  calculations  on  the  stationary  model  (d<ished  lines)  and  the  non-stationary 
one  are  represented.  In  fig.  a  is  shown  the  behaviour  of  the  complex  2d/Nu 
which  hfis  been  calculated  on  the  non-stationary  model  taking  into  acct)unt 
the  natural  convection  and  without  it.  The  relative  gas  velocity  at  the  natural 


2d/N'iijirn  a) 


2d/Nu,fim  6) 


Fig,  Dopendences  of  the  value  /  =  2(1/ Nn  oo  stationary  particle  tempera¬ 
ture  at  Ta  =  1400  K.  a)  7’,,,  -  1400  K;  —  -  taking  Into  account,  x  x  x 
wit.hoiit  taking  into  account  of  the  naf.tiral  convection;  b)  7^  =  500  K,  - 
t'ff  ^  Q,  X  -X  X  -  Fy  ~  0.  -  -  the  calculation  iiccording  to  the 

stationary  model,  ,  x-  x  x  ,  x  x  x  -  the  calculation  axxording  to  the 
non-stationary  model,  din  —  505  //m. 


convection  w;is  obtainc(i  by  tlic  Griishof  criterion  G'r  in  the  form 


Vn.c  ^y/oJar, 
a 


Cr 


gd}^(7\~T2) 
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The  natural  convection  cnlargos  the  criterion  Nu  what  leads  to  some  decreas¬ 
ing  of  the  complex  /.  The  elTcctivc  air  velocity  during  th<?  iiatural  convection 
induces  the  increase  of  mass  transfer  what  leads  to  the  growth  of  the  burning 
rate,  consequently,  to  the  decrease  of  particle  burning  time  On  the 

induction  period  tin.d  and  the  critical  diameter  de  (the  diameter  at  which  the 
particle  extinction  takes  place)  the  natural  convection  does  not  influence,  since 
the  chemical  reactions  occur  in  the  kinetic  range  (see  Table).  The  influence 
of  the  naturd  convection  more  strongly  tells  at  large  initial  diameters  of  the 
particle. 


Table 


The  influence  of  the  natural  convection  on  the  burning  characteristics  of 
carbon  particles  at  T2  ~  7^  =  HOO  K. 


diteg ,  /rm 

taking  into  account,  (-f), 
without  taking  into  account  (- ) 
of  tlie  natural  convection 

Mndf 

df.,  p;m 

100 

- 

0.069 

0.17 

21.3 

+ 

0.067 

0.16 

21.5 

50.5 

- 

0.54 

4.74 

21.3 

0.53 

4.40 

21.5 

800 

0.99 

12.23 

21.4 

+ 

0.83 

11.58 

21.5 

in  the  fig.  b  are  represented  the  results  of  calculations  according  to  the 
non-statioiiiiry  model  for  two  cases:  i)  the  forces  of  gravitation  and  friction 
act  on  the  particle,  ii)  only  the  friction  force  acts  on  the  particle.  In  the  first 
case  {Fg  ^  0)  the  particle  ignition  takes  place  as  a  result  of  heating  of  the 
particle  by  the  convective  heat  and  mass  transfer.  In  the  second  case  {  Fg  =  0) 
heat  losses  by  radiation  at  weak  heating  by  heat  conductivity  do  not  make  it 
possible  the  particle  to  ignite. 

The  extremums  on  curves  2dlNu{Ti)  define  unsteady  particle  states: 

points  ii  and  12  correspond  to  ignition,  points  ct  and  62  .  to  extinction 

of  the  particle.  They  are  diaracterized  by  the  condition  of  the  loss  of  the 
stability  of  stationary  regimes  of  heat  and  mass  transfer  ddjdTi  —  0,  that  is 
equivalent  to  the  second  Semenov  condition 

^Qch  _  ^(gm.c  +  gr) 
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or  t.o  the  (■onl.a^’J,  of  the  curve?;  (UhiTi)  and  qm.c  + 

Branches  connecting  points  and  Cj,  i2  and  ^2  describe  the  infUience  of 
the  initial  particle  temperature  on  the  critical  partitvte  diameter  and  velocity 
which  characterize  the  ignition.  In  this  range  (it  is  shaded  in  the  figure)  the 
hysteresis  character  is  realized:  two  possible  steady  stationary  states  of  the 
particle  correspond  to  one  value  of  IdjNu  —  the  high-  and  lowtemperature 
states.  The  highteinperaturc  regime  of  heat  and  mass  transfer  is  realized 
only  in  the  case  when  the  initial  particle  temperature  at  the  fixed  value  of 
the  complex  2dfNu  exceeds  the  temperature  value  that  lies  on  the  curve 
connecting  points  i  sr  e. 

Pbr  Tiu  <  T-i  at  the  increase  of  gas  tempcrat.urc  the  critical  values  of  the 
complex  /  (and  so  for  d  and  V)  corresponding  to  ignition  and  extinction  come 
close  and  at  some  gas  temperature  value  Ts  “  the  degeneracy  of  critical 
conditions  takes  place  (hero  d^dfSTf  -  U),  that  is  for  ^2  >  the  crisis-free 
transition  from  one  stationary  state  into  another  and  vice  versa  is  realized. 
Points  i  and  e  degenerate  into  the  inflection  point  7  (71  for  and  ei,  and  72 
for  i2  and  62). 

Parameters  of  the  degeneracy  of  critical  conditions  of  heterogeneous  ignition 
and  extinction  of  the  particle  {2djNu)^^  and  {2dfNu).y.^  are  obtained  from  the 
solving  of  the  system  of  equations: 

qcU  ~  qm.r  +  qr, 

ciqch  _  +  <//) 

di\  ~  ar, 

(i'^qcb.  _  +  qr) 

The  solution  of  this  system  shows  that  in  the  range  of  large  particles  the 
value  of  {2d/Ntt)y  grows  when  the  oxidizer  concentration  increases.  The 
decrease  of  wall  temperature  enlarges  {2d}Nu)^^,  since  the  existence  of  the 
hysi.eresis  loop  in  the  range  of  large  diameters  and  small  velocities  (points  h 
and  e.2)  is  connected  with  lieal;  exchange  of  the  part.icle  by  radiation  [1].  Here 
{2dlNu)^^  practically  does  not  ciiaiige  that  is  connected  with  prevailing  of 
heat  exchange  by  heat  conductivity. 
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During  fuU-sctde  evalution  of  propulsion  systems  using  an  open  engine  test  bed 
a  considerable  amount  of  hannhil  gaseous  and  condensed  species  is  released  into 
atmosphere. 

The  calculation  of  adiabatic  temperature  and  equilibrium  composition  of  combustion 
products  given  shows  that  combustion  products  contain  considerable  amounts  of 
ecologically  harmful  substances:  HCl,  Cl,  CO.  A  series  of  calculations  was  carried  out 
for  determining  the  effect  of  cooling  on  combustion  products.  On  the  basis  of  these 
calculations  one  can  draw  conclusion  that  in  cooling  up  to  SOOisT  such  impurities  as  CO, 
Cl,  disappear  in  gaseous  products,  but  then  the  concentration  of  hydrochloride  which 
appears  in  considerable  amounts  increases.  Therefore  a  question  on  neutralisation  of 
ecologically  harmful  gases  released  in  atmosphere  in  testing  solid-fuel  charges  mainly 
comes  to  the  question  on  nentraHsation  of  SCI  or  muriatic  solution. 

Consumption  of  water  solution  of  alkali  with  respect  to  the  weight  of  dry  composition 
is  100  %.  The  results  of  calculation  show  that  neutralisation  of  products  by  water 
solution  of  sodium  hydrate  leads  to  a  sharp  decrease  of  concentration  of  hydrogen 
chloride.  At  the  same  time  sodium  chloride  appears  in  products,  it  being  in  gaseous 
state  at  temperature  lOOOJT  and  at  lower  temperatures  in  condensed  state.  In  cooling 
products  lower  1500 AT  the  content  of  hydrogen  chloride  in  gases  docs  not  depend  on 
temperature.  The  concentration  of  SCI  in  exhaust  gases  is  the  same  for  300  and  1200Ar. 
At  temperature  higher  than  ISDOAT  equilibrium  content  of  hydrogen  chloride  somewhat 
rises  and  action  of  concentration  solutions  of  alkali  it  becomes  higher  than  for  less  alkali 
solutions.  It  is  explained  by  the  fact  that  at  higher  temperatures  thermal  dissociations 
NaOH  increases  and  in  gases  at  temperatures  higher  than  IBOOlir  sodium  vaporous  are 
present  and  content  of  water  stream  increases  due  to  radicals  OS, 

Summing  the  results  of  thermodynamic  analysis  a  conclusion  may  be  drawn  that 
application  of  alkali  solutions  for  neutralisation  of  gases  forming  in  combustion  of 
solid  proi>ellant  can  reduce  to  minimum  content  of  ecologically  harmful  impurities  in 
combustion  products. 

Mathematical  modelling  is  to  help  indentify  constituent  components  of  the  gas- 
particle  cloud  at  all  the  pollution  stages  including  the  formation  of  the  impurity  clond 
in  the  ground  layer  and  its  nphft  and  spreading  in  the  atmosphere.  Each  stage,  however, 
is  characterized  by  a  time-space  scale  which  differs  greatly  from  th(»e  of  other  stages. 
Therefore,  to  improve  the  efficiency  of  computational  procedure  it  apptears  expedient  to 
do  a  8tej>-by-8tep  investigation  of  impurity  spreading.  This  approach  makes  it  pK»sibie 
to  use  output  data  for  each  stage  as  input  for  the  subsequent  st^es.  Tnis  paper  presents 
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numerical  simulations  of  the  overheated  impurity  uplift  and  propagation  in  the  atmo¬ 
sphere.  The  latter  is  a  three-dimensional  process  whose  simulation  turns  out  to  be 
complicated  and  costly.  To  avoid  these  difHcultes  we  used  a  combined  model  where 
the  vertical  impurity  dispersial  was  calculated  by  means  of  a  set  of  differential  equations, 
such  as  equations  of  motion,  turbulent  diffusion  and  heat  conduction,  while  the  hoiisontal 
dispersal  was  computed  by  a  Gaussian  model,  The  proposed  model  is  complemented  with 
assigned  semi-empirical  dependences  for  altitude  variations  of  wind  vdiocity,  turbulent 
exchange  coefficient,  atmospheric  temperature  and  turbulence  statistics.  The  initial 
conditions  for  the  impurity  uplift  calculations  (cloud  size,  temperature  and  composition) 
were  determined  from  thermodynamic  computations  of  combustion  products  in  a  solid- 
rocket  engine. 

The  overheated  exhaust  gases  appear  to  be  lighter  than  the  ambient  air,  They  move 
upwards  in  the  gravity  field  entraining  gas-particle  species  to  give  a  free  convection 
flow  in  the  form  of  a  cloud.  The  most  critical  parameters  for  the  impurity  uplift  and 
transport  processes  are  wind  velocity,  atmospheric  turbulence  intensity  and  temperature 
distribution. 

Stable  stratification  was  found  to  cause  a  faster  decrease  in  the  buoyancy  force  in 
the  upward  direction  than  in  the  isothermal  atmosphere  and  hence  a  decrease  in  the 
uplift  velocity.  In  the  external  part  of  the  flow  there  forms  a  temperature  defect  which 
tends  to  increase  with  altitude  and  with  the  increase  in  the  stratification  parameter. 
The  buoyancy  force  eventually  goes  to  zero  and  then  iissumes  a  negative  value  which 
gives  rise  to  a  maximum  of  the  uplift  height  and  the  reversed  movement  of  the  cloud. 
The  cloud  does  not  rise  higher  than  a  certain  penetration  altitude.  The  trajectory  of  the 
center-of-graviiy  of  the  cloud  varies  smoothly,  first  ascending,  then  descending  under  the 
action  of  the  vertical  density  gradient. 

The  absence  of  wind  (calm)  coupled  with  the  temperature  inversion  is  the  most 
unfavorable  factor  for  the  estimation  of  the  atmospheric  pollution.  If  the  height  of  the 
boundary  of  the  inversion  layer  is  fairly  low,  then  for  light  winds  the  impurity  cloud 
seems  to  hover  over  the  pollution  source  and  the  impurity  concentration  increases 
considerably.  For  example,  as  the  wind  velocity  decreases  from  1  m/s  to  zero,  the 
maximum  impurity  concentration  near  the  Earth’s  surface  may  increase  approximately 
by  an  order  of  magnitude. 

Our  calculations  show  that  above  30  km  nitrogen  oxides  generally  thought  to  be  the 
major  ozone  destructors  are  virtually  missing  from  the  exhaust  products  and  ozone  may 
decay  basically  via  the  following  chlorine  cycle, 

For  the  solid  rocket  propellant  formulation  in  question  the  content  of  nitrogen  oxides 
above  30  km  is  much  lower  than  that  of  chlorine- bearing  compounds.  Therefore,  at  z  >  30 
km  the  major  mechanism  of  ozone  depletion  is  associated  with  chlorine  cycle.  However, 
in  the  tropopause  and  lower  stratosphere  the  impact  of  nitrogen  oiddes  found  in  exhaust 
products  is  fairly  significant  and  neglect  of  nitrogen  contribution  to  ozone  depletion  may 
introduce  large  errors.  Characteristic  distributions  of 

As  a  result  of  the  interaction  of  the  SPC  products  with  atmospheric  ozone  there 
develops  a  region  in  the  stratosphere  wherefrom  O3  is  either  missing  or  its  concentration 
is  negligible  as  compared  to  natural  ozone  (what  is  known  as  the  ozone  hole).  As 
time  elapses,  the  hole  starts  moving  relative  to  the  Earth’s  surface  together  with  the 
atmospheric  motion  of  masses,  and  its  size  and  shape  undergo  essential  changes. 
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On  the  basis  of  the  analysis  of  the  time  history  we  recognize  two  stages  in  the 
ozone  hole  life  time.  In  the  first  stage  there  occurs  mixing  and  chemical  interaction 
of  exhaust  products  with  atmospheric  ozone  to  destroy  the  ozonoephere,  causing  the 
area  of  the  ozone  hole  to  expand.  In  the  second  stage,  when  ozone  destructors  from 
exhaust  products  have  reacted  with  ozone  masses  surrounding  the  trace  combustion 
species  and  chemical  reactiorft  terminated,  the  gas  content  in  the  trace  exhaust  products 
is  determined  solely  by  diffusion  processes  which  act  to  supply  ozone  from  fairly  distant 
regions  to  the  neighborhood  of  the  ozone  hole.  As  a  result  the  area  of  the  ozone  hole  is 
decreased  and  the  original  atmospheric  condition  regained. 

The  area  of  the  ozone  hole  b  at  its  maximum  (approximately  5  *  10^  )  within  4 

hours  after  the  launch.  Notably,  due  to  advective  processes  the  center  of  the  ozone  hole 
b  shifted  with  the  wind  by  40-50  km  from  the  launch  pad.  Thus  the  40-50  km  area  in 
the  vicinity  of  the  launch  pad  b  a  potentially  hazardous  zone  because  it  b  in  thb  area 
that  most  of  the  ozone  destruction  occurs  and  ultraviolet  and  hard  cosmic  radiation  b 
enhanced. 

The  calculations  done  allowed  us  to  estimate  the  pollution  level  for  different 
meteorologies  conditions  and  pollution  sources  and  provide  recommendations  as  to  how 
to  reduce  its  dbastrous  effect  on  the  environment. 
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